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theoretically  and  experimentally  as  a  function  of  primary  frequency  directivity. 
The  broadband  parametric  array  describes  the  sound  field  produced  by  the  non¬ 
linear  interaction  of  the  frequency  components  of  a  pulsed  primary  signal. 

In  the  present  investigation,  the  secondary  sound  field  of  the  baseband  modula¬ 
tion  products  is  considered.  The  primary  transmissions  are  produced  by 
physically  realizable  sources,  and  the  broadband  parametric  array  for  nonplane 
wave  fields  is  developed.  ' 
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A  general  solution  for  the  parametric  array  for  nonplane  wave  fields  is  presented 
in  the  form  of  the  superposition  of  primary  field  components.  These  components 
are  assumed  to  be  spherically  spreading  from  the  origin  and  characterized 
throughout  the  sound  field  by  a  single,  frequency  dependent  directivity  function, 
or  off-axis  transfer  function.  The  general  solution  is  reduced  from  a  four¬ 
dimensional  integration  to  a  two-dimensional  integration  by  assuming  the 
directivity  to  be  symmetric  about  the  acoustic  axis  and  restricting  the  region  of 
validity  to  the  farfield  of  the  parametric  array.  The  farfield  solution  is 
obtained  for  a  Gaussian  pulse  primary  transmission  emitted  by  a  circular  piston 
source.  The  solution  is  evaluated  in  limiting  forms  as  a  function  of 
directivity.  When  the  primary  beam  is  extremely  narrow,  the  parametric  array  is 
absorption  limited  and  the  parametric  array  for  plane  wave  fields  is  recovered. 
The  characteristics  of  the  secondary  off-axis  transfer  function  are  determined  by 
the  absorption  at  the  primary  frequencies  and  are  monotonically  decreasing  with 
angle  off-axis  and  secondary  frequency.  When  the  Rutherford  beamwidth  is  small 
compared  to  the  primary  beamwidchs,  the  parametric  array  is  diffraction  limited. 
The  characteristics  of  the  secondary  off-axis  transfer  function  in  the 
diffraction  limit  generally  follow  that  of  the  primary  off-axis  transfer  function 
with  off-axis  angle  but  are  relatively  insensitive  to  secondary  frequency.  As 
the  parametric  array  shifts  from  absorption  limited  to  the  diffraction  limited, 
the  dependence  on  secondary  frequency  changes  from  quadratic  to  linear.  This 
change  is  characterized  by  a  shift  toward  lower  frequencies  in  the  secondary 
spectra  and  a  reduction  in  the  on-axis  frequency  at  the  spectral  peak  by  a 
factor  of  'll  for  a  Gaussian  pulse  primary  transmission. 

■  An  experimental  study  of  the  broadband  parametric  array  is  presented  which 
provides  data  with  variations  in  experimental  parameters  and  an  emphasis  on  the 
behavior  as  a  function  of  primary  directivity.  The  results  for  low  amplitude 
primary  signals  generally  confirm  the  theoretical  predictions.  An  extreme 
sensitivity  to  the  sidelobe  structure  in  the  primary  off-axis  transfer  function 
is  noted  in  the  frequency  at  the  spectral  peak  of  the  secondary  power  spectra. 

No  such  sensitivity  was  noted  in  the  level  at  the  spectral  peak  or  in  the 
secondary  power  beam  pattern.  For  high  amplitude  primary  signals,  the  results 
indicate  that  envelope  distortion  due  to  finite  amplitude  effects  in  the 
primary  signal  causes  spectral  distortion  in  the  secondary  signals.  The  spectral 
distortion  is  greatest  on-axis  where  the  frequency  at  the  spectral  peak  is 
reduced.  Far  off-axis,  the  characteristics  are  the  same  as  those  for  low  ampli¬ 
tude  primary  transmissions. 
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yields  a  peak  secondary  transmission 
magnitude  of  r 

vector  from  the  origin  to  the  field  point 
magnitude  of  rQ 

vector  from  the  origin  to  the  volume  element  of 

integration 

magnitude  of  R 

vector  from  the  volume  element  of  integration  to 

the  field  point 

nearfield  parameter  =  So/A 

entropy 

crossectional  area 

simplified  time  function 

primary  and  secondary  peak  source  level, 

respectively 

step-down  ratio 

primary  and  secondary  peak  sound  pressure  level, 

respectively 

time 

reference  time 
duration  of  a  pulse 

magnitude  of  the  sound  velocity 
xv 
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u 

Uo,Ul,U2 

sound  velocity 

static,  first  order  (primary),  and  second  order 

u(t) 

(secondary)  sound  velocity,  respectively 

gate  time  function 

V 

o 

interaction  volume 

X 

rectilinear  coordinate 

y 

rectilinear  coordinate 

z 

rectilinear  coordinate 

a 

absorption  coefficient 

a. 

(a^n+a^m)/2,  average  absorption  coefficient  for 

alc 

the  primary  wave 

absorption  coefficient  at  the  carrier  or  primary 

aln*a2n 

center  frequency 

absorption  coefficient  for  the  nth  frequency 

ad 

components  of  the  primary  and  secondary  waves, 

respectively 

absorption  coefficient  corresponding  to  the 
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difference  frequency 

1  +  B/2A 

Y 

-V 

angle  between  rQ  and  r 

6(0 

Dirac  delta  function 

A 

strength  of  a  delta  function 

e 

angle  between  the  acoustic  axis  and  the  projection 

* 

of  the  acoustic  axis  in  the  x,z  plane 

angle  of  rotation  about  the  acoustic  axis  for  the 

projection  of  r,  angle  of  rotation  about  r  (Muir), 

■+ 

angle  between  the  projection  of  rQ  on  the  x-y  plane 
and  r  (Blue) 


*o 
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projection  of  r^ 

for  the 
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angle  between  and  r  (Muir) ,  angle  between  the 
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o 

(Blue) 

angular  increment 

angle  between  the  projection  of  r  on  the  x-y  plane 
and  the  x  axis  (Blue) 

2irf,  angular  frequency 

angular  frequency  of  the  primary  and  secondary 
waves,  respectively 

the  nth  angular  frequency  components  of  the  primary 
and  secondary  waves ,  respectively 
angular  frequency  of  the  carrier  or  primary  center 
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static,  first  and  second  order  density,  respectively 

angle  between  r  and  the  acoustic  axis 
o 

angle  measured  with  respect  to  the  acoustic  axis 
angle  between  the  acoustic  axis  and  r 

o 

half  power  angle 
half  power  beamwidth 

primary  and  secondary  half  power  beamwidth, 
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total  distance,  r^+R;  angle  between  r  and  the 
projection  of  the  acoustic  axis  in  the  x-z  plane 
differential  operator 
Laplacian  operator 
d’Alembertian  operator 

complex  conjugate  as  a  superscript,  convolution 
between  two  functions 
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CHAPTER  I 
INTRODUCTION 

The  parametric  acoustic  arr  w,  as  orig'.nally  described  by 
Westervelt,  has  been  the  subject  of  numerous  theoretical  and  experi¬ 
mental  studies.  In  the  most  general  sense,  the  parametric  array  is  the 
sound  field  of  the  interaction  products  that  result  from  the  mixing  of 
two  or  more  frequency  components  in  a  common  volume  of  a  medium.  The 
sound  field  of  the  baseband  interaction  products  or  difference  fre¬ 
quencies  has  been  of  greatest  interest  since  these  components  survive 
after  the  primary  components  from  which  they  were  generated  have  been  - 
severely  attenuated  because  of  the  increasing  absorption  of  energy  in 
the  medium  at  higher  frequencies.  The  characteristics  of  the  secondary 
sound  field  are  unconventional  by  comparison  with  linear  sources  in 
that  highly  directive  low  frequency  sound  with  very  low  sidelobe  radia¬ 
tion  may  be  transmitted  with  physically  small  apertures. 

The  most  widely  studied  parametric  array  is  the  absorption 
limited  narrowband  parametric  array  for  plane  wave  fields,  which  was  the 
subject  of  Westervelt 's^  study.  This  parametric  array  requires  a 
physically  unrealizable  primary  source  configuration  for  implementation 
and  therefore  may  only  be  approximated  in  practice.  This  array  is 
termed  absorption  limited  since  the  off-axis  properties  are  determined 
by  the  absorption  at  the  primary  frequencies  in  the  medium.  The 
narrowband  parametric  array  for  nonplane  wave  fields,  which  allows  for 
realizable  primary  source  configurations,  has  been  addressed  in 
additional  treatments.  This  type  of  parametric  array  is  difficult  to 
express  in  closed  form  except  for  a  few  limiting  cases.  This  array  may 
be  dominantly  absorption  limited  in  some  configurations;  but  in  most 
cases,  the  off-axis  characteristics  are  influenced  by  the  off-axis 
properties  of  the  primary  source  and  the  parametric  array  is  termed 
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diffraction  limited.  The  broadband  parametric  array  is  described  by  the 
secondary  field  produced  by  a  pulsed  primary  transmission  where  a 
continuum  of  primary  frequencies  interact  to  produce  a  continuum  of 
secondary  frequencies.  The  broadband  parametric  array  for  plane  wave 
fields,  which  is  the  absorption  limited  broadband  parametric  array,  has 
been  described  in  a  few  investigations  which  will  be  discussed.  The 
broadband  parametric  array  for  nonplane  wave  fields  results  from  physi¬ 
cally  realizable  primary  source  configurations  and  may  be  dominantly 
absorption  or  diffraction  limited.  The  broadband  parametric  array  for 
nonplane  wave  fields  is  the  physically  realizable  broadband  parametric 
array  and  is  the  principal  subject  of  this  dissertaiton. 

In  the  second  chapter,  a  theoretical  review  is  presented  that 
describes  the  derivation  of  an  inhomogeneous  wave  equation  character¬ 
izing  the  secondary  sound  field  for  an  unbounded  medium.  The  various 
types  of  parametric  arrays  are  discussed  which  include  narrowband  and 
broadband  parametric  arrays  for  both  plane  wave  and  nonplane  wave 
fields.  Included  in  this  discussion  are  references  and  descriptions  of 
many  of  the  principal  treatments  which  contribute  to  the  understanding 
of  parametric  arrays. 

Detailed  theoretical  considerations  of  broadband  parametric 
arrays  are  addressed  in  the  third  chapter.  These  considerations  include 
a  detailed  study  of  the  broadband  parametric  array  for  plane  wave 
fields.  Off-axis  characteristics  are  related  to  on-axis  characteristics 
through  a  linear  systems  representation  of  an  off-axis  impulse  response 
and  transfer  function.  The  properties  of  the  secondary  field  are 
evaluated  in  both  the  time  and  frequency  domains.  A  source  strength 
density  is  derived  that  is  generally  applicable  for  the  broadband  para¬ 
metric  array  for  plane  wave  fields.  Solutions  for  the  secondary 
pressure  derived  from  this  source  strength  density  are  expressed  in 
general  and  for  a  Gaussian  pulse  primary  transmission  from  a  circular 
piston  source.  The  limiting  forms  of  these  expressions  are  also  con¬ 
sidered. 
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The  results  of  a  comprehensive  experimental  evaluation  of  the 
broadband  parametric  array  for  nonplane  wave  fields  are  presented  in 
Chapter  IV.  Both  the  on-axis  and  off-axis  properties  of  broadband 
parametric  arrays  are  examined  for  low  amplitude  primary  waves  trans¬ 
mitted  with  a  Gaussian  pulse  from  circular  pistons  of  different  radii. 
Data  are  obtained  for  directivity  ratios  (ratio  of  the  primary  half¬ 
power  beamwidth  to  the  Rutherford  beamwidth  from  0.54  to  3.69.  This 
span  of  parameter  values  transition  from  dominantly  absorption  limited 
behavior  to  diffraction  limited  behavior.  These  data  are  compared  with 
theoretical  results  in  the  farfield  of  the  parametric  array  and  with 
absorption  and  diffraction  limited  forms  of  the  solution.  The  perform¬ 
ance  of  the  broadband  parametric  array  is  also  evaluated  experimentally 
for  high  amplitude  primary  waves. 

The  summary  and  conclusions  regarding  these  investigations  of 
the  broadband  parametric  array  are  considered  in  Chapter  V.  Several 
appendices  are  included  that  contain  supportive  theoretical  and 
experimental  details. 


CHAPTER  II 


THEORETICAL  BACKGROUND 

The  physical  basis  of  the  parametric  acoustic  array  (and  of 
nonlinear  acoustics  in  general)  is  the  variability  of  the  speed  of  an 
acoustic  wave.  This  variability  results  from  the  nonlinearity  of  the 
medium  as  described  by  an  equation  of  state  and  Trom  the  convection  of 
a  sound  wave  with  the  speed  of  the  part: *  -:.a  it  sets  in  motion. 

As  sound  waves  propagate,  amplitude  dep*?- 'rects  occur  which 
include  harmonic  and  intermodulation  disr,-” .  in  the  frequency  domain 
and  waveform  distortion  in  the  time  Aotr^y.-  '^Tie  parametric  acoustic 
array  is  a  special  case  of  the  interm-odulation  distortion  effect  which 
accompanies  the  propagation  of  an  acoustic  signal  in  a  nonlinear  medium. 

The  theoretical  development  of  the  parametric  acoustic  array 
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was  developed  from  the  work  of  Lighthill  ’  who  obtained  an  exact 
equation  describing  the  radiation  of  sound  produced  by  turbulence. 
Westervelt  *  used  Lighthill’ s  results  in  studying  the  scattering  of 
sound  by  sound  and  in  the  original  development  of  the  parametric  acous¬ 
tic  array. ^  In  this  development,  Lighthill1 s  equation  was  reduced  to 
an  inhomogeneous  wave  equation  describing  the  secondary  pressure  field 
that  is  generated  by  the  interaction  of  acoustic  signals. 

In  the  next  section  a  derivation  of  the  inhomogeneous  wave 
equation  is  given  which  begins  with  the  equations  of  fluid  mechanics 
and  an  isentropic  equation  of  state.  Following  this  derivation, 
specific  types  of  parametric  arrays  are  reviewed. 

A.  The  Inhomogeneous  Wave  Equation 

The  equations  which  describe  the  continuity  of  mass,  the 
conservation  of  momentum,  and  the  state  of  a  fluid  are  the  fundamental 
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relations  governing  acoustics.  For  a  lossless  fluid  with  no  real 
sources  present,  the  continuity  and  momentum  equations  in  Eulerian 
coordinates  may  be  written  as  follows: 

mass  continuity  equation 


Is-  +  V 

at 


(pu)  =  0 


and 


(2.1) 


momentum  balance  equation 


■—  (pu)  +  Vp  +  pu(V*u)  +  (u*V)pu  =  0 


(2.2) 


where 


p  is  total  density, 
t  is  time, 

7  is  a  spatial  differential  operator, 
u  is  particle  velocity,  and 
p  is  total  pressure. 


For  a  lossless  fluid,  the  entropy  remains  constant  and  the 
fluid  remains  at  equilibrium  during  the  passage  of  a  sound  wave.  Under 
these  conditions  the  state  of  the  fluid  may  be  described  by  an  isen- 
tropic  equation  of  state  relating  pressure  and  density  given  by 


P  =  P(P) 


S=constant 


(2.3) 


where  S  is  entropy. 

This  equation  may  be  expanded  in  a  Taylor  series  about  the 
ambient  or  static  density  of  the  fluid.  When  the  first  three  terms  of 
this  expansion  are  retained,  the  isentropic  equation  of  state  may  be 
expressed  as 
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'  w 


p 


p  +  c  (p-p  )  + 
o  o  o 


^(5)K)J 


(2.4) 


where 

Pq  is  the  ambient  or  static  pressure, 
is  the  ambient  or  static  density, 
c^  is  the  small  signal  sound  speed,  which  is  defined  by  the 
relation 


and 


2 

c 

o 


(If) 

P=P  ,  S=constant 
o 


(2.5) 


1 

A 


S=constant 


(2.6) 


The  quantity  B/A  is  known  as  the  parameter  of  nonlinearity 
and  must  be  determined  experimentally  for  specific  fluids, 
mental  measurements  have  indicated  a  value  for  B/A  of  5.2 
degassed  fresh  water  at  30°C. 


for  fluids 
Experi- 
for 


The  acoustic  field  may  be  described  more  concisely  through 
a  combination  of  the  mass  continuity  equation  (2.1),  the  momentum 
balance  equation  (2.2),  and  the  isentropic  equation  of  state  (2.4). 

This  is  done  using  a  perturbation  technique  in  a  manner  similar  to 

7  8  1  4  5  9. 

Muir,  Truchard,  and  Westervelt.  *  ’  ’  This  technique  consists 

of  replacing  the  dependent  variables  appearing  in  Eqs.  (2.1),  (2.2), 

and  (2.4)  by  the  sum  of  their  equilibrium  and  higher  order  variational 

components,  and  then  forming  separate  equations  which  must  be  satisfied 

by  the  variables  of  each  order.  If  only  static,  first,  and  second 

order  relations  are  of  interest,  the  field  variables  may  be  written 
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u  =  u  +  U-  +  u„ 
ol  2 

p  =  p0  +  Pi  +  p2  (2>7) 

P  =  po  +  P1  +  p2  • 

It  is  assumed  that  the  fluid  is  at  rest,  and  therefore  uq=0.  When 
Eq.  (2.7)  is  substituted  into  Eqs.  (2.1),  (2.2),  and  (2.4),  the  first 
and  second  order  acoustic  relations  are  found.  The  first  order  rela¬ 
tions  are  those  of  linear  acoustics  and  describe  what  is  referred  to 
as  the  primary  field  in  the  context  of  the  parametric  acoustic  array. 
These  relations  are  given  by  the  following: 

mass  continuity  equation 


(2.8) 


momentum  balance  equation 


3U1  . 

D  '  *  '  + 

Mo  at 


isentropic  equation  of  state 


(2.9) 


P 


1 


2 

Co  P1 


(2.10) 


wave  equation  for  p^ 


2  1  3*  _2 

5  pl  -  1  ~2  P1  "  D  Pi  '  0 
c  at 
o 


where  V  is  the  Laplacian  operator,  and 


(2.11) 
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is  the  d 'Alembertian  or  wave  operator. 


The  second  order  relations  are  extensively  studied  in 
nonlinear  acoustics  and  describe  the  secondary  field  in  the  context  of 
the  parametric  acoustic  array.  These  relations  are  given  as  follows: 

mass  continuity  equation 

3p2  ^ 

—  +  PqV  •  u2  +  V  •  (p^)  =  0  ,  (2.12) 


momentum  balance  equation 


—  (P0VpA>  +  ^P2  +  P0^1(^*u1>  +  (u^o^  =  0  ,  (2.13) 


isentroplc  equation  of  state 


P2  CoP2  '  2A  2 


(2.14) 


c  p 
o  o 


wave  equation  for  p2 


□  P- 


_1 _  B_ 

4  2A 

P  c 


)+(ux 


(2.15) 


The  first  order  wave  equation  (2.11)  is  homogeneous  because 
of  the  initial  assumption  of  a  source  free,  lossless  medium.  However, 
the  second  order  wave  equation  (2.15)  is  inhomogeneous.  The  source 
term  in  Eq.  (2.15)  exists  because  of  the  presence  of  the  primary  field 
and  represents  a  virtual  source.  Further  simplification  of  this  source 
term  is  necessary  to  allow  a  detailed  analysis  of  the  secondary 
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pressure  field.  With  the  assumption  that  u  is  irrotational , 
Eq.  (2.15)  may  be  written 


P  c 
o  o 


(2.16) 


When  the  mass  continuity  equation  (2.8),  the  momentum  balance 
equation  (2.9),  the  isentropic  equation  of  state  (2.10),  and  the  wave 
equation  for  particle  velocity  in  the  form  of  Eq.  (2.11)  are  combined, 
Eq.  (2.16)  becomes 


1 _  B_ 

4  2A 


P 


c 

o  o 


(2.17) 


-  P, 


2  4 
P  c 

l_o  o 


IZj 

at 


P  C 
o  o 


55T*T“!| 


2  2 
D  c 
o  o 


(vPl-vPl)| 


Simplification  of  Eq.  (2.17)  is  accomplished  with  identities  which  are 
obtained  by  utilizing  the  momentum  balance  equation  (2.9)  and  the  wave 
equation  for  pressure  (2.11).  These  identities  are  given  by 


When  Eqs.  (2.18)  and  (2.19)  are  substituted  into  Eq.  (2.17),  the  result 
is 
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(2.20) 


□  P2  = 


4 

P  c 
o  o 


/(' 

V1  +  2A  / 

3  ■ 

(pl)‘ 


3t 


S  (3//Pldt)  +  2°  l“ll 

_  o  o 


- 2  pl 

2P  c2  , 

o  O  J 
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Westervelt  ’  ’  interpreted  the  d 'Alembertian  operation  in 
the  source  term  of  Eq.  (2.20)  as  a  local  contribution  to  the  secondary 
pressure  field.  This  local  pressure  does  not  exist  external  to  the 
region  determined  by  the  primary  pressure  field  and  gives  no  contribu¬ 
tion  to  the  radiated  secondary  pressure  field.  Muir2  analyzed  the 
magnitude  of  the  local  secondary  pressure.  His  results  indicate  that 
for  plane  harmonic  primary  waves  there  is  no  local  secondary  pressure, 
and  for  spherical  harmonic  waves  the  local  secondary  pressure  vanishes 
asymptotically  with  range.  After  the  local  secondary  pressure  terms 
are  dropped,  Eq.  (2.20)  becomes 
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□  P2  = 


4 

P  c 
o  o 


2(pl)J 


3t 


(2.21) 


where  8=1+(B/2A).  Westervelt  derived  this  equation  initially  by 

assuming  the  small  signal  plane  wave  impedance  relation.  He  re- 
18 

derived  this  equation  in  the  manner  which  has  been  presented,  thereby 
extending  the  validity  of  the  original  result  to  arbitrary  wave  con¬ 
figurations.  Equation  (2.21)  may  be  rewritten  in  the  form  adopted  by 
Westervelt, 


_  2 

□  P2  "  -P 


1SL 
o  3t 


(2.22) 


where  q=(B/poCo)(9(pl)2  /3tj  is  the  source  strength  density.  This 
inhomogeneous  wave  equation  (2.22)  has  the  same  form  as  the  small  signal 
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wave  equation  for  a  medium  in  which  there  are  real  mass  sources  of 
source  strength  density  q.  The  source  strength  density  in  Eq.  (2.22) 
gives  rise  to  the  generation  of  secondary  acoustic  waves  which  result 
from  the  nonlinear  interaction  of  primary  acoustic  waves.  The  proper¬ 
ties  of  the  source  strength  density  are  those  associated  with  a  virtual 
source  that  is  distributed  in  space  over  the  region  occupied  by  the 
primary  waves  and  delayed  in  time  by  the  propagation  speed  in  the 
medium.  The  secondary  waves  are  sometimes  referred  to  as  the  scattered 
waves,  by  analogy  with  scattering  theory  where  each  elemental  volume 
or  scatterer  in  the  interaction  volume  acts  as  a  secondary  sound  source. 
The  general  solution^  of  Eq.  (2.22)  for  a  distributed  source  of 
source  strength  density  q  is  given  in  the  time  domain  by 


where 


tq  is  a  vector  from  the  origin  to  the  volume  element  of 
integration, 

r  is  a  vector  from  the  origin  to  the  field  point, 

R  =  r-r  , 

$  , 


R  -  I r-r0 1 » 

v  is  the  interaction  volume, 
o 

t  is  an  initial  reference  time, 
o 

6(*)  is  the  Dirac  delta  function,  and 

[•]  is  the  time  dependent  Green's  function  for  a  three  dimen¬ 
sional  unbounded  medium. 
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If  the  source  strength  density  is  a  periodic  function  of  time,  then  the 
solution  of  Eq.  (2.22)  in  the  time  domain  may  also  be  written 


P2(r»t> 


a  Iff 


t) 


-jk2nR 


dv  ,  (2.24) 
o 


where 

q^(r^,t)  is  the  nth  harmonic  component  of  the  source  strength 
density, 

k„  =  /c  ,  and 
2n  2n  o 

is  the  angular  frequency  of  the  nth  component  of  the  secondary 
zn  ‘ 

signal. 

The  general  solution  of  Eq.  (2.22)  may  also  be  expressed  in  the 
frequency  domain  using  Fourier  transform  relationships  given  by 


and 


g(t)  =  _1{G(f)  > 


oo 

/  G (f )  ej2lTft  df 

—  CO 


(2.25) 


G(f )  =  fF  { g  (t) }  = 


/ 

—  CO 


g(t) 


e-j2TTft 


dt 


(2.26) 


where 

{ • }  is  the  Fourier  transform, 

(F  {•}  is  the  inverse  Fourier  transform, 
f  =  ^  ,  and 

J  =  ^  • 

The  general  solution  of  Eq.  (2.22)  in  the  frequency  domain  is  written 

dv  ,  (2.27) 
o 
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P2(r,f)  =  P0 


Iff  T{fr«<vv} 


r  -^2h-| 

a 

4ttR  " 


l 


where  [•]  is  the  frequency  domain  Green's  function  f or  a  three- 
dimensional  unbounded  medium. 


The  assumptions  under  which  Eqs .  (2.23)  and  (2.27)  were 
derived  include  the  assumption  of  lossless  fluid  medium  with  no  real 
sources.  This  assumption  may  now  be  relaxed  without  seriously  com¬ 
promising  the  validity  of  the  solutions  for  the  secondary  field  com¬ 
ponents.  The  existence  of  a  primary  field  which  gives  rise  to  a 
secondary  field  through  the  source  density  function  implies  the  presence 
of  a  real  source  which  generates  the  primary  field.  This  real  source 
can  be  included  through  the  specification  of  both  the  spatial  and  the 
signal  characteristics  of  the  primary  radiation.  Losses  may  be  included 
in  an  ad  hoc  manner  by  noting  the  average  attenuation  over  frequency 
bands  or  the  attenuation  of  each  frequency  component  in  the  primary 
and  secondary  radiation  and  applying  the  principle  of  superposition. 

The  solutions  of  Eqs.  (2.23)  and  (2.27)  have  been  studied  for  a  variety 
of  source  strength  density  functions  which  are  derived  through  varia¬ 
tions  in  the  spatial  configuration  of  the  interaction  volume  and  in  the 
primary  acoustic  signal.  It  is  the  differences  in  characteristics  of 
the  source  strength  density  function  which  give  rise  to  different  types 
of  parametric  arrays.  Several  kinds  of  parametric  arrays  have  been 
considered  and  include  both  narrowband  and  broadband  transmissions  for 
plane  and  nonplane  wave  fields. 

B.  The  Narrowband  Parametric  Array  for  Plane  Wave  Fields 

The  narrowband  parametric  array  for  plane  wave  fields  is  the 
subject  of  the  original  study  by  Westervelt^  and  of  numerous  subsequent 
developments.  Westervelt  assumed  a  small  amplitude  primary  signal 
which  is  the  linear  combination  of  two  monochromatic  plane  waves.  The 
first  order  pressure  is  given  by 

-ar 

p.  (r  ,0,t)  =  P  e  0  [cos  (id.  ,  t— k. ,  r  )  +  cos  (ui, -t-k,  »r  )].  (2.28) 
1  o  O  11  II  O  iz  iz  o 
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where 


Pq  is  the  amplitude  of  each  wave  at  the  origin, 

a  =  (a^-t-a^)^  is  the  average  absorption  coefficient  for  the 

primary  waves,  and 

r  =  |r  | 
o  o 

The  primary  waves  are  collimated  along  the  x  axis  in  a  region  of 
circular  cross  section  Sq,  and  are  attenuated  as  a  function  of  range 
by  absorption.  The  geometry  for  this  configuration  is  shown  in 
Fig.  2.1. 


The  source  strength  density  may  be  found  by  substituting 
Eq.  (2.28)  into  Eq.  (2.22).  Through  the  dependence  on  the  square  of 
the  primary  pressure,  the  second  harmonic  at  each  primary  frequency  and 
modulation  products  given  by  the  sum  and  difference  of  each  primary 
frequency  are  generated  in  the  source  strength  density.  The  absorption 
coefficient  varies  generally  as  the  square  of  the  frequency;  therefore, 
the  difference  frequency  component  will  experience  less  rapid  attenua¬ 
tion  as  a  function  of  range  than  the  other  modulation  components.  With 
the  retention  of  only  the  difference  frequency  component,  the  source 
strength  density  may  be  expressed  in  complex  form  as 


q(ro,0,t) 


j  CD  .  6P 

J  d  o 

2  4 
P  c 
o  o 


-2-o  ^VVo5 


(2.29) 


where 


Re{q}  or  Im{q}  must  be  taken  to  give  physical  meaning,  and 


wd  = 
kd  = 


I“l2-“ll> 
I  k12_kll I 


The  source  strength  density  is  periodic;  therefore,  the 
solution  in  the  time  domain  for  the  pressure  at  the  difference  fre¬ 
quency  is  given  by  Eq.  (2.24)  and  written 
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FIGURE  2.1 

GEOMETRY  FOR  ANALYSIS  OF  A  PARAMETRIC 
ACOUSTIC  ARRAY  FOR  PLANE  WAVE  FIELDS 


It  is  assumed  that  the  field  point  P  is  far  from  the  zone  of  inter¬ 
action,  or  r>>r  .  With  this  assumption,  the  Fraunhofer,  or  farfield 
o 

approximation  may  be  invoked,  which  amounts  to  representing  R  by  the 
law  of  cosines  and  performing  a  binomial  expansion,  giving 

R  =  r  -  r  cos9  +  ...  ,  (2.31) 

o 

where 

r  =  | r |  ,  and 

0  is  the  angle  between  r  and  the  acoustic  axis. 

In  Eq.  (2.30),  the  slowly  varying  part  of  the  Green's  function  is 
replaced  by  the  first  term  of  the  expansion  and,  in  the  rapidly  varying 
part  or  phase  of  the  Green's  function,  R  is  replaced  by  the  first  two 
terms  of  the  expansion.  The  interaction  volume  has  simple  cylindrical 
geometry  which  physically  may  be  attained  only  approximately.  The 
volume  integration  in  Eq.  (2.30)  reduces  to  a  one-dimensional  integral 
along  the  x  axis  given  by 


Pd(r,0,t) 


2  2 
wTgP  S 
d  o  o 

Atto  c^r 
o  o 


j(udt— kdr) 


exp [-j (k ,-2ja-k ,cos0) r  ]dr 
q  a  o  o 


(2.32) 


This  integration  must  be  carried  out  to  a  distance  £,  such  that  the 
source  strength  density  is  attenuated  to  a  negligible  value  and  the 
generation  of  pressure  at  the  difference  frequency  is  effectively 
completed.  The  effective  length  of  the  parametric  array,  given  by 
r=l/(2a),  is  the  distance  at  which  the  amplitude  of  the  source  strength 
density  drops  to  1/e  of  its  original  value.  Operationally,  the  upper 


17 


limit  of  the  integral  may  be  extended  to  infinity,  which  yields  the 
solution  for  the  pressure  at  the  difference  frequency  in  the  farfield 
of  the  parametric  array  as 


Pd(r,0,t) 


“d6PoSo  JV-V’ 

- -  e 

8ttd  c  r 
o  o 


-“dr 

e 


(2.33) 


where  attenuation  as  a  function  of  range  due  to  vis  :ous  absorption  at 
the  difference  frequency  has  been  included  ad  hoc. 


The  amplitude  directivity  factor  for  the  secondary  pressure 
from  Eq.  (2.33)  is  given  by 


Pd(r,0,t) 

pd(r,0,t) 


(2.34) 


which  yields  a  half-power  beamwidth  of 


0 


HP 


4  sin 


-1 


(2.35) 


The  directivity,  given  by  Eq.  (2.34),  is  referred  to  as  the  Rutherford 
directivity  because  of  its  similarity  in  form  to  Rutherford  scattering 
in  atomic  physics. 


The  narrowband  parametric  array  for  plane  wave  fields  is  an 
end-fire  volumetric  array  of  virtual  sources  which  i.s  exponentially 
shaded  in  length  by  absorption.  As  a  result  of  this  shading,  the 
directional  response  for  the  parametric  array  is  absorption  limited 
and  is  practically  devoid  of  sidelobe  radiation.  In  addition,  the 
dependence  of  the  half-power  beamwidth  on  a  and  kd  allows  highly  direc¬ 
tive  radiation  at  relatively  low  frequencies  from  physically  small 

apertures.  Experimental  investigations  to  verify  Westervelt's  theoreti- 

11  12 

cal  results  were  performed  by  Bellen  and  Beyer,  Berktay  and  Smith, 
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Hobaek,  ^  Zverev  and  Kalachev,^  and  Muir  and  Blue.^  Westervelt ' s^ 
analysis  was  generally  confirmed  in  these  investigations;  however, 
several  of  the  experiments  indicated  greater  directivity  at  the  differ¬ 
ence  frequency  than  predicted.  To  account  for  this  increase  in 

16 

directivity,  Westervelt’s  treatment  was  extended  by  Naze  and  Tj«5tta, 

17  12 

Berktay,  and  Berktay  and  Smith  to  allow  for  an  aperture  factor. 

This  factor  is  simply  the  directivity  function  at  the  difference  fre¬ 
quency  that  is  associated  with  the  finite  cross  section  of  the  inter¬ 
action  volume.  The  aperture  factor  is  included  as  a  multiplier  to  the 
Rutherford  directivity  and  can  dominate  the  directivity  of  the  paramet¬ 
ric  array  if  the  cross  section  of  the  interaction  volume  is  large  in 
terms  of  wavelengths  at  the  difference  frequency. 

There  are  several  limitations  associated  with  the  analytical 
development  of  the  narrowbeam  parametric  array  for  plane  wave  fields. 

It  is  assumed  that  the  primary  signal  consists  of  plane  waves  which  are 
collimated  in  a  narrow  beam  of  uniform  cross-sectional  area.  These 
spatial  characteristics  are  unrealistic  in  practice  and  may  only  be 
approximated  by  an  appropriate  selection  of  parameters  such  that  all 
the  interaction  takes  place  in  the  nearfield,  i.e.,  the  Fresnel  region, 
of  a  primary  source.  Another  limitation  is  that  narrowband  development 
does  not  account  for  the  initiation  and  the  termination  of  the  primary 
transmission,  which  must  occur  in  all  physical  situations.  Finally,  the 
analytical  results  are  valid  only  far  from  the  interaction  region,  i.e., 
in  the  farfield  of  the  parametric  array.  Some  of  the  limitations  which 
have  been  discussed  are  overcome  in  the  theoretical  models  which  are 
presented  in  the  next  section. 

C.  The  Narrowband  Parametric  Array  for  Nonplane  Wave  Fields 

The  narrowband  parametric  array  for  nonplane  wave  fields 
includes  arrays  with  realizable  spatial  distributions  of  the  primary 
waves.  These  parametric  arrays  differ  from  those  for  plane  wave  fields 
in  that  much  of  the  interaction  takes  place  in  the  farfield  or 
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Fraunhofer  region  of  the  primary  source  and  the  parametric  array  is 

diffraction  limited.  One  of  the  earliest  investigations  of  nonplane 

18 

wave  fields  was  by  Lauvstad,  Naze,  and  Tj^tta  who  obtained  an  approxi¬ 
mate  solution  for  the  radiation  at  the  difference  frequency  in  the  far- 

field  of  the  parametric  array  produced  by  diverging  primary  sound  beams 

19 

with  directivity  associated  with  a  circular  piston.  Lauvstad  and 
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Berktay  '  *  also  derived  farfield  solutions  for  both  cylindrical 

22 

and  spherical  primary  waves.  Mellen  and  Moffett  noted  that  with 
bounded  aperture  primary  sources,  the  radiation  at  the  difference  fre¬ 
quency  resulted  from  plane  wave  interaction  in  the  nearfield  of  the 
source  and  from  nonplane  wave  interaction  in  the  farfield  of  the  source. 

Matched  plane  wave  and  nonplane  wave  farfield  solutions  for  the  para- 
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metric  array  were  studied  by  Moffett  and  Mellen  for  circular  apertures 

24 

and  by  Moffett,  Mellen,  and  Konrad  for  rectangular  apertures.  Berktay 
25 

and  Leahy  examined  the  farfield  performance  of  the  parametric  array 
for  primary  waves  emitted  from  circular  and  rectangular  apertures  by 
computing  correction  factors  to  Westervelt's  plane  wave  solution. 

Investigations  restricted  to  the  nearfield  of  the  parametric  array  were 

26  27 

conducted  by  Rolleigh,  Novikov,  Rudenko,  and  Soluyan,  and  Novikov, 

28 

Rybachek,  and  Timoshenko.  In  each  of  these  studies,  small  signal 

absorption  was  ignored  in  order  to  obtain  theoretical  results.  Other 

treatments  of  the  narrowband  parametric  array  for  nonplane  wave  fields 

include  an  investigation  of  the  parametric  array  in  air  by  Bennett  and 
29 

Blackstock  and  several  studies  based  on  approximate  differential 
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equations  by  Fenlon,  Fenlon  and  McKendrie,  and  Tj«5tta  and  Tjdtta. 

The  experimental  data  which  were  presented  with  many  of  the  theoretical 
developments  establish  the  range  of  validity  for  each  theoretical  treat¬ 
ment  and  the  effect  of  various  approximations. 

Closed  form  solutions  for  the  pressure  at  the  difference 
frequency  for  nonplane  wave  fields  are  generally  not  tractable  when  a 
bounded  aperture  is  considered  as  the  primary  source.  Muir  analyzed 
the  narrowband  parametric  array  for  spherical  wave  fields.  In  his 

analysis,  no  approximations  were  made  regarding  position  vectors  and  the 
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scattering  integral  was  evaluated 
amplitude  primary  signal  which  is 
chromatic  spherical  waves  emitted 
baffle.  The  pri:  .ry  signal  which 


numerically.  ie  assumed  a  low 
a  linear  combination  of  two  mono- 
from  a  circular  piston  in  an  infinite 
he  assumed  has  the  form 


R  -nc.r 

P,  (r  ,0,t)  =  P  ,  — -  D(k  a  sino)  e  °  cos (w  t-k, , r  ) 
1  o  ol  r,  il  11  11  o 


Ro  _ai9r 

+  P  0  —  D(k,„a  sino)  e  “  cos (w, -t-k, „r  ) 
oZ  r  1Z  1Z  1Z  o 

o 


(2.36) 


where 

P  ,  ~is  the  amplitude  of  each  wave  at  R  . 
ol,o2  r  o 

Rq  is  a  near field  parameter, 

a  is  the  radius  of  the  circular  piston, 

a  is  the  angle  between  rQ  and  the  acoustic  axis,  and 

D(*)  is  the  directivity  function  for  a  circular  piston  in  an 

infinite  rigid  baffle. 

The  primary  waves  described  by  Eq .  (2.36)  propagate  as  spherical  waves 
and  are  attenuated  as  a  function  of  range  by  viscous  absorption.  The 
interaction  volume  is  determined  by  the  common  region  of  the  farfield 
beam  patterns  at  each  primary  frequency.  The  geometry  for  this  con¬ 
figuration  is  shown  in  Fig.  2.2.  Using  the  form  given  in  Eq.  (2.36) 
for  the  primary  waves,  the  source  strength  density  may  be  computed 
using  Eq.  (2.22).  Retaining  only  the  difference  frequency  component, 
the  resulting  source  strength  density  is  a  periodic  function  of  time; 
and  the  solution  in  the  time  domain  for  the  pressure  at  the  difference 
frequency  is  given  by  Eq.  (2.24).  With  the  inclusion  of  attenuation  as 
a  function  of  range  due  to  viscous  absorption  at  the  difference  fre¬ 
quency,  the  solution  for  the  secondary  pressure  is  written 
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FIGURE  2.2 

GEOMETRY  FOR  ANALYSIS  OF  A  PARAMETRIC 
ACOUSTIC  ARRAY  FOR  NONPLANE  WAVE  FIELDS 


(Allat  Muit  and  WiHatta) 


ARL  UT 
AS-80-966-P 
JVM . GA 
3.25-80 


9  9 


Pd(r,e,t) 


r  <P  +0  77 

ol^Po2RoV9  f  f  f  "-l 

oCokllk12a  J  J  J 


[k^a  sino]  sino] 

—  o 

sin  o 
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x  —  exp  -  {[jk  +  (a ,,+a  )  ]  r  +  [jk  -‘■a  ]R}  sin^d^d^dr 

K  a  il  iZ  o  a  a 


(2.37) 


where 


2  ?  2  1/2 

sino  =  { [sin(0-<(i)  +  sin<t>  cos0 (1-cosi+O  ]  +  sin"^  sin  , 

/  2  2  \ly/2 
R  =  (r  +ro-2rrocos<Ji  j  ,  and 

<f>e  is  an  additional  angular  increment  needed  to  adequately  account 
for  the  main  lobe  of  the  primary  beam  patterns. 
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Muir  and  Willette  derived  the  expression  given  in  Eq.  (2.37) 
and  evaluated  it  numerically.  Experimental  results  were  presented  which 
support  agreement  with  theory.  They  also  studied  the  effects  of  various 
farfield  approximations  for  the  position  vector  by  comparison  with  the 
exact  position  vector.  Equation  (2.37)  completely  accounts  for  the 
effects  of  diffraction  in  the  farfield  over  the  main  lobe  of  the  pri¬ 
mary  radiation  and  is  valid  in  the  farfield  of  the  primary  waves  and 
throughout  the  nearfield  and  the  farfield  of  the  parametric  array. 


A  solution  similar  to  that  of  Muir  and  Willette  was  derived 
43 

by  Blue  utilizing  the  geometry  shown  in  Fig.  2.3.  He  assumed  an 
equation  in  the  form  of  Eq.  (2.36)  for  the  primary  waves  and  of 
Eq.  (2.24)  for  the  secondary  pressure.  Blue  also  made  the  usual 
farfield  approximation  and  wrote  the  solution  as 
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X 


d11(<J'»4>)D-L2  )  cosi^d4)d^ 


\  (a^-t-o^)  +  jk-dsin2  jy  cos  1  [cos (V-i^) cos  (4>-<}>)  ]  j 


where  D(.)  is  a  primary  directivity  function.  This  solution  provides 
good  insight  into  the  directivity  of  the  parametric  array  in  general. 
Equation  (2.38)  is  a  double  convolution  integral  of  the  cosine  weighted 
product  of  the  primary  beam  functions  with  the  Rutherford  directivity. 
This  realization  implies  that  the  directivity  of  the  parametric  array 
follows  the  Rutherford  directivity  when  one  or  both  of  the  primary  beam 
functions  can  be  treated  as  spatial  impulse  functions.  Conversely, 
when  the  Rutherford  directivity  acts  as  a  spatial  impulse  function,  the 
product  of  the  beam  functions  at  the  primary  frequencies  determines  the 

directivity  of  the  parametric  array.  Experimental  results  have  been 
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presented  by  Smith,  Bennett  and  Slack,  and  Huckabay  to  confirm 

this  behavior. 


The  aperture  factor,  or  the  directivity  function  at  the 
difference  frequency  associated  with  the  effective  cross  section  of  the 
interaction  volume,  can  dominate  the  directivity  of  the  parametric  array 
for  plane  wave  fields.  The  aperture  factor  associated  with  the  para¬ 
metric  array  for  nonplane  wave  fields  was  studied  by  Muir,^  Moffett 
23  38  39 

and  Mellen,  ’  and  Hobaek  and  Vestrheim.  The  results  of  these 
studies  indicate  that  the  aperture  factor  may  be  included  as  a  simple 
multiplier  of  the  convolution  integral  in  Eq.  (2.38).  If  the  convolu¬ 
tion  integral  is  dominated  by  the  Rutherford  directivity,  then  the  aper¬ 
ture  factor  can  be  important;  however,  if  the  primary  product  direc¬ 
tivity  dominates  the  convolution  integral,  then  the  aperture  factor 
becomes  negligible. 


D.  The  Broadband  Parametric  Array  for  Plane  Wave  Fields 

The  broadband  parametric  array  for  plane  wave  fields  was 
first  studied  by  Berktay^2  who,  like  Westervelt,  assumed  a  primary 
signal  collimated  in  a  narrowbeam  of  cross-sectional  area  Sq,  as  shown 
in  Fig.  2.1.  In  contrast  to  the  narrowband  parametric  array  for  plane 
wave  fields,  Berktay  assumed  a  small  amplitude  primary  signal  in  which 
a  single  frequency  carrier  wave  is  modulated  by  an  envelope  function. 
This  signal  is  written 


pl(ro,0,t)  =  Po  6 


-a  r 
lc  o 


fit  -  —  cos  (a),  t-k..  r  ) 
\  cl  lc  lc  o 


.  (2.39) 


It  is  assumed  that  the  envelope  function  f(t-r  /c  )  varies  slowly 

o  o 

compared  to  cos (u^t-k^r^)  t  and  that  the  absorption  coefficient  at  the 
carrier  frequency  may  be  used  to  represent  the  absorption  of  the  band- 
limited  signal.  Only  the  lowest  frequency  components  are  retained  from 
the  source  density  function,  which  is  found  by  substituting  Eq.  (2.39) 
into  Eq.  (2.22).  The  source  density  function  is  not  a  periodic  function 
of  time;  therefore,  the  solution  in  the  time  domain  for  the  secondary 
pressure  is  found  from  Eq.  (2.23).  After  introducing  the  farfield 
approximation  given  by  Eq.  (2.31),  the  solution  is  written 
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which,  on  axis  (8=0),  reduces  to 
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Again,  the  upper  limit  of  integration  may  be  extended  to  r  =“  since 

viscous  absorption  at  the  carrier  frequency  will  diminish  the  integrand 
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and  effectively  truncate  the  secondary  generation.  Moffett  derived 
the  result  given  in  Eq.  (2.40)  and  compared  it  with  experimental 
data^’^  for  various  envelope  functions.  Additional  theoretical  and 
experimental  comparisons  were  reported  by  Moffett,  Westervelt,  and 
Beyer,  ^>^3  Merklinger,  Mellen  and  Browning,  ^  and  Davy  and 

Hixson. 


Berktay  analyzed  the  frequency  response  of  the  broadband 
parametric  array  for  plane  wave  fields  by  substituting  the  source 
strength  density  into  Eq.  (2.27),  which  gives  the  solution  for  the 
secondary  pressure  in  the  frequency  domain.  After  introducing  the  far- 
field  approximation,  Eq.  (2.31),  the  solution  is  written 


P2(r,0,f) 
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Westervelt  offered  a  description  of  the  on- axis  and  off-axxs 
behavior  in  the  time  domain  of  the  broadband  parametric  array  for  plane 
wave  fields.  He  returned  to  the  solution  for  the  narrowband  parametric 
array  for  plane  wave  fields,  Eq.  (2.33),  and  noted  the  primary  pressure 
dependence  given  by 


Pd(r,9,t)  « 


.  2 
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on-axis,  this  dependence  may  be  written 


and  off-axis,  where  k,>>a,  Che  dependence  may  be  written 
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Westervelt  generalized  the  results  of  Eq.  (2.44)  and  (2.45)  and  noted 
that  the  on-axis  secondary  pressure  of  the  broadband  parametric  array 
for  plane  wave  fields  should  be  proportional  to  the  second  time  deriva¬ 
tive  of  the  envelope  function,  and  the  off-axis  secondary  pressure 
should  be  proportional  to  the  first  time  derivative  of  the  envelope 
function. 

The  broadband  nature  of  the  parametric  array  for  plane  wave 
fields  was  analyzed  by  Muir,^’^  and  Muir  and  Blue^  by  applying  the 
principle  of  superposition  to  the  steady  state  solution  obtained  by 
Westervelt.'*'  In  the  frequency  domain,  the  solution  for  the  secondary 
pressure  at  a  discrete  difference  frequency  may  be  written 


F*(f)  F(f+f .)  df  ,  (2.46) 

a 

where 

*  denotes  complex  conjugate, 

A(r,e,fd)  <$(f-fd)  =  JF  j[pd(r,e,t))/p2  j  , 

^d  =  ’  anc* 

F(f)  =  fF { f ( t ) >  ,  the  amplitude  spectrum  of  the  primary  pressure 
waveform. 

The  validity  of  Eq.  (2.46)  is  subject  to  the  same  assumptions  under 
which  Westervelt’s  steady  state  solution  was  derived.  These  assumptions 
require  the  frequencies  in  the  primary  band  to  be  much  greater  than  the 
difference  frequency;  therefore,  the  envelope  function  given  by 
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P2^r,0’fd)  = 


A(r,9,fd) 


/ 
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Eq.  (2.39)  must  be  slowly  varying  with  respect  to  the  carrier  frequency. 
The  autocorrelation  theorem  may  be  used  to  transform  the  integral  in 
Eq.  (2.46)  to  give 

P2(r,0,fd)  =  A(r,0,fd)  6(f-fd)  fF { f2(t) >  ,  (2.47) 

2 

where  only  the  slowly  varying  component  of  f  (t)  is  retained.  The 
secondary  pressure  may  be  expressed  as  a  continuous  function  of  fre¬ 
quency  through  the  superposition  of  each  difference  frequency  component 
as  given  in  the  frequency  domain  by 

P2(r,8,f)  =  y* p2(r,0,fd)  dfd 

=  A(r,6,f)  £F{f2(t) }  .  (2.48) 

The  waveform,  or  the  solution  for  the  secondary  pressure  in 
the  time  domain,  is  obtained  from  the  inverse  Fourier  transform  of 
Eq.  (2.48),  which  is  written 

P2(r,e,t)  =  3r_1{A(r,0,f ) }  *  f2(t)  ,  (2.49) 

where  (F  ^{A(r,0,f)}  is  the  impulse  response  of  the  parametric  array  for 
plane  wave  fields  and  *  between  two  functions  denotes  convolution.  The 
solutions  obtained  by  superposition,  Eqs.  (2.48)  and  (2.49),  are 
straightforward  and  in  forms  which  are  suitable  for  computer  analysis. 

E.  The  Broadband  Parametric  Array  for  Nonplane  Wave  Fields 

The  broadband  parametric  array  for  nonplane  wave  fields  is  an 
extension  of  Berktay's12  broadband  parametric  array  development  where 
most  of  the  interaction  occurs  in  the  farfield  or  Fraunhofer  region  of 
the  primary  source  radiation.  The  primary  signal  has  features  asso¬ 
ciated  with  both  the  narrowband  array  for  nonplane  wave  fields. 
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Eq.  (2.36),  and  the  broadband  array  for  plane  wave  fields,  Eq.  (2.39). 

A  small  amplitude  primary  signal  which  conforms  to  the  geometry  in 
Fig.  2.2  is  given  by 

R  -a  (r  -R  )  /  r  \ 

p^(rQ,t)  =  ~  Pq  e  °  °  D(o)  fft  -  ~  I  cos(w^t-k1ro)  .  (2.50) 

This  form  for  the  primary  signal  has  been  used  in  previous 
studies’*^’ and  is  applicable  when  it  is  assumed  that  the  primary 
signal  is  sufficiently  narrowband  that  parameters  characterising  the 
primary  signal  are  constant  over  the  bandwidth  of  the  transmission. 
Under  this  assumption,  the  directivity  function  simply  determines  the 
amplitude  dependence  of  the  primary  signal  as  a  function  of  angle  off- 
axis.  The  source  density  function  is  found  by  substituting  Eq.  (2.50) 
into  Eq.  (2.22)  and  retaining  only  the  slowly  varying  term. 

The  resulting  source  density  function  is  not  a  periodic 
function  of  time;  therefore,  the  solution  for  the  secondary  pressure  in 
the  time  domain  is  given  by  Eq.  (2.23).  This  solution  is  written 


P2(r,0,t) 


„_2_2 

BP  R 
o  o 

8irp  c2 
o  o 


2a,  R 
1  o 


fffi < 


-2a, r  „ 

1  °  ~2,  , 
e  D  (o) 


(2.51) 


In  the  frequency  domain,  the  solution  is  found  from  Eq.  (2.27)  and  is 
given  by 
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2  2  2 

gP  R  201^ 
o  o  2  lo 
P  (r,0,f)  =  - 7-  e 

8irp  c 
o  o 


fffi? 


-2a, r  0 
1  o  -2,  . 
e  D  (a) 


(2.52) 


The  similarities  between  the  forms  of  the  narrowband  and 
broadband  parametric  arrays  for  nonplane  wave  fields  are  evident.  Much 
of  the  detailed  theoretical  and  experimental  work  has  concentrated  on 
narrowband  parametric  arrays  and  the  broadband  parametric  array  for 
plane  wave  fields.  There  have  been  few  treatments  devoted  to  the 
broadband  parametric  array  for  nonplane  wave  fields. 

Moffett  and  Mello  presented  a  model  in  which  the  secondary 
pressure  in  the  frequency  domain  was  found  by  superposition  of  solutions 
for  the  narrowband  parametric  array  for  nonplane  wave  fields.  They 
assumed  only  amplitude  dependent  directivity  functions  and  restricted 
the  details  of  the  analysis  to  the  on-axis  behavior  of  the  broadband 
array  for  plane  wave  fields.  This  treatment  is  similar  to  that  of  Muir 
who  obtained  the  broadband  parametric  array  for  plane  wave  fields  by 
superposition  of  Westervelt's  result  over  the  primary  frequency 
components . 

51  52 

Rolleigh  and  Mize  et  al.  considered  the  broadband 
parametric  array  for  nonplane  wave  fields  as  a  linear  system.  A 
primary  signal  similar  to  that  of  Eq.  (2.50)  and  a  solution  in  the  form 
of  Eq.  (2.51)  were  assumed.  Spherical  geometry  was  assumed,  a  change 
of  variables  was  made,  and  Eq.  (2.51)  was  put  in  the  form  of  a  spatial 
convolution  integral  given  by 
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(2.53) 


o 


where 

C  =  rQ  +  R, 

S(t-C/cQ)  involves  the  double  time  derivative  of  the  square  of  the 
envelope  function,  and 

h(£,<f>,4<)  is  the  spatial  impulse  response  for  the  parametric 
array  with  real  primary  beam  functions. 


The  spatial  impulse  response  h(£,<}!,\j))  involves  a  double  integration 
over  angles.  Rolleigh  reduced  this  double  integration  to  a  single 

integration  with  the  assumption  of  a  Gaussian  beam  function  for  the 
primary  signal.  The  solution  in  the  frequency  domain  for  the  secondary 
pressure  is  given  by 


P2(r,9,f) 


H(r,fl,f) 


(2.54) 


where 

H(r,6,f)  =  JF  {h(£,4>,ij/) }  ,  and 

h(£»<J>>i|')  =  0  for  C  <  R  . 

o 

Rolleigh  showed  that  JF{h(C,4>»'l') }  was  amenable  to  numerical  analysis  due 
to  the  well  behaved  properties  of  the  integrand,  and  that  on  the  acous¬ 
tic  axis  JF  {h(?,<j>,ij>) }  could  be  evaluated  in  closed  form.  The  off- 
axis  behavior  was  examined  in  the  farfield  of  the  parametric  array  for 
angles  near  the  acoustic  axis.  The  results  of  this  analysis  indicated 
that  spherical  spreading  and  diffraction  change  the  off-axis  spectral 
distortion  relative  to  that  for  the  broadband  parametric  array  for 
plane  wave  fields. 
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F. 


The  Parametric  Acoustic  Array  with  Finite  Amplitude  Primary  Waves 

The  various  parametric  acoustic  arrays  which  have  been 
discussed  were  derived  for  small  amplitude  primary  signals.  For  large 
amplitude  primary  waves,  finite  amplitude  effects  must  be  considered 
in  order  to  adequately  describe  the  behavior  of  a  parametric  array. 

Finite  amplitude  effects  are  the  result  of  variations  in 
sound  speed  with  the  particle  velocity  in  the  direction  of  sound  propa¬ 
gation.  In  a  fluid,  these  variations  are  described  by 

c’  =  c  +  u  +  |~  u  ,  (2.55) 

o  2A 


where 

c1  is  the  propagation  speed  of  any  particular  point  on  the  acoustic 
wave, 

— V  -V 

u  =  u*n  is  the  particle  velocity  in  the  direction  of  sound 
propagation,  and 

n  is  a  unit  vector  in  the  direction  of  sound  propagation. 
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The  varying  propagation  speed,  as  discussed  by  Blackstock, 
is  a  function  of  the  small  signal  sound  speed  in  the  medium  and  the 
particle  velocity.  The  second  term  in  Eq.  (2.55)  results  from  con¬ 
vection  and  the  third  term  results  from  the  nonlinearity  of  the 
pressure-density  relation.  According  to  Eq.  (2.55),  the  portions  of  an 
acoustic  wave  in  a  compression  phase  and  in  an  expansion  phase  propa¬ 
gate  faster  and  slower,  respectively,  than  the  small  signal  sound  speed. 
This  behavior  results  in  cumulative  waveform  distortion  with  propagation 
distance  and  the  eventual  formation  of  a  shock  wave.  Larger  amplitude 
waves  form  shock  waves  at  shorter  distances  than  lower  amplitude  waves. 
For  very  weak  waves,  small  signal  dissipation  is  dominant  over  cumula¬ 
tive  waveform  distortion  and  the  formation  of  a  shock  wave  is  precluded. 
The  cumulative  effects  of  waveform  distortion  with  propagation  distance 


have  been  studied  by  Blackstock  *-n  t^ie  1 requency  domain.  He 

obtained  a  general  solution  in  the  form  of  a  Fourier  series  for  the 
distortion  with  propagation  distance  for  an  originally  sinusoidal 
pressure  wave.  The  fundamental  component  of  this  series  decreases 
monotonically  with  propagation  distance.  Each  harmonic  component  first 
increases,  and  then  decreases,  with  range.  This  behavior  suggests  that 
as  an  originally  sinusoidal  wave  propagates  with  finite  amplitude, 
energy  is  converted  from  the  fundamental  to  the  harmonic  components 
during  the  formation  of  a  shock  wave.  An  increase  in  the  level  of  the 
original  wave  causes  an  increase  in  the  harmonic  component  growth  rate. 
Eventually,  the  amplitude  of  the  fundamental  component  saturates  at  an 
upper  limit  which  cannot  be  increased  with  a  further  increase  in  level 
of  the  original  wave. 

The  parametric  acoustic  array  with  finite  amplitude  primary 
levels  has  been  the  subject  of  many  theoretical  and  experimental 
studies.  As  large  amplitude  primary  waves  propagate  coaxially,  wave¬ 
form  distortion  and  the  growth  of  harmonic  components  will  occur.  Each 
frequency  component  will  interact  with  each  other  component  to  produce 
modulation  products,  as  indicated  by  Brinkmann. The  higher  order 
interactions  are  very  weak;  therefore,  in  the  context  of  the  parametric 
array,  the  secondary  signal  or  difference  frequency  component  is  the 
most  important.  The  primary  effect  which  results  from  finite  amplitude 
primary  levels  is  that  the  conversion  of  energy  from  the  fundamental  to 
harmonic  components  acts  as  a  range  dependent  attenuation,  in  addition 
to  viscous  absorption.  The  effects  of  finite  amplitude  attenuation  on 

the  parametric  array  have  been  studied  by  Bartram  and  Westervelt , ^ 
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Bartram,  and  Zverev,  Kalachev,  and  Stepanov.  The  ratio  of  initial 

power  at  the  fundamental  to  the  final  power  of  the  secondary  signal,  or 

the  conversion  efficiency  of  the  parametric  array,  is  lower  for  finite 

amplitude  primary  waves  than  for  small  signal  primary  waves.  As 

the  primary  level  of  the  parametric  array  is  increased,  the  dependence 

of  the  secondary  level  on  the  primary  level  changes  from  quadratic  to 

linear  because  of  the  onset  of  finite  amplitude  effects.  The 
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performance  of  the  parametric  array  for  saturated  primary  signals  has 
been  evaluated  by  Mellen  and  Browning, ^  Mellen,  Browning,  and  Konrad, ^ 
and  Willette  and  Moffett. ^  The  length  of  the  parametric  array  becomes 
shorter  as  the  total  attenuation  is  increased  with  the  inclusion  of 
finite  amplitude  attenuation.  With  parametric  arrays  for  plane  wave 
fields,  the  increase  in  attenuation  results  in  an  increase  in  the 
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half-power  beamwidth,  given  by  Eq.  (2.35),  as  shown  by  Merklinger. 

In  considering  parametric  arrays  for  nonplane  wave  fields,  the  amplitude 
response  of  the  primary  beam  functions  varies  such  that  the  effects  of 
finite  amplitude  attenuation  are  greatest  on  the  acoustic  axis.  This 
results  in  a  blunting  of  the  secondary  radiation  pattern  relative  to 
the  small  signal  case  with  an  increase  in  the  half-power  beamwi'dth. 

The  effects  which  finite  amplitude  primary  waves  have  on  parametric 

arrays  for  nonplane  wave  fields  are  discussed  by  Muir,^  Moffett  and 

22  23  24 

Mellen,  ’  and  Moffett,  Mellen,  and  Konrad.  The  increase  in  half¬ 
power  beamwidth  for  the  parametric  array  is  a  direct  result  of  the 
truncation  of  the  virtual  end-fire  parametric  array  at  shorter  ranges 
due  to  finite  amplitude  attenuation.  Additional  treatments  of  the 
parametric  array  with  finite  amplitude  primary  waves  are  given  by 

v  i  62»63 
Fenlon. 
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CHAPTER  III 

THEORETICAL  CONSIDERATIONS  FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

A.  Radiation  in  the  Farfield  of  the  Broadband  Parametric  Array  for 
Plane  Wave  Fields  with  Low  Amplitude  Primary  Waves 

The  absorption  limited  broadband  parametric  array  for  plane 
wave  fields  was  reviewed  in  Chapter  II.  Several  on-axis  and  off-axis 
properties  were  described.  In  this  section,  the  on-axis  and  off-axis 
analysis  is  extended  with  the  assumption  of  a  particular  primary  signal 
envelope  function.  The  properties  of  the  secondary  transmission  are 
evaluated  in  terms  of  the  primary  transmission  and  expressed  in  both  the 
time  and  frequency  domains.  This  study  is  preparatory  to  the  develop¬ 
ment  of  the  broadband  parametric  array  for  nonplane  wave  fields  which 
reduces  to  the  broadband  parametric  array  for  plane  wave  fields  in  the 
absorption  limit. 

1.  General  Solution 


The  general  solutions  for  the  pressure  field  of  the  broadband 
parametric  array  for  plane  wave  fields  are  summarized  in  Eqs.  (2.39), 
(2.40),  (2.41),  and  (2.42).  An  analysis  of  the  spectral  distortion  of 
a  broadband  signal  which  is  transmitted  by  the  parametric  array  for 
plane  wave  fields  as  a  function  of  off-axis  angle  may  be  considered  by 
treating  the  on-axis  and  the  off-axis  properties  of  Eq.  (2.42)  as  the 
respective  input  and  output  of  a  linear  system.  This  relationship  is 
given  by 


P2(r,0,f)  =  P2(r,0,f)  H20(f) 


(3.1) 


where  HOQ  is  the  off-axis  transfer  function  for  the  second  order  field. 
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From  Eq.  (2.42)  the  off-axis  transfer  function  for  ;he  broadband 
parametric  array  for  plane  wave  fields  is  given  by 


H2e(f)  - 

which  is  closely  related  to  the  amplitude  directivity  factor  for  the 
parametric  array  for  plane  wave  fields  given  by  Eq.  (2.34).  The  off- 
axis  power  transfer  function  I H2q (^ ) I  shown  in  Fig.  3.1.  Continuing 
the  linear  system  approach,  the  solution  for  the  secondary  pressure 
field  in  the  time  domain  is  expressed  by 

P2(r,0,t)  =  p2(r,0,t)  *  h20(t)  ,  (3.3) 

where  h20(t)  is  the  off-axis  impulse  response  function,  which  is  written 

h20(t)  =  T'1  [h20 (f)j  .  (3.4) 

The  linear  system  approach  separates  the  radiation 
characteristics  of  the  parametric  array  into  on-axis  and  off-axis  com¬ 
ponents.  Analyses  of  off-axis  spectral  distortion  of  specified  or 
measured  signals  are  developed  from  the  on-axis  signal  characteristics 
to  yield  off-axis  solutions  in  both  the  time  and  frequency  domains. 

2.  On-Axis  Radiation 


1  +  3 


sin 


2/6 


lc 


(3.2) 


The  on-axis  pressure  in  the  farfield  of  the  broadband 
parametric  array  for  plane  wave  fields  is  given  in  the  time  domain  by 
Eq.  (2.41).  In  the  frequency  domain,  the  on-axis  solution  is  given  by 
Eq.  (2.42)  evaluated  at  8=0,  and  expressed  by 
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P2(r,0,f) 


6P2S 
o  o 


If, ip  C  <*- 
o  o  lc 


-jk2r  2  , 

<**2  5f(t) 


n.  1-. 


The  analysis  of  specific  signal  properties  associated  with  th 
broadband  parametric  array  for  plane  wave  fields  requires  the  specif ica 
tion  of  a  pulse  envelope  function  for  the  primary  transmission.  The 
assumption  of  a  Gaussian  pulse  envelope  is  mathematically  tractable  and 
is  a  reasonable  approximation  to  pulse  envelopes  which  may  be  attained 
in  practice.  The  particular  characteristics  of  a  Gaussian  pulse  trans¬ 
mission  are  detailed  in  Appendix  A. 


A  Gaussian  envelope  function  is  assumed  t)  have  the  form 


f(t) 


2  2 

-n  t 


(A. I?) 


which  has  a  Fourier  transform  given  by 


ff|f(t)|  =  F(f ) 


(A. 13) 


where  n  is  a  constant  which  determines  the  effective  duration  of  the 
pulse  and  therefore  the  effective  width  of  the  frequency  spectrum.  The 
spectral  characteristics  of  the  Gaussian  pulse  transmission  may  be 
specified  by  defining  the  constant  n  in  the  frequency  domain  given  by 


n  =  2.668 


'lc 


(A. 18) 


where  is  the  quality  factor  for  the  primary  transmission  determined 
by  the  center  frequency  divided  by  the  half-power  bandwidth.  If  the 
pulse  duration  is  defined  as  the  inverse  of  the  half -power  band¬ 
width,  then  the  pulse  duration  is  the  interval  in  time  between  the 
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0.17  amplitude  levels  of  the  pulse  envelope.  Using  Eq .  (A. 12)  in 
Eq.  (2.41),  the  on-axis  solution  for  the  farfield  secondary  pressure 
with  a  Gaussian  pulse  transmission  is  given  in  the  time  domain  by 


P2(r,0,t) 


2 

BP  S 
o  o 


16ttp  c  ra, 
o  o  lc 


4n  4n 


1 


x  exp 


-2n 


(3.6) 


This  result  was  derived  by  Berkay;  the  shape  of  the  pressure  pulse  is 
shown  in  Fig.  3.2.  The  on-axis  solution  in  the  frequency  domain  for  the 
farfield  secondary  pressure  with  a  Gaussian  pulse  transmission  is  found 
from  Eqs.  (A. 13)  and  (3.5)  and  is  given  by 


P2(r,0,f) 


BP^S  -jk?r 

o  o  ^ 


16ttp  c  ret. 
o  o  lc 


2  1  /if 
u»2  exp 


(3.7) 


For  a  constant  value  of  n,  the  frequency  at  which  the  on-axis  secondary 
pressure  is  a  maximum  is  found  by  maximizing  Eq.  (3.7)  with  respect  to 
f?.  Differentiating  Eq.  (3.7)  and  setting  the  result  to  zero  gives 


i  C-)‘ 


f2+  1 


f2  exp 


0 


(3.8) 


This  expression  indicates  a  maximum  when 


f 


2p 


n 


(3.9) 


This  result  was  derived  by  Findeisen . ^  Using  the  definitions  of  n 
given  in  Eqs.  (A. 16)  and  (A. 18),  the  peak  frequency  for  the  on-axis 
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FIGURE  3.2 

ON-AXIS  PULSE  SHAPE  FOR  THE  BROADBAND 
PARAMETRIC  ARRAY  FOR  PLANE  WAVE  FIELDS 
WITH  A  GAUSSIAN  PULSE  PRIMARY  TRANSMISSION 
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secondary  transmission  can  be  related  to  the  half-power  bandwidth  and 
Q^.  The  frequency  at  which  the  on-axis  secondary  pressure  is  a  maximum 
may  also  be  expressed  by 


f  =  i  2  If  -f  ! 

2p  • r 11  121 


=  1.2 


'lc 


(3.10) 


where  Fj.  and  F^  are  the  lower  and  upper  half-power  frequencies, 
respectively,  for  the  primary  transmission.  The  step-down  ratio  (SDR), 
the  ratio  of  primary  center  frequency  to  secondary  frequency  at  which 
the  on-axis  secondary  pressure  is  a  maximum,  is  defined  by  Eq.  (3.10) 
and  is  written 


SDR 


(3.11) 


As  Q1  diminishes,  the  SDR  becomes  small  and  the  frequency  at  which  the 
on-axis  secondary  pressure  is  a  maximum  approaches  the  primary  center 
frequency.  This  situation  violates  the  assumption  tha  the  Gaussian 
pulse  envelope  function  varies  slowly  compared  to  cos(w^ct)  and  repre¬ 
sents  an  unrealistic  case  for  the  broadband  parametric  array.  If  the 
SDR  is  restricted  to  values  of  2  or  greater,  then  is  restricted  to 
values  of  2.4  or  greater.  The  frequency  at  which  the  on-axis  secondary 
pressure  is  a  maximum  as  a  function  of  the  primary  center  frequency  is 
shown  in  Fig.  3.3. 

A  representation*^  of  the  power  spectral  density  function  for 
finite  duration  signals  may  be  derived  from  the  energy  spectral  density 
function  by  weighting  the  energy  spectral  density  function  with  respect 
to  the  duration  of  the  signal,  or 
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PRIMARY  CENTER  FREQUENCY,  -  kHz 
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FIGURE  3.3 

FREQUENCY  AT  WHICH  THE  ON-AXIS  SECONDARY  PRESSURE  IS  A  MAXIMUM 
AS  A  FUNCTION  OF  THE  PRIMARY  CENTER  FREQUENCY  WITH  THE  Q  OF 
A  GAUSSIAN  PULSE  PRIMARY  TRANSMISSION  AS  A  PARAMETER 
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where 

P  (f)  is  the  power  spectral  density  function, 

E  is  the  energy  spectral  density  function,  and 

T  is  the  duration  of  the  signal. 

The  normalized  power  spectral  density  function  for  the  on-axis  trans¬ 
mission  of  the  broadband  parametric  array  for  plane  wave  fields  with  a 
Gaussian  pulse  primary  transmission  is  derived  from  Eq.  (3.7)  and  given 
by 


P2xxN(f) 


Po(r,0,f)  Po(r,0,f) 
P2(r,0,f2  )  P2(r,0,f  ) 


(3.13) 


If  the  power  spectral  density  function  for  the  on-axis  secondary  signal 
is  normalized  to  the  on-axis  signal  produced  by  a  Gaussian  pulse  primary 
transmission  with  a  quality  factor  of  Q^p>  then  Eq.  (3.13)  becomes 


P-  vt  ( f ) 
2xxN 


✓2(2.668) 


exp 


(3.14) 


The  power  spectral  density  function  for  the  on-axis  secondary  signal  is 
shown  in  Fig.  3.4  for  several  values  of  of  the  Gaussian  pulse  primary 
transmission  and  a  Q  of  2.4.  The  highest  frequency  at  a  spectral  peak 
and  relative  spectral  power  level  occurs  for  the  lowest  value  of  Q^.  As 
the  value  of  is  increased,  both  the  frequency  at  the  spectral  peak 
and  the  relative  power  levels  decrease.  The  asymmetry  of  the  spectra 
and  the  reduction  in  relative  power  levels  result  from  f?  weighting  in 
the  power  spectral  density  function. 
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The  half-power  bandwidth  of  the  on-axis  secondary 
transmission  is  determined  by  identifying  the  frequencies,  F^  and  F^, 
at  which  the  normalized  power  spectral  density  function,  given  by 
Eq.  (3 . 14) ,  equals  0.5.  Using  Eq.  (A. 18)  in  Eq.  (3.9),  the  con¬ 
dition  which  must  be  satisfied  is  written 


P2xxN(f) 


0.5 


(3.15) 


This  equation  is  transcendental  and  must  be  solved  numerically.  The 
solutions  of  Eq.  (3.15)  are  combined  to  give  a  half-power  bandwidth  of 


F21_F22 I  =  |0.873-2.039|^- 


=  1.166  - 

IT 


(3.16) 


Using  Eq.  (A. 16)  for  n  in  terms  of  primary  bandwidth,  Eq.  (3.15)  becomes 


lF2rF22l  =  °*990|F11-F 


11  12 1 


S  F  -F  1 
1  11  12' 


(3.17) 


This  result  indicates  the  approximate  equivalence  of  the  half-power 
bandwidths  for  the  Gaussian  pulse  primary  transmission  and  the  signal 
produced  on  axis  with  the  broadband  parametric  array  for  plane  wave 
fields.  This  is  the  same  result  which  Muir^  obtained  qualitatively 
for  the  narrowband  parametric  array  for  plane  wave  fields. 

The  quality  factor  of  the  on-axis  secondary  transmission 
may  be  put  in  the  form  of  Eq.  (A. 17)  which,  when  combined  with 
Eqs.  (3.10)  and  (3.17),  gives 
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S  1.2 


(3.18) 


On  axis,  Q2  is  a  constant  which  is  approximately  equal  to  1.2 
for  a  Gaussian  pulse  primary  signal . 

3.  Off-Axis  Radiation 


The  off-axis  pressure  in  the  farfield  of  the  broadband 
parametric  array  for  plane  wave  fields  is  given  in  the  time  domain  as 

before  by  the  convolution  of  the  on-axis  solution  with  the  off-axis 
impulse  response  function.  Similarly  in  the  frequency  domain,  the  off- 
axis  solution  is  given  by  the  product  of  the  on-axis  solution  with  the 
off-axis  transfer  function. 


With  a  Gaussian  pulse  primary  transmission,  the  off-axis 
solution  in  the  time  domain  is  found  by  substituting  Eqs.  (3.4)  and 
(3.6)  into  Eq.  (3.3),  which  gives 


p2(r,0,t) 


r  2 

BP  S 

r  .  ,2  -| 

0  0 

4n2 

4n  (C  “  “  1 

16irp  c  roc. 
o  o  lc 

\  Co> 

n 

r  .  2n 

r  .  .  “ 

x  exp 

-2"  ('  -  r) 

*  K  Ke(f>| 

x  o 

1—  w  — 

(3.19) 


In  the  frequency  domain  the  off-axis  solution  is  found  by  substituting 
Eqs.  (3.2)  and  (3.7)  in  Eq.  (3.1),  which  gives 
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P2(r,0,f) 


8P2S 
o  o 


16irp  c  ra, 
o  o  lc 


-jk2r 

e 


2  1  nr 
“2  nV2 


(3.20) 


Evaluation  of  the  spectral  distortion  .is  a  function  of 
off-axis  angle  involves  the  computation  of  the  normalized  power  spectral 
density  function  and  its  behavior  as  a  function  of  off-axis  angle.  The 
normalized  power  spectral  density  function  is  given  by 


P2xxN<£) 


P2(r,e,f)  p*(r,e,f) 
P2(r,0,f2p)  P*(r,0,f2p) 


(3.21) 


With  the  use  of  Eqs.  (A. 18),  (3.9),  and  (3.20),  the  normalized  power 
spectral  density  function  becomes 


For  a  given  primary  transmission,  the  on-axis  secondary  signal 
has  a  spectral  peak  which  occurs  at  a  secondary  frequency  given  by 
Eq.  (3.10)  and  a  half-power  spectral  width  given  by  Eqs.  (3.16)  or 
(3.17).  As  the  angle  off-axis  increases,  the  spectral  shape  changes 
with  a  downward  shift  in  frequency  at  the  spectral  peak  and  a  reduction 
in  relative  spectral  power  level. 
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V 


The  frequency  at  the  spectral  peak  as  a  function  of  angle  off 
axis  may  be  found  by  maximizing  the  normalized  power  spectral  density 
function  with  respect  to  frequency.  Retaining  the  dependence  on  f2  in 
Eq.  (3.22)  and  differentiating  with  respect  to  frequency  gives 


Setting  Eq.  (3.23)  equal  to  zero  and  solving  for  f2,  the  frequency  at 
the  spectral  peak  is  given  by 


Ve>  ■ 


m 


2 


+ 


(2.668  f 


+  6 


lc 


(3.24) 
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This  equation  is  valid  for  angles  off  axis  only.  On  the  acoustic  axis 
the  frequency  at  the  spectral  peak  is  given  by  Eqs.  (3.9)  and  (3.10). 


B .  Radiation  in  the  Farfield  of  the  Broadband  Parametric  Array  for 
Nonplane  Wave  Fields  with  Low  Amplitude  Primary  Waves 

Some  preliminary  theoretical  considerations  for  the  broadband 
parametric  array  for  nonplane  wave  fields  were  discussed  in  Chapter  II. 
The  analytical  description  for  the  primary  field,  given  by  Eq.  (2.50), 
included  a  narrowband  directivity  function  as  a  simple  multiplier.  This 
analytical  form  only  determines  the  off-axis  amplitude  dependence  of  the 
field  and  does  not  account  for  frequency  dependence  over  the  bandwidth 
of  the  primary  transmission.  If  a  parametric  array  is  formed  with  a 
conventional  linear  source,  then  the  primary  field  will  be  frequency 
dependent  off-axis  over  the  bandwidth  of  the  primary  transmission.  The 
off-axis  properties  of  the  parametric  array  for  nonplane  wave  fields  can 
be  influenced  by  the  frequency  dependence  of  the  primary  field. 

In  this  chapter,  the  broadband  parametric  array  for  nonplane 
wave  fields  will  be  developed  for  a  primary  field  with  frequency  depen¬ 
dent  directivity  functions.  The  limiting  forms  of  this  development  will 
be  considered  and  the  off-axis  spectral  properties  quantified. 

1.  General  Solution 


The  primary  field  associated  with  the  radiation  from  a 
conventional  linear  source  is  discussed  in  Appendix  A  for  the  radiation 
from  a  circular  piston.  This  type  of  field  may  be  expressed  in  the  time 
domain  as 


pi ■  pi(V°'t)  *  hie(t) 


(A.  5) 
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where  p^(ro,0,t)  is  the  on-axis  solution  and  h^(t)  is  the  off-axis 
impulse  response.  An  on-axis  primary  signal  is  assumed  which  is  spheri 
cally  spreading  from  the  origin.  This  signal  is  a  good  representation 
for  the  on-axis  primary  signal  in  the  farfield  of  a  conventional  source 
i.e.,  at  ranges  rQ>Ro,  but  is  a  poor  representation  close  to  the  source 
The  limitations  of  this  signal  representation  in  the  nearfield  is  mini¬ 
mal  in  terms  of  the  analytical  representation  of  the  parametric  array 
for  nonplane  wave  fields  provided  the  dimensions  of  the  conventional 
source  are  not  too  large  in  terms  of  wavelength  at  the  secondary  fre¬ 
quency  and  provided  low  amplitude  primary  levels  are  dominant.  This 
assumption  will  be  discussed  further  elsewhere  in  this  section.  The 
assumed  signal  is  expressed  by 

R  -a1  (r  -R  )  /  r  \ 

pl(ro*0>t)  =  ^  Po  e  °  °  f(C  "  %  jCOS(u,lct"klcro)  *  {3-25) 

where  Pq  is  the  peak  pressure  amplitude  at  Rq.  It  is  also  assumed  that 
the  envelope  function  f(t-rQ/co)  is  slowly  varying  compared  to 
cos(a)^ct-k^crQ)  and  that  the  absorption  coefficient  at  the  primary 
center  frequency  or  carrier  frequency  may  be  used  to  represent  the 
absorption  of  the  bandlimited  signal. 

The  source  strength  density  q,  which  is  a  representation  of 
the  strength  of  the  virtual  sources  distributed  in  the  interaction 
volume,  is  found  from  Eq.  (2.22)  with  Eqs.  (3.25)  and  (A. 5),  and  is 
expressed  by 


‘*(ro*0»to) 


2  2 
6R  P 
o  o 

2  4  2 
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o  o  o 


-2al(ro-Ro) 
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)  cos(“lcto_klCro)  *  h10(to) 


(3.26) 
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In  the  frequency  domain,  the  source  strength  density  becomes 


^  | ^ <ro  * 6  *  t0) 


jw-eR2P2  -2a. (r  -R  ) 
2  o  o  loo 

2  4  2  6 
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o  o  o 


x  fF  j  f  t - ]  cos  (to.  t  -k.  r  )  *  h1Q(t  ) 

\  o  cq J  lc  o  lc  o  10  o 


*  (Fjflt - I cos(to.  t  -k..  r  )  *  h  (t  ) 

\  o  cq  j  lc  o  lc  o  10  o 


(3.27) 


The  frequency  domain  representation  of  the  source  strength  density 
function  involves  the  self -convolution  of  the  primary  pressure  in  the 
frequency  domain.  Using  the  shift,  convolution,  and  modulation 
theorems,  Eq.  (3.27)  may  be  written 


^h(r0*9>to)  = 


jui  8R2P2  -2a.  (r  -R  ) 

2  o  o  loo 

2  4  2  e 
per 
o  o  o 


*  Te  (F<fl-£lc>  +  F<£l+flc)1  »10<fl> 


*  2  e  ^F££i-fic^  +  F(fi+,1C)I  Hie(fiJ 


(3.28) 


In  this  equation  F(f)  is  the  Fourier  transform  of  the  primary  envelope 

function  given  by  £F{f(t)}  and  H10(f)  is  the  off-axis  transfer  function 

for  the  primary  transmission  given  by  «F{h--(t)}. 

W  10 


When  written  as  a  convolution  integral,  the  source  strength 
density  in  Eq.  (3.28)  may  be  expressed  as  the  sum  of  four  terms  which 
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are  the  modulation  products  resulting  from  the  nonlinear  process.  Each 
of  the  four  terms  is  a  convolution  integral  which  may  be  expressed  as  a 
correlation  integral  utilizing  the  properties  of  the  Fourier  transform 
of  real  signals.  Equation  (3.28)  then  becomes 
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o  o  o 
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■/—  00 

+  f  f  (f1+flc)  f(f2+fi+fic)  Hie(fi)  Hie(fi+f2)  dfi 
•'—00 

(3.29) 

Each  of  these  correlation  integrals  is  a  function  of  frequency  shift  or 
frequency  difference  which  is  a  secondary  frequency  component.  The  mod¬ 
ulation  products  which  each  of  the  integrals  represent  may  be  determined 
by  evaluation  of  the  integrand  and  recalling  that  F(f)  is  symmetric  and 
has  a  peak  at  f=0  and  that  H^0(f)=l  on-axis.  The  first  integral  yields 
a  baseband  secondary  spectrum  which  results  from  the  interaction  of 
positive  primary  frequency  components.  The  second  and  third  integrals 
yield  the  negative  and  positive  secondary  spectra,  respectively,  cen¬ 
tered  at  the  second  harmonic  of  the  carrier  frequency.  The  last 
integral  yields  another  baseband  secondary  spectrum,  which  results  from 


54 


interaction  of  negative  primary  frequency  components.  Retaining  the 
lowest  frequency  components  of  the  source  strength  density  and 
restricting  the  analysis  to  positive  frequencies,  the  source  strength 
density  in  the  frequency  domain  reduces  to  a  function  of  the  first  cor¬ 
relation  integral  only  and  is  given  by 


q(rQ,e,to) 


J“268oPo 

,242 
2p  c  r 
o  o  o 


-2a, (r  -R  ) 
loo 
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*  f  F*(fr£lc)  F<fl+f2-flc>  “l6(fl+f2>  d£l  • 

—  CO 

(3.30) 

The  general  solution  for  the  secondary  pressure  field  produced 
by  the  interaction  of  primary  field  components  is  given  in  the  frequency 
domain  by  Eq.  (2.27).  Using  the  representation  of  the  source  strength 
density  in  the  frequency  domain  given  by  Eq.  (3.30),  the  general  solu¬ 
tion  is  written 


P2(r,0,f)  = 
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o 

This  four-dimensional  solution  is  a  superposition  of  first  order 
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interaction  components  determined  by  the  primary  frequency  spectra.  The 

amplitude  dependence  of  each  interaction  component  results  in  part  from 

the  volumetric  integration  over  the  virtual  sources.  The  representation 

of  the  secondary  field  given  by  Eq.  (3.31)  is  general  for  nonplane  wave 

fields  and  similar  in  form  to  Eq.  (2.46),  which  is  a  representation  of 

7  48 

the  broadband  parametric  array  for  plane  wave  fields  used  by  Muir  * 

49 

and  Muir  and  Blue. 

Equation  (3.31)  was  derived  assuming  primary  radiation  which 
is  spherically  spreading  from  the  origin.  This  assumption  is  valid  if 
the  effects  of  the  aperture  factor  can  be  ignored  in  the  nearfield  of 
the  primary  radiation.  The  aperture  factor  which  was  discussed  in 
Chapter  II  is  the  directivity  function  at  the  secondary  frequency 
associated  with  the  effective  cross  section  of  the  interaction  volume. 
The  aperture  factor  which  results  in  the  farfield  of  the  primary  radia¬ 
tion  is  a  diffraction  effect  for  diverging  waves  and  is  included  in  the 
formalism  of  Eq.  (3.31).  An  implicit  assumption  which  is  included  in 
the  derivation  is  placed  on  the  off-axis  transfer  function.  This  func¬ 
tion  must  be  valid  throughout  the  interaction  volume.  Since  the  aper¬ 
ture  factor  is  ignored  in  the  nearfield  of  the  primary  radiation,  the 
specification  of  an  off-axis  transfer  function  which  Is  valid  only  in 
the  farfield  of  the  primary  radiation  is  required. 

2.  Farfield  Solution 

The  volumetric  integration  in  Eq.  (3.31)  is  fundamental  to 
solutions  for  the  parametric  array  for  nonplane  wave  fields.  It  has 
been  evaluated  for  the  narrowband  parametric  array  by  numerical  integra¬ 
tion  and  by  making  various  approximations  and  simplifying  assumptions 
in  the  integrand  to  fit  specific  physical  situations.  For  the  current 
analysis  a  solution  to  Eq.  (3.31)  is  sought  in  the  farfield  of  the  para¬ 
metric  array.  A  farfield  solution  is  indicative  of  many  trends  which 
occur  in  the  nearfield  of  the  parametric  array.  These  trends  remain 
after  the  energy  associated  with  the  primary  frequencies  has  been 
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removed  by  absorption  in  the  medium  and  the  farfield  of  the  parametric 
array  is  attained. 

Conventional  spherical  geometry  is  assumed  and  indicated  in 
Fig.  3.5  where  the  acoustic  axis  is  along  the  z-axis.  The  subscripted 
variables  indicate  position  of  virtual  volumetric  sources  in  the  pri¬ 
mary  field  and  unprimed  variables  indicate  position  relative  to  the 
field  points.  It  is  assumed  that  the  field  point  is  far  from  the  inter¬ 
action  volume  or  r>>r  .  Under  this  assumption  the  Fraunhofer  or  far- 

o 

field  approximation  is  justified  where  R  is  given  by  two  terms  of  the 
binomial  expansion  of  the  law  of  cosines  approximate  representation,  or 

R  =  r  -  r  cosy  +  ...  .  (3.32) 

o 

In  this  relation  y  is  the  angle  between  the  vector  to  the  source  point 
and  the  vector  to  the  field  point  and  is  defined  by 

cosy  =  cos9  cos9  +  sin9  sin0  c.os(4>->!’  )  .  (3.33) 

o  o  o 

In  the  numerator  of  the  Green's  function,  R  is  replaced  by  the  first  two 
terms  in  Eq.  (3.33)  and  in  the  denominator  by  only  the  first  term. 

After  making  the  farfield  approximations,  including  the  differential 
volume  in  spherical  coordinates,  and  collecting  terms,  Eq.  (3.31) 
becomes 
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FIGURE  3.5 

SPHERICAL  GEOMETRY  FOR  ANALYSIS  OF  A  PARAMETRIC  ACOUSTIC  ARRAY 
FOR  NON  PLANE  WAVE  FIELDS  WITH  THE  ACOUSTIC  AXIS  ALONG  THE  z-AXIS 
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where  ct^c^-c^.  As  a  result  of  the  use  of  the  farfield  approximations, 
the  range  integral  is  trivial  and  may  be  expressed  in  closed  form.  The 
secondary  pressure  in  the  farfield  is  written 
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where  the  dependence  of  H^0(f)  on  angles  has  been  included  explicitly. 

The  surface  integral  in  Eq.  (3.35)  is  similar  to  a  form  derived  by 
25 

Berktay  and  Leahy  for  the  narrowband  parametric  array.  It  has  been 

referred  to  as  the  convolution  integral  for  the  parametric  array  after 
34 

Blue  assumed  different  geometry,  made  small  angle  approximations,  and 
derived  the  result  given  by  Eq.  (2.38),  which  involves  the  double  con¬ 
volution  of  the  weighted  product  of  the  primary  beam  functions  with  the 
Rutherford  directivity. 

The  solution  for  the  secondary  pressure  on  the  acoustic  axis 
is  given  by  P2(r,0,f)  and  obtained  by  setting  0=0  in  Eq.  (3.35).  In 
this  case,  cosy  reduces  to  cosQ  .  The  off-axis  transfer  function  H  (f) 
may  be  computed  from  Eq.  'v3.1)  and  the  normalized  power  spectral  density 
function  Pox  jj(f)  is  given  by  Eq.  (3.21). 

A  simplification  of  the  general  farfield  solution  for  the 
broadband  parametric  array,  Eq.  (3.35),  is  obtained  when  the  directivity 
or  off-axis  transfer  functions  have  circular  symmetry  about  the  acoustic 
axis  or  an  independence  with  respect  to  <{i.  Referring  to  Fig.  3.5  for 
this  case,  y  is  given  by 

cosy  =  cos©  cosQ  +  sin0  sin0  cosiF  ,  (3.36) 

o  o  o 

where  an  arbitrary  rotation  of  coordinates  about  the  acoustic  axis 
removes  the  ip  dependence.  Under  these  conditions  the  integral  over  tp 
may  be  performed  and  Eq.  (3.35)  becomes 
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(3.37) 


The  space  integral  in  this  equation  is  similar  to  a  result  which  Berktay 
25 

and  Leahy  derived  for  the  narrowband  parametric  array.  Equation  (3.37) 
is  the  general  solution  in  the  farfield  of  the  broadband  parametric 
array  for  nonplane  wave  fields  with  directivity  or  off-axis  transfer 
functions  which  are  symmetric  about  the  acoustic  axis.  With  the  far¬ 
field  assumption  and  the  assumption  of  circular  symmetry,  the  four¬ 
dimensional  integral  in  the  general  solution  given  by  Eq.  (3.31)  is 
reduced  to  a  two-dimensional  integral  in  Eq.  (3.37).  Further  analytical 
reduction  for  cases  of  practical  interest  is  difficult  but  may  be 
attempted  by  appropriate  selections  of  F(f)  and  H  (f  )  .  A  more 
straightforward  approach  is  evaluation  by  double  numerical  integration. 
This  is  accomplished  by  separating  Eq.  (3.37)  into  real  and  imaginary 
parts,  performing  separate  double  integrations  and  computing  the  magni¬ 
tude  of  the  result. 


Physical  insight  into  the  characteristics  of  the  broadband 

parametric  array  for  nonplane  wave  fields  may  be  gained  by  examining 

several  limiting  forms  of  Eq.  (3.35)  which  occur  as  a  result  of  dominant 

directional  characteristics.  The  characteristics  of  a  parametric  source 

which  is  formed  by  radiation  from  a  conventional  linear  source  are 

always  a  combination  of  effects  due  to  the  primary  directivity  or  off- 

axis  transfer  functions  and  the  Rutherford  directivity.  Limiting  forms 
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of  the  convolution  integral  have  been  evaluated  by  Berktay  and  Leahy. 
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There  are  two  limiting  conditions  associated  with  Eq.  (3.35)  which  are 
important.  The  first  condition  occurs  when  the  parametric  array  is 
diffraction  limited  or  strongly  dominated  by  the  product  of  the  direc¬ 
tivity  or  off-axis  transfer  functions  at  the  primary  frequencies.  The 
second  condition  occurs  when  the  parametric  array  is  absorption  limited 
and  the  Rutherford  directivity  function  is  strongly  dominant. 

If  the  primary  directivity  functions  are  extremely  narrow, 

then  0  is  small  and 
o 


cosy  -  cos0 


(3.38) 


Under  this  condition  Eq.  (3.35)  becomes 


i4eR2P2  2a.  R  -(ou+jk  )r  * 

p„(r,e,f)  -  .  1  °  e  2  2  J_j  (W  F(f1+f2-fic> 


8^P  CnraT 
o  o  i 


1  lc' 


2tt  v/2 


1  +  j  —  (1-COS0)  o  o 


/  f  Hie(fi’9o,'l’o)  Hie(fi+f2,0o’'i’o)  sin0odWfr 


(3.39) 


The  term  which  emerges  from  the  convolution  integral  is  simply  the 
Rutherford  directivity  function  which  was  identified  in  Eq.  (2.34).  The 
primary  beams  act  as  spatial  impulse  functions  with  strength  given  by 
the  remaining  integrals  over  the  interaction  volume  in  Eq.  (3.39).  The 
strength  of  the  impulse  function  vanishes  as  the  primary  beamwidths  go 
to  zero.  This  trend  leads  to  a  physically  unrealizable  situation  which 
is  supported  by  operational  considerations  of  the  parametric  array.  If 
a  parametric  array  is  operated  with  constant  radiated  acoustic  power, 
then  the  on-axis  source  level  at  the  primary  frequencies  becomes 
infinite  as  the  primary  beamwidths  go  to  zero.  Conversely,  if  a 
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parametric  array  is  operated  with  constant  on-axis  source  levels  at  the 
primary  frequencies,  then  the  radiated  acoustic  power  vanishes  as  the 
primary  beamwidths  go  to  zero.  This  limiting  form  is  simply  a  spherical 
wave  representation  of  the  broadband  parametric  array  for  plane  wave 
fields  which  was  discussed  in  detail  in  the  previous  section.  In  this 
limit  of  Eq.  (3.35),  the  assumption  of  primary  waves  which  are  spheri¬ 
cally  spreading  from  the  origin  results  in  physical  inaccuracies.  If 
radiation  is  considered  from  a  conventional  source  and  the  primary  beams 
are  extremely  narrow,  then  the  primary  source  has  a  long  nearfield 
characterized  by  plane  wave  behavior  which  is  generally  confined  to  the 
volume  in  front  of  and  determined  by  the  area  of  the  source.  In  this 
case,  most  of  the  nonlinear  interaction  takes  place  in  the  Fresnel 
diffraction  region  of  the  source,  the  parametric  array  is  absorption 
limited,  and  the  theoretical  development  for  the  broadband  parametric 
array  for  plane  wave  fields  is  applicable. 

The  other  limiting  case  of  interest  is  identified  when  the 
half-power  beamwidth  associated  with  the  Rutherford  directivity  is  small 
compared  to  the  half -power  beamwidths  of  the  primary  radiations.  The 
evaluation  of  this  special  case  is  aided  by  the  adoption  of  the  geometry 
shown  in  Fig.  3.6.  This  geometry  differs  from  that  of  Fig.  3.5  in  that 
the  field  point  rather  than  the  acoustic  axis  lies  on  the  z-axis.  If 
the  half-power  beamwidth  of  the  Rutherford  directivity  is  small,  the 
convolution  integral  in  Eq.  (3.35)  may  be  written 


I  = 


2it  it/ 2 

If 


H10 (f l’°)  H19(fl+f2’a)  Sin0o 


1  +  j 


-5-  (1-COS0  ) 

32  0 

HP 


d9  dt|i 
o  o 


(3.40) 


With  the  spherical  geometry  in  Fig.  3.5,  the  primary  directivity 
functions  or  off-axis  transfer  functions  are  simple  functions  of  angle 
while  the  Rutherford  directivity  is  a  complicated  function  of  angle. 
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FIGURE  3.6 

SPHERICAL  GEOMETRY  FOR  ANALYSIS  OF  A  PARAMETRIC  ACOUSTIC  ARRAY 
FOR  NON  PLANE  WAVE  FIELDS  WITH  THE  FIELD  POINT  ON  THE  z  AXIS 
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The  new  geometry  in  Fig.  3.6  reverses  this  functional  dependence.  The 
directivity  for  the  parametric  array  is  now  a  function  of  £  and  f. ,  which 
are  angles  that  determine  the  position  of  the  acoustic  axis.  The  direc¬ 
tivity  function  at  the  primary  frequencies  is  a  function  of  o,  given  by 


coso  =  cos©  cos£  +  sin©  sin£  cos(ii>  -e)  .  (3.41) 

o  o  o 


As  ©  becomes  small,  the  convolution  integral  of  Eq.  (3.40) 
HP 

has  significant  value  only  when  0q  is  small  also,  when  the  elemental 
volume  is  near  the  vector  to  the  field  point.  Under  these  conditions  a 
may  be  approximated  by 


coso  -  cos£ 


(3.42) 


and  the  value  of  H10(f)  which  contributes  to  the  integral  is  just  the 
value  near  the  fixed  angle  £.  The  convolution  integral  becomes 


1  =  Hie^fi’^  Hie^fi+f2’^ 


2tt  tt/2 

//; 


J  J  k2 

o  o  1  +  j  —  (1-cosS  ) 

J  aT  o' 
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(3.43) 


After  performing  the  integrals  in  Eq.  (3.43),  the  general  result  for 
this  limiting  case  is  obtained  from  Eq.  (3.35)  and  written 


2a.  R  -(a,+jk,)r 
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df  •  (3-44) 


When  the  half-power  beamwidth  corresponding  with  the  Rutherford 
directivity  is  small  in  general  and  also  small  compared  with  the  half¬ 
power  beamwidths  at  the  primary  frequencies,  the  Rutherford  directivity 
acts  as  a  spatial  impulse  function  and  the  parametric  array  is  diffrac¬ 
tion  limited.  Most  of  the  nonlinear  interaction  takes  place  in  the 
Fraunhofer  region  of  the  primary  source  and  the  directivity  of  the  para¬ 
metric  array  is  determined  by  the  product  of  the  directivity  or  off- 
axis  transfer  functions  at  the  primary  frequencies.  For  this  limiting 
form,  the  assumption  of  spherically  spreading  primary  waves  from  the 
origin  does  not  compromise  validity  of  the  result  given  by  Eq.  (3.44). 

In  deriving  Eq.  (3.44)  it  was  assumed  that  the  half -power 
beamwidth  was  small  in  general.  According  to  Eq.  (2.35),  this  assump¬ 
tion  is  valid  as  the  absorption  at  the  primary  center  frequency  goes 
to  zero  and  as  the  secondary  frequency  goes  to  infinity.  Both  of  these 
conditions  are  physically  unreasonable.  The  other  assumption,  which 
requires  that  the  half-power  beamwidth  associated  with  the  Rutherford 
directivity  be  small  compared  to  the  half-power  beamwidths  associated 
with  the  primary  frequencies,  is  physically  reasonable.  It  is  this  con¬ 
dition  for  which  the  limiting  form  given  by  Eq.  (3.44)  should  be 
considered . 

It  should  be  noted  that  the  limiting  form  for  broadband 
radiation  at  the  primary  frequencies  is  linearly  dependent  on  the  secon¬ 
dary  frequency  in  contrast  to  the  quadratic  dependence  on  secondary  fre¬ 
quency  for  the  broadband  parametric  array  for  plane  wave  fields.  This 
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trend  was  noted  by  Berktay  and  Leahy  and  by  Blue,  and  has  strong 
influence  on  the  broadband  parametric  array.  In  going  from  an  absorp¬ 
tion  limited  parametric  array  to  a  diffraction  limited  parametric  array, 
a  shift  in  the  power  spectra  toward  lower  frequencies  will  occur  for  a 
given  primary  transmission. 
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3. 


Farfield  Solution  with  a  Gaussian  Pulse  Primary  Transmission 
from  a  Circular  Piston 


The  general  solutions  for  the  secondary  pressure  fields  which 
result  with  the  broadband  parametric  array  for  nonplane  wave  fields  are 
difficult  to  evaluate.  Analyses  of  the  limiting  forms  of  these  solu¬ 
tions  are  helpful  especially  when  the  primary  directivity  acts  as  a 
spatial  impulse  function  and  the  broadband  parametric  array  for  plane 
wave  fields  is  recovered.  The  results  for  this  case  are  described  in 
Section  III. A.  When  the  Rutherford  directivity  acts  as  a  spatial 
impulse  function,  another  limiting  form  is  recovered  where  the  product 
of  the  off-axis  transfer  functions  is  dominant  in  determining  the  secon¬ 
dary  field.  The  detailed  characteristics  of  the  secondary  sound  field 
with  the  parametric  array  for  nonplane  wave  fields  may  be  described  more 
concisely  by  considering  a  particular  primary  pulse  transmission  from  a 
given  conventional  linear  source. 

The  characteristics  of  a  Gaussian  pulse  primary  transmission 
from  a  circular  piston  is  considered  in  detail  in  Appendix  A.  This 
transmission  is  mathematically  tractable  and  is  characterized  by  a  spec¬ 
trum  with  a  single  frequency  peak  on-axis.  The  transmission  from  a  cir¬ 
cular  piston  source  is  symmetric  about  the  acoustic  axis;  therefore  the 
solution  for  secondary  pressure  in  the  farfield  of  the  broadband  para¬ 
metric  array  is  obtained  from  Eq.  (3.37).  Utilizing  the  properties  of 
the  Gaussian  pulse  transmission  given  by  Eqs.  (A. 13)  and  (A. 18)  and  the 
circular  piston  given  by  Eq.  (A. 4),  the  solution  is  written 
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(3.45) 


In  this  form  the  secondary  pressure  depends  on  the  quality  factor  of 
the  primary  transmission  and  the  primary  center  frequency  f  .  For  a 
given  f^c,  the  secondary  pressure  is  directly  proportional  to  Q  .  This 
behavior  is  a  fundamental  result  of  signal  analyses  which  indicates 
that,  if  the  total  energy  in  the  primary  transmission  is  constant,  then 
energy  will  be  concentrated  for  a  high  transmission  over  a  narrow 
frequency  band  for  both  the  primary  and  secondary  transmissions.  As 
is  decreased,  the  energy  is  spread  over  a  broader  band  and  the  level  at 
a  given  frequency  component  is  reduced.  For  a  constant  transmission, 
the  secondary  pressure  is  inversely  proportional  to  the  step-down  ratio 
(ratio  of  primary  center  frequency  to  secondary  frequency). 


Equation  (3.45)  is  well  suited  for  numerical  analysis  and  was 
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evaluated  using  a  double  numerical  integration  routine  which  employs 
an  adaptive  variable  grid  iterative  technique.  The  results  of  the  eval¬ 
uation  are  given  in  Chapter  4  where  computations  with  Eq.  (3.45)  are 
compared  with  experimental  results. 
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Examination  of  the  limiting  forms  of  Eq .  (3.45)  offer  physical 
insight  into  the  characteristics  of  the  secondary  field.  When  the 
Rutherford  directivity  is  dominant,  the  absorption  limited  parametric 
array  for  plane  wave  fields  is  recovered.  When  the  directivity  or  off- 
axis  transfer  functions  at  the  primary  frequencies  are  dominant,  the 
solution  in  this  limit  is  given  by  Eq.  (3.44)  with  appropriate  substitu¬ 
tions,  and  written 
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(3.46) 


On  the  acoustic  axis  where  9=0,  Eq.  (3.46)  reduces  to 
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This  on-axis  solution  for  the  extreme  case  of  the  diffraction  limited 
broadband  parametric  array  with  a  Gaussian  pulse  primary  transmission 
may  be  compared  with  Eq.  (3.7)  which  is  the  on-axis  solution  for  the 
broadband  parametric  array  for  plane  wave  fields  with  a  Gaussian  pulse 
transmission.  The  terms  in  Eq.  (3.47)  involving  are  slowly 

varying  with  respect  to  secondary  frequency  for  values  of  ^/“.p^l, 
which  are  of  practical  importance.  Comparison  of  Eqs.  (3.47)  and  (3.7) 
indicates  that  the  spectral  shapes  for  the  two  limiting  cases  are  simi¬ 
lar  except  for  the  shift  from  quadratic  to  linear  dependence  on  secon¬ 
dary  frequency  in  going  from  absorption  to  diffraction  limited 
parametric  arrays.  Equation  (3.47)  yields  a  spectral  peak  at  a  fre¬ 
quency  given  by 


2.668  flc 

TT  Qj 


(3.48) 


This  result  is  the  peak  frequency  for  the  on-axis  secondary  transmission 
produced  in  the  extreme  case  of  a  diffraction  limited  broadband  para¬ 
metric  array.  The  peak  frequency  given  by  Eq.  (3.48)  may  be  compared 
with  the  result  given  by  Eq.  (3.9)  for  the  broadband  parametric  array 
for  plane  wave  fields.  This  comparison  indicates  that  the  peak 
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secondary  frequency  on-axis  is  higher  for  the  absorption  limited 
braodband  parametric  array  by  a  factor  of  /2. 

The  off-axis  transfer  function  H„.(f)  for  the  extreme 
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diffraction  limited  parametric  array  is  obtained  by  using  Eqs.  (3.46) 
and  (3.47)  in  Eq.  (3.1).  The  off-axis  transfer  function  is  given  by 


J^[k^a  sin0]  [  (k^+k^a  sin0] 

[k^a’sin0]  [(k^+k2)a  sin0]  dfl  '  (3.49) 

The  off-axis  transfer  function  describes  the  spectral  distortion  which 
occurs  at  angles  with  respect  to  the  acoustic  axis.  The  square  of  the 
magnitude  of  the  off-axis  transfer  function  is  a  representation  of  the 
secondary  power  spectrum  if  on-axis  this  spectrum  is  uniform  over  all 
secondary  frequencies. 

T  e  off-axis  transfer  function  was  evaluated  for  a  fixed  value 

of  and  f^.  In  Fig.  3.7,  the  frequency  dependence  is  examined  as  a 

function  of  off-axis  angle.  The  off-axis  angle  is  included  in  the 

dimensionless  parameter  (a  sin0  )/c  .  On  the  acoustic  axis,  H0.(f)  is 

oo  Z0 

uniform  as  a  function  of  frequency.  Off-axis  1^(0  decreases  at>  all 
frequencies  with  only  a  slightly  larger  decrease  at  the  higher  fre¬ 
quencies.  This  trend  persists  until  a  value  of  off-axis  parameter  of 
1.05,  where  H2g(f)  shows  a  decrease  and  then  an  increase  in  level  as  a 
function  of  frequency.  At  the  next  value  of  off-axis  parameter  1.31, 
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FIGURE  3.7 

OFF-AXIS  POWER  TRANSFER  FUNCTION  venui  FREQUENCY 
FOR  THE  EXTREME  DIFFRACTION  LIMITED  BROADBAND 
PARAMETRIC  ARRAY  FOR  NON  PLANE  WAVE  FIELDS 

f (t)  -  GAUSSIAN,  Ht0«)  -  CIRCULAR  PISTON 
Q.  =  10.91  f.  =600  kHz 
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H 2 0 ( ^ )  a8ain  exhibits  the  original  pattern.  This  behavior  of  the 
off-axis  transfer  function  is  better  understood  in  Fig.  3.8,  where 
H^^Cf)  is  plotted  as  a  function  of  the  off-axis  parameter.  There  are 
two  curves  on  the  plot.  For  zero  secondary  frequency  H20(f)  has 
characteristics  of  a  conventional  off-axis  transfer  function  exhibited 
by  primary  and  secondary  lobes.  For  a  secondary  frequency  of  100  kHz, 
the  curves  are  almost  the  same  with  a  slight  divergence  around  the 
amplitude  null  and  a  difference  of  about  3  dB  at  the  common  secondary 
peak.  This  off-axis  transfer  function  for  the  broadband  parametric 
array  in  the  extreme  diffraction  limit  is  very  unconventional.  It  is 
acute  in  angular  response  but  is  very  insensitive  in  frequency  response. 
In  this  limit  the  broadband  parametric  array  exhibits  very  stable 
directional  characteristics  which  are  almost  frequency  independent. 

One  feature  which  may  be  noted  in  the  100  kHz  curve  is  the  small  peak 
near  a  value  of  1.31  for  the  off-axis  parameter.  This  peak  occurs  in 
integrating  across  the  nulls  of  two  off-axis  transfer  functions  in  the 
primary  band,  which  differ  in  frequency  by  100  kHz. 

A  representation  of  Eq.  (3.46)  which  is  applicable  over  the 
main  beam  of  the  off-axis  transfer  function  may  be  obtained  by  approxi¬ 
mating  the  principal  maximum  of  H  (f)  for  the  circular  piston  with  a 
Gaussian  off-axis  transfer  function.  The  approximate  H1Q(f)  is  given 
by 


H18 l)  =  eXp  “fklaK  sin0ol2  > 

where  K  is  a  constant  which  is  chosen  to  fit  the  principal  maximum  at  a 
given  amplitude. 

If  it  is  required  that  the  Gaussian  approximation  fit  the 
principal  maximum  at  the  one-half  amplitude  points,  then  A  must  satisfy 
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Equation  (3.51)  is  satisfied  when  K=0.33.  Substituting  this  function  in 
Eq.  (3.44),  the  secondary  field  for  the  extreme  diffraction  limited 
parametric  array  with  a  Gaussian  primary  pulse  transmission  and  a 
Gaussian  off-axis  transfer  function  is  expressed  in  closed  form  by 


P1(r,9,f) 
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where 

c1  =  ~  2TraK  sin9  >  and 
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C2  2.668  f,  ' 

lc 

The  on-axis  solution  P2(r,0,f)  with  the  Gaussian  approximation  is  equiv- 
lent  to  Eq.  (3.47)  which  is  the  on-axis  solution  previously  derived  for 
this  limit  assuming  the  off-axis  transfer  function  for  the  circular 
piston. 


The  off-axis  transfer  function  H„  (f)  with  the  Gaussian 

4o 

approximation  is  computed  from  Eq.  (3.1)  using  Eqs.  (3.52)  and  (3.47), 
and  is  given  by 
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H20(f)  =  C2(C1+C2) 

2C22f2lc  .  (3.53) 

Equation  (3.53)  was  evaluated  for  the  same  parameters  which 
were  used  to  evaluate  Eq.  (3.49).  The  frequency  dependence  of  H-.(f)  in 
the  Gaussian  approximation  is  examined  in  Fig.  3.9  as  a  function  of  off- 
axis  angle.  A  comparison  of  Figs.  3.7  and  3.9  indicates  similar  trends 
for  both  transfer  functions  with  the  exception  that  there  is  no  case 
shown  in  Fig.  3.9,  where  the  transfer  function  deviates  from  smooth 
regular  behavior  as  it  did  in  Fig.  3.7.  The  performance  of  H20(f)  in 
the  Gaussian  approximation  is  examined  as  a  function  angle  in  Fig.  3.10 
for  secondary  frequencies  of  0  and  100  kHz.  A  comparison  of  this 
figure  with  Fig.  3.8  again  indicates  similar  trends  to  a  value  of  off- 
axis  parameters  of  about  0.79.  At  this  value  the  trends  diverge,  with 
H  (f)  for  the  Gaussian  approximation  continuing  to  decrease  monotoni- 
cally  while  H20(f)  for  the  circular  piston  decreases  rapidly  and  then 
increases  in  response  to  the  sidelobes  in  H^g(f).  The  agreement  between 
the  two  functions  for  H  (f)  is  reasonably  good  to  the  -30  dB  power 
level  over  the  principal  maximum. 

The  analyses  which  have  been  presented  have  addressed  the 
broadband  parametric  array  for  nonplane  wave  fields  in  the  extreme 
diffraction  limit.  As  the  characteristics  of  the  parametric  transmis¬ 
sion  shift  from  the  absorption  limit  to  the  diffraction  limit,  the 
dependence  on  secondary  frequency  shifts  from  quadratic  to  linear.  With 
this  shift  in  frequency  dependence  an  accompanying  shift  toward  lower 
frequencies  in  the  power  spectra  will  occur.  This  shift  is  character¬ 
ized  by  a  reduction  in  the  frequency  at  the  spectral  peak  by  a  factor  of 
^2  on-axis  for  a  Gaussian  transmission.  In  the  diffraction  limit,  the 
specific  primary  off-axis  transfer  function  H^(f)  is  dominant  in 
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FIGURE  3.9 

OFF-AXIS  POWER  TRANSFER  FUNCTION  versus  FREQUENCY 
FOR  THE  EXTREME  DIFFRACTION  LIMITED  BROADBAND 
PARAMETRIC  ARRAY  FOR  NON  PLANE  WAVE  FIELDS 

f  (t)  -  GAUSSIAN,  H,0(f)  -  GAUSSIAN  APPROXIMATION 
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FIGURE  3.10 

OFF-AXIS  POWER  TRANSFER  FUNCTION  versus  OFF-AXIS  ANGLE 
FOR  THE  EXTREME  DIFFRACTION  LIMITED  BROADBAND 
PARAMETRIC  ARRAY  FOR  NON  PLANE  WAVE  FIELDS 

f(t)  -  GAUSSIAN,  H10(f)  -  GAUSSIAN  APPROXIMATION 
Q1  =  10.91  f1c  =  600  kHz 


ARL:UT 
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CHAPTER  IV 


EXPERIMENTAL  STUDY  OF  THE  BROADBAND 
PARAMETRIC  ARRAY 

The  objective  of  the  experimental  study  is  to  measure  the 
characteristics  of  the  broadband  parametric  array  as  a  function  of 
various  parameters.  Particular  emphasis  is  directed  toward  the 
behavior  of  the  broadband  parametric  array  as  a  function  of  pri¬ 
mary  directivity  or  off-axis  transfer  functions. 

The  experimental  approach  includes  the  use  of  a  variable 
aperture  circular  piston  source.  This  source  provides  variable  primary 
directivity  and  allows  parameter  values  such  as  primary  center  frequency 
and  transducer  transfer  function  to  be  held  constant.  Highly  directive 
primary  radiation  with  most  of  the  nonlinear  interaction  occurring  in 
the  Fresnel  diffraction  region  of  the  primary  field  represents  a  physi¬ 
cal  approximation  of  the  absorption  limited  broadband  parametric  array 
for  plane  wave  fields.  As  the  primary  directivity  decreases,  the  tran¬ 
sition  to  the  diffraction  limited  broadband  parametric  array  for  non¬ 
plane  wave  fields  occurs  and  the  nonlinear  interaction,  which  takes 
place  in  the  Fraunhofer  diffraction  region,  increases.  The  transition 
from  plane  wave  to  nonplane  wave  fields  occurs  most  rapidly  in  going 
from  values  of  0.5  to  1.5  in  the  directivity  ratio  (ratio  of  primary 
half-power  beamwidth  to  Rutherford  bearawidth) . 


Experimental  parameter  values  were  selected  that  allow  the 
formation  of  broadband  parametric  arrays  which  transition  from  a  domi¬ 
nance  in  the  absorption  limit  to  a  dominance  in  the  diffraction  limit. 
The  aperture  factor,  which  was  briefly  discussed  in  Chapters  II  and  III, 
has  greater  influence  on  the  off-axis  properties  of  a  broadband  para¬ 
metric  transmission  as  the  directivity  of  the  primary  source  increases. 
Experimental  parameter  values  were  selected  to  minimize  the  effects  of 
the  aperture  factor  and  restrict  the  effects  to  the  -15  dB  level  or  less 
in  the  secondary  beam  pattern.  Other  constraints  that  affect  the 
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selection  of  experimental  parameters  include  the  restriction  of  step- 
down  ratio  to  a  value  of  10  or  less  and  the  determination  of  primary 
center  frequency  so  that  the  nearfield  extent  of  the  parametric  array 
is  mainly  confined  to  measurement  ranges  that  are  accessible  at  the 
fresh  water  test  site. 

The  experimental  parameter  values  that  were  selected  are  given 
in  Table  IV. 1.  The  center  frequency  for  the  primary  transmission  was 
chosen  to  be  600  kHz.  With  this  frequency,  the  parametric  array  length 
is  from  35.7-47.6  m  based  on  absorption  in  fresh  wa^er  from  Schulkin  and 
Marsh^  and  the  change  in  water  temperature  during  the  test  period.  A 
circular  piston  with  a  selectable  radius  is  specified  for  the  source 
transducer.  Four  source  configurations  and  directivities  are.  indicated 
in  Table  IV. 1.  The  values  of  primary  pulse  length,  primary  pulse  qual¬ 
ity  factor,  primary  center  frequency,  and  step-down  ratio  are  inter¬ 
related  and  were  selected  after  evaluation  of  the  transmit  performance 
of  the  selectable  radius  circular  piston  projector.  The  Rutherford 
beamwidth  at  66  kHz  and  the  primary  product  beamwidths  are  given  for 
comparison. 

A.  Equipment  and  Procedures 

The  experimental  study  was  conducted  in  fresh  water  at 
Lake  Travis  Test  Station  (LTTS) ,  an  instrumented,  multiple  barge  facil¬ 
ity  moored  in  about  30  m  of  water.  The  physical  configuration  for  the 
experiments  consisted  of  a  mechanically  instrumented  column  on  which  the 
selectable  radius  circular  piston  projector  was  mounted  at  a  depth  of 
4.48  m,  and  two  portable  columns  on  which  the  hydrophones  were  mounted 
at  a  depth  of  4.48  m  also.  The  acoustic  paths  were  primarily  underneath 
the  floats  on  which  the  barges  rest. 

A  block  diagram  of  the  experiments  is  shown  in  Fig.  4.1.  The 
transmit  section  consists  of  an  oscillator  whose  output  is  multiplied  by 
a  Gaussian  envelope  in  the  shaped  pulse  generator.  The  Gaussian 
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TABLE  IV. 1 

EXPERIMENTAL  PARAMETER  VALUES 


Primary  Center  Frequency,  600  kHz 

Water  Temperature,  test  period  12-17°C 

Absorption  Coefficient,  600  kHz,  Fresh  Water,  0.0105-0.0140  Np/m 


Speed  of  Sound,  Fresh  Water,  20°C,  cq 
Parametric  Array  Length,  l/2a^c 
Source  Transducer 

Radius 

(cm) 

8.89 

5.08 

2.54 

1.27 

Primary  Pulse  Shape 

Primary  Pulse  Lengths  (0.17  amplitude 
Primary  Pulse  Quality  Factors, 
Secondary  Center  Frequencies,  £^c 
Step-Down  Ratios,  SDR 
Rutherford  Beamwidth  at  66  kHz,  a,  =0 
Primary  Product  Beamwidths 

Radius 

(cm) 

8.89 

5.08 

2.54 

1.27 


1480  m/sec 

35.7-47.6  m 

selectable  radius, 
circular  piston 

Half -Power  Beamwidth 
at  600  kHz  (deg) 

0.75 

1.33 

2.80 

5.90 

Gaussian 

level),  T  14.2-16.6  usee 

8.5-10.0 
66-68  kHz 
~9 

0105  1.40° 


Half -Power  Beamwidth 
_ (deg) _ 

0.53 

0.94 

1.98 

4.17 
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OSCILLATOR 


is  stored  in  PROM,  which  may  be  clocked  at  a  selectable  rate  to  control 
the  pulse  length.  The  digital  representation  of  the  Gaussian  envelope 
is  converted  to  an  analog  representation  and  passed  through  a  low  pass 
filter  before  multiplication  by  the  carrier  signal.  The  Gaussian  pulse 
is  passed  through  a  bandpass  filter  centered  at  the  carrier  to  remove 
unwanted  frequency  components  generated  in  the  shaped  pulse  generator 
and  input  to  a  2000  V*A  power  amplifier.  The  amplified  signal  is 
coupled  to  the  projector  through  impedance  matching  and  tuning  networks 
that  maximize  the  power  transfer  for  each  configuration  of  the  pro¬ 
jecting  array.  The  shaped  pulse  generator  provides  the  master  sync  sig¬ 
nal  on  which  all  timing  is  based.  Additional  details  regarding  the 
transmit  section  and  the  shaped  pulse  generator  are  given  in  Appendix  B. 
The  fixed  column  that  supports  the  projector  is  instrumented  to  provide 
calibrated  electromechanical  train,  tilt,  and  depth  positioning. 

The  transmitting  array  is  the  selectable  radius  circular 
piston  projector  (ARL  384-1)  shown  in  Fig.  4.2.  The  transducer  consists 
of  12  concentric  circular  ceramic  rings  and  a  center  piston,  which  are 
driven  in  the  thickness  mode  and  may  be  selectively  summed  to  form  pis¬ 
tons  of  different  radii.  The  directional  properties  of  the  projector 
at  600  kHz,  which  were  measured  and  predicted  for  a  piston  set  in  an 
infinite  rigid  baffle,  are  summarized  in  Table  IV. 2.  Details  of  the 
transducer  design,  construction,  and  performance  are  discussed  in 
Appendix  B.  Only  4  of  the  13  pistons  were  selected  for  detailed  study 
of  broadband  parametric  arrays. 

The  receive  section  consists  of  two  probe  hydrophones  which 
sampled  the  sound  field.  Their  orientation  was  such  that  they  were 
slightly  displaced  in  range  and  azimuth.  The  outputs  of  these  hydro¬ 
phones  were  input  to  separate  channels  of  passive  filters  and  signal 
amplifiers  and  summed  for  transmission  over  the  same  cable.  The  hydro¬ 
phone  outputs  were  calibrated  at  the  cable  output,  which  was  coupled  to 
another  group  of  filters  and  signal  amplifiers.  The  conditioned  signals 
were  input  to  two  analysis  channels  where  the  signals  that  were  received 
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FIGURE  4.2 

SELECTABLE  RADIUS  CIRCULAR  PISTON  PROJECTOR  (ARL  384-1) 


TABLE  IV. 2 


DIRECTIONAL  PROPERTIES  OF  THE  SELECTABLE 
RADIUS  CIRCULAR  PISTON  PROJECTOR  (ARL  384-1)  AT  600  kHz 


Radius 

(cm) 

Measured  Half- 
Power  Beamwidth 
(deg) 

Theoretical  Half- 
Power  Beamwidth 
(defi) 

*1.27 

5.90 

5.81 

1.91 

3.90 

3.81 

*2.54 

2.80 

2.91 

3.18 

2.14 

2.32 

3.81 

1.18 

1.94 

4.45 

1.44 

1.66 

*5.08 

1.33 

1.45 

5.72 

1.24 

1.29 

6.35 

1.12 

1.16 

6.99 

1.02 

1.06 

7.62 

0.93 

0.97 

8.26 

0.86 

0.89 

*8.89 

0.75 

0.83 

Indicates  pistons  that  were  selected  for  experimental  study 
with  the  broadband  parametric  array. 


at  each  hydrophone  were  separated  by  a  receive  gate.  The  gated  signals 
on  each  channel  were  then  passed  through  a  bandpass  filter  and  coupled 
to  the  analysis  equipment.  The  analysis  equipment  included  an  oscillo¬ 
scope  to  view  the  signal  waveforms,  an  energy  detector  with  both  polar 
and  rectilinear  plotters  to  record  directivity  patterns,  and  a  dual¬ 
channel  FFT  processor.  The  bulk  of  the  analysis  capability  was  provided 
by  the  Nicolet  660A  dual-channel  FFT  processor,  which  was  used  for  both 
time  domain  and  frequency  domain  signal  processing.  Detail  descriptions 
of  the  receive  systems  are  given  in  Appendix  B. 

The  probe  hydrophones  that  were  used  to  measure  the  sound 
field  characteristics  included  the  ARL-239-1,  which  was  used  for  all 
single  channel  measurements  at  both  the  primary  and  secondary  frequen¬ 
cies,  and  the  USRD-E27-130P,  which  was  used  in  the  cross-signal 
analyses.  Each  of  these  hydrophones  was  operated  well  below  resonance 
in  order  to  achieve  reasonably  uniform  frequency  responses  in  the  secon¬ 
dary  band,  which  are  needed  for  direct  comparisons  of  measured  and 
theoretical  results.  In  addition,  these  hydrophones  are  small  so  that 
the  angle  subtended  by  the  hydrophone  cross  section  is  small  and  there¬ 
fore  the  ability  to  determine  subtle  changes  in  signal  characteristics 
as  a  function  of  angle  is  not  compromised. 

A  number  of  steps  were  required  to  reduce  the  experimental 
data  which  is  presented  in  this  chapter.  The  measurement  error  in 
determining  absolute  level  was  estimated  to  be  about  2  dB.  This  meas¬ 
urement  error  results  from  a  number  of  system  and  environmental  con¬ 
siderations  and  is  indicative  of  the  long  term  total  system  stability. 
The  short  term  measurement  error  in  determining  relative  levels  is 
generally  much  less  than  the  long  term  or  absolute  measurement  error  and 
is  estimated  to  be  about  0.5  dB.  There  is  also  error  in  determining  the 
peak  frequency  of  a  broad  smooth  spectral  peak.  This  error  will  be  dis¬ 
cussed  where  applicable  in  the  presentation  of  experimental  results. 


B.  The  Broadband  Parametric  Array  With  Low  Amplitude  Primary  Waves 


The  performance  of  the  broadband  parametric  array  with  low 
amplitude  primary  waves  refers  to  the  operating  conditions  where  the 
behavior  of  the  primary  signal  is  predominately  linear.  This  operating 
domain  is  identified  by  the  harmonic  characteristics  of  the  primary  sig¬ 
nal  as  a  function  of  measurement  range.  There  is  a  somewhat  gradual 
transition  from  linear  to  finite  amplitude  behavior,  and  identification 
of  these  regions  is  discussed  in  Section  C. 

Experimental  data  that  are  presented  in  this  section  are 

compared  with  theoretical  results,  which  are  valid  in  the  farfield  of 

the  broadband  parametric  array.  The  experimental  data  were  obtained  at 

a  measurement  range  of  23.6  m.  This  distance  is  from  50%  to  66%  of  the 

farfield  range,  but  is  also  in  a  region  where  most  of  the  significant 

nonlinear  interaction  has  already  occurred.  Estimates  of  propagation 
37 

curves  for  the  cases  that  are  discussed  indicate  that  data  taken  at 
the  measurement  range  differ  in  amplitude  about  1-1.5  dB  from  farfield 
behavior  of  the  parametric  array;  therefore  theoretical  results  for 
performance  in  the  farfield  of  the  parametric  array  should  be  comparable 
to  experimental  results.  Inputs  to  the  theoretical  models  include  both 
assumed  and  measured  parameters.  The  theoretical  results  are' therefore 
subject  to  errors  that  are  implicit  in  measurement  errors  associated 
with  the  input  parameters. 

1*  High  Primary  Directivity 

The  initial  experimental  effort  involved  a  primary  source 
configuration  that  used  the  full  active  face,  a  *  8.89  cm,  of  the  selec¬ 
table  radius  circular  piston  projector.  In  this  configuration,  the 
primary  half-power  beamwidth  is  0.75°  from  Table  IV. 2,  and  the  Rayleigh 

distance  R  is  10.07  m  or  about  43%  of  the  measurement  range.  This  case 
o 

represents  a  physical  approximation  of  the  broadband  parametric  array 


for  plane  wave  fields  with  a  directivity  ratio  of  about  0.54.  The  water 
temperature  during  the  collection  of  experimental  data  was  about  17°C. 

The  on-axis  acoustic  waveforms  are  shown  in  Fig.  4.3.  The 
narrowband  primary  signal.  Fig.  4.3(a),  is  predominantly  Gaussian 
shaped;  however,  the  rise  time  of  the  waveform  is  slightly  less  than  the 
decay  time.  This  result  is  due  to  response  characteristics  of  the 
tuning  networks  and  the  projecting  array.  The  pulse  length  at  the  17% 
peak  amplitude  level  is  about  17-19  ysec  and  the  measured  peak  sound 
pressure  level  is  195.2  dB  re  1  yPa.  This  level  is  near  the  transition 
from  linear  to  finite  amplitude  behavior  at  the  primary  frequency,  but 
is  still  predominantly  linear.  The  secondary  signal  output,  Fig.  4.3(b), 
has  a  finite  duration  over  about  60  ysec  hut  only  the  first  20-22  ysec 
is  indicative  of  the  secondary  signal  of  Fig.  3.2.  The  remaining  por¬ 
tions  of  the  output  result  from  the  impulse  response  of  the  hydrophone, 
a  situation  which  was  also  noted  by  Ceen  and  Pace.67 

The  off-axis  behavior  of  the  secondary  signal  is  presented  in 

Fig.  4.4.  The  waveforms  are  gated  so  that  only  the  secondary  signal  is 

retained  and  shown  in  0.25°  increments  to  an  off-axis  angle  of  2.75°. 

As  the  sequence  progresses,  the  second  major  positive  peak  in  the  signal 

is  progressively  diminished  and  the  zeroes  and  points  of  inflection  are 

spread  in  time.  This  sequence  indicates  the  trend  described  by 

g 

Westervelt,  where  the  secondary  signal  on-axis  is  proportional  to  the 
second  derivative  of  the  square  of  envelope  function  and,  off-axis,  the 
proportionality  is  with  respect  to  the  first  derivative.  It  should  be 
noted  that  the  display  gain  was  adjusted  to  aid  in  viewing  each  wave¬ 
form,  which  accounts  for  the  increase  in  noise  level  with  increasing 
angle.  An  analysis  of  ungated  secondary  outputs  is  given  in  Appendix  C. 

The  autocorrelation  of  the  secondary  signals  at  each  off-axis 
angle  is  shown  in  Fig.  4.5.  On-axis,  the  autocorrelation  has  three 
positive  and  two  negative  peaks.  As  the  angle  off-axis  increases,  the 
first  and  last  positive  peaks  diminish  and  the  absolute  time  lag  between 
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FIGURE  4.3 

ON-AXIS  ACOUSTIC  WAVEFORMS  FOR 
THE  BROADBAND  PARAMETRIC  ARRAY 
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FIGURE  4.4 

OFF-AXIS  ACOUSTIC  WAVEFORMS  FOR  THE  BROADBAND  PARAMETRIC  ARRAY 
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FIGURE  4.5 

AUTOCORRELATION  OF  THE  OFF-AXIS  ACOUSTIC  SIGNALS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 
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the  two  negative  peaks  increases.  The  relative  amplitude  of  the 
autocorrelation  at  each  angle  measured  at  t=0  is  noted  in  the  figure. 
These  relative  amplitudes  trace  the  energy  beam  pattern  for  the 
radiation. 

A  more  detailed  analysis  of  the  properties  of  the  broadband 
parametric  array  may  be  considered  in  the  frequency  domain.  The  nor¬ 
malized  power  density  spectra  for  the  measured  on-axis  transmission  are 
shown  in  Fig.  4.6.  This  spectrum  has  250  Hz  resolution  and  was  computed 
from  the  incoherent  average  of  100  individual  spectra  ove.  successive 
transmission  sequences  with  a  repetition  rate  of  about  4  msec.  The 
spectrum  was  calibrated  at  65  kHz  by  substituting  a  pure  tone  at  the 
calibrated  receive  system  output  and  adjusting  the  amplitude  to  produce 
an  equivalent  signal  input.  The  result  was  adjusted  to  compensate  for 
the  resolution  bandwidth  of  the  FFT  processor,  and  the  relative  duration 
of  the  measured  signal  and  the  calibration  signal.  Using  this  proce¬ 
dure,  the  peak  sound  pressure  level,  SPI^,  of  the  66  kHz  spectral  compo¬ 
nent  was  determined  to  be  158.0  dB  re  1  pPa  at  23.6  m. 

Theoretical  results  were  computed  for  comparison  with  measured 
results  by  numerical  integration  of  Eq.  (3.45).  This  equation  is  the 
solution  for  the  frequency  components  of  the  secondary  pressure  in  the 
farfield  of  the  broadband  parametric  array  assuming  a  Gaussian  pulse 
primary  transmission  from  a  circular  piston.  Since  the  primary  trans¬ 
mission  is  not  exactly  Gaussian,  the  specific  characteristics  of  the 
transmission  were  established  from  both  the  experimental  data  and  itera¬ 
tions  of  Eq.  (3.45).  The  parameters  of  the  primary  transmission  are 
given  by  a  of  10  and  a  pulse  duration  of  16.7  usee  at  the  171  peak 
amplitude  level. 

With  these  parameter  values,  an  on-axis  peak  secondary  sound 
pressure  level  at  66  kHz  of  162.5  dB  re  1  uPa  at  23.6  m  is  predicted 
from  Eq.  (3.45),  which  is  about  4.5  dB  greater  than  the  measured  value. 
About  2  dB  of  this  difference  results  from  the  incompleteness  of  the 


parametric  array  with  the  measurement  range  at  the  end  of  the  nearfield. 
The  remaining  2.5  dB  is  due  to  the  cumulative  effects  of  errors  in  meas¬ 
ured  values  and  assumed  values.  The  normalized  power  spectral  density 
function  was  computed  from  Eqs.  (3.21)  and  (3.45)  and  is  shown  in 
Fig.  4.6  with  the  experimental  results.  Comparison  of  the  two  curves 
indicates  excellent  agreement  between  experimental  and  theoretical 
results  above  the  -15  dB  power  level.  Below  this  level  at  low  fre¬ 
quencies,  the  curves  slightly  diverge  with  more  power  being  present  in 
the  experimental  data.  This  behavior  results  from  several  considera¬ 
tions  including  the  fact  that  the  pulse  shown  in  Fig.  4.3  is  not  exactly 
Gaussian  as  assumed  in  the  theoretical  development.  Also,  the  effect  of 
the  frequency  response  of  the  receiving  hydrophone  was  greatly  reduced 
(but  not  totally  eliminated)  by  gating  the  receive  signal  as  discussed 
in  Appendix  C.  Finally,  as  the  S/N  declines  at  lower  frequencies,  the 
influence  of  the  background  noise  level  is  increased. 

The  normalized  power  density  spectra  for  the  secondary 
transmissions  at  off-axis  angles  from  0.0°  to  3.0°  are  presented  at 
angular  increments  of  0.25°  in  Figs.  4.7  and  4.8.  As  the  angle  off-axis 
is  increased,  the  level  and  the  frequency  at  the  spectral  peak  are 
reduced,  which  is  the  trend  predicted  for  the  broadband  parametric  array 
for  plane  wave  fields  and  which  is  under  investigation  for  the  broadband 
parametric  array  for  nonplane  wave  fields.  Quantitative  comparisons  of 
the  normalized  power  spectra  for  the  measured  and  predicted  data  are 
awkward  in  this  presentation  format.  A  format  better  suited  for  such 
comparisons  is  presented  in  Fig.  4.9,  where  the  frequency  and  level  at 
the  spectral  peak  are  plotted  as  a  function  of  normalized  off-axis 
angle.  The  normalization  is  with  respect  to  the  measured  half-power 
angle  8^  .  This  data  presentation  removes  the  effects  of  absolute 

angle  and  preserves  the  differences  between  data  that  are  dominantly 
absorption  or  diffraction  limited.  Each  of  the  power  spectra  is  char¬ 
acterized  by  a  broad  peak.  In  determining  the  frequency  and  level  at 
the  spectral  peak  from  the  experimental  measurements,  the  peaks  were 
assumed  to  be  symmetric  above  the  -1.5  dB  normalized  level.  This 
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RELATIVE  LEVEL  AT  THE  FREQUENCY  AT  THE 

SPECTRAL  PEAK  -  dB  SPECTRAL  PEAK  -  kHz 


FIGURE  4.9 

FREQUENCY  AND  LEVEL  AT  THE  SPECTRAL  PEAK  OFF-AXIS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

f(t)  -  GAUSSIAN,  f,c  =  600  kHz,  a  =  8.89  cm,  R0  =  10.07  m,  r  =  23.6  m 
03dB  =0-70°  (MEASURED) 
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assumption  leads  to  an  uncertainty  of  about  2%  in  determining  the 
frequency  at  the  spectral  peak  and  has  a  negligible  effect  on  the  deter¬ 
mination  of  level.  For  the  data  shown  in  Fig.  4.9,  the  power  spectrum 
on-axis  has  a  peak  at  about  66  kHz  and  a  half -power  width  of  about 
60  kHz.  The  theoretical  result/  shown  in  Fig.  4.9  are  computed  from  the 
farfield  solution  given  by  numerical  integration  of  Eq.  (3.45)  and 
several  limiting  forms.  The  limiting  forms  include  the  absorption 
limit,  Eq.  (3.20),  and  the  diffraction  limit  for  a  circular  piston, 

Eq.  (3.46),  and  the  Gaussian  approximation  to  a  circular  piston, 

Eq.  (3.52). 


Experimental  data  are  shown  in  Fig.  4.9  for  four  different 

primary  peak  sound  pressure  levels,  SPL^,  at  23.6  m.  These  data  are 

predominantly  low  level  in  the  primary  transmission.  As  the  angle  off- 

axis  is  increased,  the  trend  toward  decreasing  frequency  at  the  spectral 

peak  is  observed.  The  trend  is  characterized  by  a  relatively  constant 

frequency  for  the  spectral  peak  at  small  angles  and  transitions  to  a 

steady  downward  shift  in  frequency  for  the  spectral  peak  at  angles 

greater  than  about  0,  .  The  overall  shift  was  measured  to  be  from 

J  do 

about  66  kHz  on-axis  to  about  41  kHz  at  an  angle  of  4.6(9,  )  or  2.8°. 

j  do 

The  limitation  on  the  acquisition  of  data  at  larger  angles  was  due  to 
deteriorating  S/N.  It  should  be  noted  that  S/N's  at  off-axis  angles 
greater  than  3(0g  dg  )  were  extremely  poor  for  the  data  taken  with  an 
SPL^  of  187.4  dB.  As  a  result,  the  low  frequency  content  of  the  nor¬ 
malized  power  density  spectrum  may  have  been  contaminated  by  ambient 
noise  in  the  medium,  which  would  contribute  to  lower  indications  of  fre¬ 
quency  at  the  spectral  peak  for  these  data.  The  relative  level  at  the 
spectral  peak  is  indicative  of  the  energy  beam  pattern.  These  data  are 
monotonically  decreasing  to  a  level  of  -20  dB  and  yield  a  half-power 
beamwidth  of  about  1.4°  for  the  secondary  transmission. 

The  data  are  dominantly  linear  with  the  exception  of  the  curve 
taken  at  the  highest  on-axis  peak  sound  pressure  level,  SPL^.  For  these 
data,  the  trend  toward  downward  shift  in  frequency  at  the  spectral  peak 
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is  slower  and  indicates  the  onset  of  finite  amplitude  effects  in  the 
primary  transmission.  This  behavior  is  also  noted  in  the  relative  level 
at  the  spectral  peak;  however,  the  trend  is  not  as  pronounced.  These 
results  indicate  that  the  frequency  content  of  the  secondary  signal  is  a 
more  sensitive  indicator  of  the  onset  of  finite  amplitude  effects  in  the 
primary  transmission  than  the  relative  level.  The  performance  of  the 
broadband  parametric  array  during  transition  from  dominantly  linear 
behavior  in  the  primary  transmission  to  strong  finite  amplitude  behav¬ 
ior  is  discussed  further  in  Section  IV. B. 

Theoretical  results  are  presented  in  Fig.  4.9  for  comparison 
with  experimental  data.  In  the  absorption  limit  (from  Eq.  (3.20)),  the 
peak  frequency  on-axis  is  about  6  kHz  higher  than  observed.  The  trend 
toward  a  downward  shift  in  frequency  at  the  spectral  peak  is  slow  for 
small  angles,  increases  to  a  rapid  shift  in  the  region  from  0.5(0^  ^ ) 
to  1.5(6_  ),  and  becomes  relatively  constant  at  large  angles.  The 

J  UD 

level  at  the  spectral  peak  in  the  absorption  limit  is  monotonically 
decreasing  and  is  slightly  lower  in  level  than  the  experimental  data. 

In  the  diffraction  limit  (from  Eq.  (3.46)),  the  peak  on-axis 
frequency  is  about  15  kHz  lower  than  observed.  The  frequency  at  the 
spectral  peak  is  constant  to  about  0,  ,  where  a  rapid  decrease  occurs. 

At  larger  angles,  lobes  in  the  frequency  at  the  spectral  peak  are  noted, 
which  result  from  lobes  in  the  off-axis  transfer  function  shown  in 
Fig.  3.8.  No  theoretical  data  are  shown  between  lobes  in  the  diffrac¬ 
tion  limit.  In  these  regions  the  frequency  at  the  spectral  peak  is  not 
well  behaved;  however,  these  regions  are  in  the  nulls  of  the  off-axis 
transfer  function  and  little  energy  is  contained  in  the  frequency  com¬ 
ponents.  The  level  at  the  spectral  peak  decreases  very  rapidly  and 
indicates  a  half-power  beamwidth  of  about  0.3°.  In  the  diffraction 
limit  with  a  Gaussian  approximation  (from  Eq.  (3.52))  for  the  off-axis 
transfer  function,  the  frequency  at  the  spectral  peak  has  the  same  on- 
axis  value  as  observed  with  the  diffraction  limit  for  the  circular  pis¬ 
ton  transfer  funct?  m.  At  angles  off-axis,  the  frequency  behavior 
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differs  from  that  of  Eq.  (3.46).  The  frequency  at  the  spectral  peak 

decreases  monotonically  with  a  3.5  kHz  difference  at  34.  .  In  con- 

1  do 

trast,  the  relative  level  at  the  spectral  peak  for  both  diffraction 
limited  representations  are  in  agreement  to  the  -30  dB  level.  At  some 
level  below  -30  dB  the  curves  diverge  as  the  sidelobe  structure  in  the 
off-axis  transfer  function  for  the  circular  piston  becomes  dominant. 

The  limiting  forms  of  the  broadband  parametric  array  for 
nonplane  wave  fields  do  not  fit  the  experimental  data;  but  these 
limiting  forms  do  bracket  the  trends  that  are  noted  in  the  data.  The 
relative  values  of  the  on-axis  frequency  at  the  spectral  peak  may  be 
evaluated  in  determining  the  exponent  of  the  overall  dependence  of  the 
secondary  frequency.  Recall  that  the  dependence  is  quadratic  in  the 
absorption  limit  and  linear  with  secondary  frequency  in  the  diffraction 
limit.  The  overall  dependence  of  the  experimental  results  on  secondary 
frequency  is  given  by  w*.88 

The  solution  in  the  farfield  of  the  broadband  parametric  array 
with  a  Gaussian  pulse  primary  transmission  from  a  circular  piston  is 
given  by  Eq.  (3.45),  which  was  evaluated  numerically  for  comparison  with 
the  experimental  data.  These  theoretical  results  indicate  the  predicted 
frequency  at  the  spectral  peak  on-axis  and  are  in  excellent  agreement 
with  the  experimental  data  for  angles  off-axis  to  0^  •  At  larger 

angles,  the  theoretical  results  are  in  general  agreement  but  predict 
slightly  lower  frequencies  at  the  spectral  peak.  Off-axis,  the  char¬ 
acteristics  of  the  broadband  parametric  array  for  nonplane  wave  fields 
are  strongly  influenced  by  the  sidelobe  structure  in  the  off-axis  trans¬ 
fer  function.  The  difference  between  the  theoretical  predictions  and 
the  measured  results  is  due  to  the  difference  in  the  assumed  off-axis 
transfer  or  directivity  function  for  a  piston  in  an  infinite  rigid  baf¬ 
fle  and  the  real  off-axis  transfer  function  measured  for  the  source  and 
shown  in  Appendix  A.  The  measured  off-axis  transfer  or  directivity 
functions  are  complex  in  sidelobe  structure.  The  sidelobes  are  gen¬ 
erally  4-5  dB  higher  than  predicted  for  a  piston  in  an  infinite  rigid 


baffle  and  the  presence  of  spacing  lobes,  which  result  from  the 
particular  design  of  the  source  transducer,  are  noted.  The  higher  side- 
lobes  for  the  source  tend  to  retard  the  downward  shift  in  frequency  at 
the  spectral  peak  in  the  measured  data.  This  trend  will  be  discussed 
further  in  relation  to  the  secondary  transmissions  produced  by  primary 
radiation  from  other  source  configurations.  For  the  parameters  of  the 
broadband  parametric  array  considered  in  £!■  experiment,  the  convolu¬ 
tion  integral  is  dominated  by  Rutherford  .  •' Bering.  This  dominance 

tends  to  reduce  and  average  the  effete**  *•*  sidelobes  in  the  primary 
off-axis  transfer  function.  At  large  .>n..  the  influence  of  the  pri¬ 
mary  transfer  function  shifts  fron-\  nu.  u  lobe  to  the  sidelobes,  and 
is  noted  in  the  theoretical  results  from  Eq.  (3.45). 

The  relative  level  at  the  spectral  peak  is  in  excellent 
agreement  with  the  theoretical  results  predicted  by  Eq.  (3.45).  At  the 
half-power  angle,  the  theoretical  results  in  the  absorption  limit  and 
from  Eq.  (3.45)  are  in  close  agreement  with  the  experimental  data.  This 
agreement  was  noted  in  earlier  studies  as  verification  of  the  Westervelt 
model  for  the  parametric  array  for  plane  wave  fields.  This  agreement  is 
good  at  0^  jg  ;  however,  at  larger  angles  the  normalized  level  in  the 
absorption  limit  has  a  more  rapid  rate  of  reduction.  This  result  sup¬ 
ports  the  conclusion  that  the  parametric  array  formed  with  the  given 
set  of  experimental  parameters  is  only  approximated  by  the  analytical 
description  of  the  broadband  parametric  array  for  plane  wave  fields. 

A  power  beam  pattern  for  the  secondary  transmission  is  shown 
in  Fig.  4.10.  This  pattern  has  a  half-power  width  of  about  1.2°  and  the 
characteristic  appearance  of  the  parametric  directivity  with  no  discern¬ 
ible  sidelobes  over  about  25  dB  of  S/N.  This  measurement  was  performed 
with  an  energy  detector.  A  comparison  of  power  beam  pattern  representa¬ 
tions  is  shown  in  Fig.  4.11.  In  this  figure,  normalized  levels  are 
plotted  as  a  function  of  angle  for  the  output  of  an  energy  detector,  the 
amplitude  of  the  autocorrelation  for  time  shift  equal  to  zero,  and  the 
relative  level  at  the  spectral  peaks.  The  close  agreement  of  the  data 
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FIGURE  4.10 

SECONDARY  POWER  BEAM  PATTERN  FOR 
THE  BROADBAND  PARAMETRIC  ARRAY 

f(t)  -  GAUSSIAN.  f,c  =  600  kHz,  SPL,  =  195.2  dB  re  1  fiPa 
f2c  66  kHz,  a  =  8.89  cm,  RQ  =  10.07  m,  r  =  23.6  m 
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FIGURE  4.11 

COMPARISON  OF  SECONDARY  POWER  BEAM  PATTERN  REPRESENTATIONS 

f(t)  -  GAUSSIAN,  f,e  =  600  kHz,  SPL,  =  195.2  dB  re  1  *zPa,  f,c  S  66  kHz 
a  =  8.89  cm,  RQ  =  10.07  m,  r  =  23.7  m,  03  dB  =  0.70* 
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indicates  that  each  of  these  methods  is  indicative  of  the  total  energy 
in  the  secondary  signal  for  the  parametric  transmissions  and  may  be 
used  to  trace  the  power  beam  pattern. 

2.  Performance  as  a  Function  of  Primary  Directivity 

As  the  primary  directivity  is  decreased,  a  transition  in  the 
characteristics  of  the  broadband  parametric  array  occurs  from  that  of 
plane  wave  fields  in  the  absorption  limit  (physical  approximation)  to 
that  of  nonplane  wave  fields  in  the  diffraction  limit.  In  conjunction 
with  this  trend,  the  dominance  of  the  off-axis  transfer  function  shifts 
from  that  associated  with  Rutherford  scattering  to  the  square  of  the 
off-axis  primary  transfer  function.  The  off-axis  characteristics  of 
the  secondary  field  as  a  function  of  primary  directivity  were  experi¬ 
mentally  evaluated  for  directivity  ratios  of  0.89,  1.74,  and  3.69. 

These  experimental  results  are  representative  of  intermediate  direc¬ 
tivity  for  the  broadband  parametric  array  and  may  be  compared  with  the 
results  from  the  preceding  section  where  a  physical  approximation  for 
the  broadband  parametric  array  in  the  absorption  limit  was  formed  with 
a  directivity  ratio  of  0.54. 

The  broadband  parametric  array  with  a  directivity  ratio  of 
0.89  was  formed  with  the  selectable  radius  circular  piston  set  to  a 
radius  of  5.08  cm  and  evaluated  in  water  at  a  temperature  of  about  15°C. 
The  parameters  of  the  primary  transmission  are  given  by  a  of  9.8,  a 
pulse  duration  of  16.3  ysec,  and  an  on-axis  primary  peak  sound  pressure 
level,  SPL^,  of  189.0  dB  re  1  yPa  at  the  measurement  range.  The  fre¬ 
quency  and  level  at  the  spectral  peak  were  extracted  from  measured  power 
density  spectra  and  are  shown  in  Fig.  4.12  with  theoretical  results  pre¬ 
sented  for  comparison.  An  on-axis  secondary  peak  sound  pressure  level, 
SPL^,  of  147.9  dB  re  1  yPa  at  23.6  m  was  measured  at  66  kHz.  This 
measurement  is  about  4  dB  lower  than  the  value  of  152.0  dB  re  1  uPa 
predicted  from  Eq.  (3,45).  The  difference  is  attributed  to  the  incom¬ 
pleteness  of  the  parametric  array  and  the  cumulative  effects  of  errors 
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FIGURE  4.12 

FREQUENCY  AND  LEVEL  AT  THE  SPECTRAL  PEAK  OFF-AXIS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

f(t)  -  GAUSSIAN,  f1c  =  600  kHz,  SPL,  =  189.0  dB  re  1  juPa 
a  =  5.08  cm,  RQ  =  3.29  m,  r  =  23.7  m 
6 2  dg  =0.80°  (MEASURED) 
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in  measured  and  assumed  values.  The  trends  in  these  data  are  similar  to 
those  for  a  directivity  ratio  of  0.54;  however,  the  rate  of  downward 
shift  in  frequency  at  the  spectral  peak  with  angle  off-axis  is  reduced. 
The  frequency  shift  with  angle  is  not  monotonically  decreasing  as  indi¬ 
cated  by  a  minor  frequency  peak  at  about  5.3(02^).  The  relative  level 
at  the  spectral  peak  is  monotonically  decreasing  with  no  apparent  minor 
peaks  or  sidelobes. 

The  theoretical  results  for  the  limiting  forms  of  the 
broadband  parametric  array  are  shown  for  the  absorption  limit,  the 
diffraction  limit,  and  the  Gaussian  approximation  in  the  diffraction 
limit.  These  limiting  forms  bracket  the  trends  in  the  experimental 
data.  The  relative  values  of  the  on-axis  peak  spectral  frequencies 
indicate  that  the  overall  dependence  of  the  experimental  data  on  secon¬ 
dary  frequency  is  indicating  a  slight  shift  toward  the  diffrac¬ 

tion  limit  relative  to  the  data  presented  for  a  directivity  ratio  of 
0.54.  The  theoretical  results  from  the  model,  Eq.  (3.45),  are  presented 
in  the  figure.  Good  agreement  between  the  experimental  data  and  the 
theoretical  results  is  noted  in  the  frequency  at  the  spectral  peak  for 
off-axis  angles  to  9.  .  For  larger  angles,  the  theoretical  results 

indicate  frequencies  at  the  spectral  peak  that  are  slighcly  2 ower  than 
observed.  This  is  the  same  trend  noted  in  previous  analyses,  which  is 
attributed  to  the  higher  sidelobe  levels  in  the  real  piimary  off-axis 
transfer  functions  than  assumed  in  the  theoretical  model.  Both 
Eq.  (3.45)  and  the  experimental  results  indicate  minor  frequency  peaks 
that  result  from  the  cumulative  effects  of  secondary  lobes  in  the  pri¬ 
mary  off-axis  transfer  function.  The  minor  frequency  peaks  do  not  occur 
at  the  same  angle,  which  is  another  indication  that  the  real  and  the 
assumed  primary  transfer  functions  differ. 

Close  agreement  over  all  angles  is  noted  in  the  level  at  the 
spectral  peak,  with  the  results  from  Eq.  (3.45)  only  slightly  higher 
than  the  measured  data.  There  is  no  apparent  influence  of  the  secondary 

lobes  of  the  primary  transfer  function  on  the  level  at  the  spectral  peak 
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over  a  25  dB  range.  The  energy  beam  pattern  for  this  parametric  array 
is  shown  in  Fig.  4.13.  The  measured  half -power  beamwidth  is  about  1.6°, 
and  there  are  no  sidelobes  apparent  in  the  pattern. 


The  broadband  parametric  array  with  a  directivity  ratio  of 
1.74  was  formed  with  a  radius  of  2.54  cm  and  evaluated  in  water  at  about 
12°C.  The  parameters  of  the  primary  transmission  are  given  by  a  of 
8.9,  a  pulse  duration  of  14.9  psec ,  and  an  on-axis  primary  peak  sound 
pressure  level  of  187.8  dB  re  1  pPa  at  the  measurement  range.  The  fre¬ 
quency  and  level  at  the  spectral  peak  extracted  from  both  measured  and 
theoretical  results  are  shown  in  Fig.  4,14.  The  on-axis  secondary  peak 
sound  pressure  level,  SPL^,  is  153.6  dB  re  1  uPa  at  23.6  m  measured  at 
67  kHz.  The  predicted  value  for  the  secondary  sound  pressure  level  com¬ 
puted  from  Eq.  (3.45)  is  156.5  dB  re  1  yPa  at  23.6  m.  The  predicted 
value  is  higher  than  the  measured  value  by  about  3  dB,  which  again  is 
attributed  to  the  incompleteness  of  the  parametric  array  and  errors  in 
measured  and  assumed  values.  A  further  decrease  in  the  rate  of  downward 
shift  in  frequency  as  a  function  of  off-axis  angle  is  indicated  relative 
to  the  results  for  smaller  directivity  ratios,  and  a  minor  frequency 
peak  is  noted  at  about  4.4(0^^).  The  relative  level  at  the  spectral 
peak  is  monotonically  decreasing  and  there  is  some  indication  of  a 
secondary  peak  at  about  4.3(0,  ).  It  is  not  surprising  that  sidelobes 

should  appear  in  the  secondary  beam  patterns  as  the  diffraction  limit  is 

approached.  Sidelobes  have  been  observed  in  the  beam  patterns  of  the 
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narrowband  parametric  array  for  high  directivity  ratios. 

The  theoretical  results  for  the  limiting  forms  for  the 

broadband  parametric  array  are  shown  in  Fig.  4.14.  The  transition 

toward  the  diffraction  limit  is  noted  in  the  overall  dependence  on 

1.37 

secondary  frequency  given  by  and  the  increased  influence  of  the 

primary  off-axis  transfer  functions.  In  comparing  the  experimental 

data  with  the  results  from  Eq.  (3.45),  the  behavior  noted  is  similar  to 

the  cases  for  lower  directivity  ratios.  Good  agreement  is  observed  in 

the  frequency  at  the  spectral  peak  for  off-axis  angles  to  0,  .  At 

j  do 
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FIGURE  4.13 

SECONDARY  POWER  BEAM  PATTERN  FOR 
THE  BROADBAND  PARAMETRIC  ARRAY 

f(t)  -  GAUSSIAN,  f1c  =  600  kHz,  SPL,  =  189.0  dB  re  1  juPa 
*2c  ~  66  kHz,  a  =  6.08  cm'  R0  =  ^-29  m'  r  =  23.6  m 
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larger  angles  the  theoretical  results  are  slightly  low  because  the 
sidelobe  levels  in  the  real  off-axis  transfer  functions  are  higher  than 
assumed  in  the  model.  In  these  data,  however,  the  minor  frequency  peaks 
occurring  in  both  the  experimental  and  theoretical  data  are  approxi¬ 
mately  aligned  over  off-axis  angles  from  about  3.5(8-  )  to  5(0,  ,D). 

J  dB  3  dB 

This  behavior  is  observed  even  though  the  sidelobe  levels  are  different. 
Fewer  sidelobes  in  the  primary  off-axis  transfer  function  are  combined 
to  give  the  general  result.  As  the  directivity  of  a  linear  source  is 
decreased,  the  lower  order  sidelobes  become  broader  and  agreement 
between  theoretical  and  actual  orientation  of  the  lobes  generally 
improves . 


The  agreement  between  the  theoretical  and  experimental  results 
is  excellent  except  where  the  sidelobe  is  noted  in  the  experimental 
data.  No  discernible  lobe  is  predicted  in  this  region  from  Eq.  (3.45). 
The  energy  beam  pattern  is  shown  in  Fig.  4.15  where  a  half-power  beam- 
width  of  about  2.9°  is  indicated. 

The  final  broadband  parametric  array  was  formed  with  a 
directivity  ratio  of  3.69  and  a  piston  radius  of  1.27  cm  and  was  eval¬ 
uated  in  water  at  about  12°C.  The  parameters  of  the  primary  transmis¬ 
sion  are  given  by  a  of  8.5,  a  pulse  duration  of  14.2  psec,  and  an 
on-axis  primary  peak  sound  pressure  level  of  182.0  dB  re  1  pPa  at  the 
measurement  range.  The  measured  and  theoretical  results  for  the  fre¬ 
quency  and  level  at  the  spectral  peak  are  shown  in  Fig.  4.16.  The  on- 
axis  secondary  peak  sound  pressure  level,  SPI^,  for  this  case  is 
149.6  dB  re  1  pPa  at  23.6  m  measured  at  68  kHz.  The  secondary  sound 
pressure  level  computed  from  Eq.  (3.45)  is  predicted  to  be 
151.0  dB  re  1  pPa  at  the  measurement  range.  The  predicted  value  differs 
from  the  measured  value  by  about  1.4  dB  due  to  errors  in  measured  and 
assumed  values  and  to  the  incompleteness  of  the  parametric  array.  The 
rate  of  downward  shift  in  frequency  has  further  decreased  with  a  broad 
minor  peak  in  frequency  at  about  4.3(02^).  The  relative  level  at  the 
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spectral  peak  is  monotonically  decreasing  with  a  broad,  yet  well  defined 

sidelobe  indicated  near  4,0(9,  ._  ). 

3  do 

Analyses  of  the  limiting  forms  for  the  broadband  parametric 

array  indicate  a  further  transition  toward  the  diffraction  limit.  The 

overall  dependence  on  secondary  frequency  is  determined  from  the  fre- 

1  3 

quency  at  the  spectral  peak  and  is  given  by  '  •  The  lobes  in  fre¬ 
quency  at  the  spectral  peak  correspond  very  well  with  the  lobes  in  the 
diffraction  limit.  The  level  at  the  spectral  peak  is  also  better  esti¬ 
mated  in  the  diffraction  limit  than  in  the  absorption  limit. 

The  experimental  data  are,  as  previously  noted,  higher  in 
frequency  at  the  spectral  peak  than  the  theoretical  results  from 
Eq.  (3.45)  for  off-axis  angles  greater  than  0^  ^g  •  However,  the  broad 
peak  observed  in  the  data  is  predicted  very  well  by  the  theoretical 
results.  As  the  diffraction  limit  is  approached,  the  Rutherford  direc¬ 
tivity  acts  as  a  sampling  function  in  the  convolution  integral  and  the 
product  of  the  primary  off-axis  transfer  functions  is  indicated. 

The  theoretical  and  experimental  results  are  in  close 
agreement  for  the  relative  level  at  the  spectral  peak  to  about  -20  dB  or 
an  off-axis  angle  of  about  3.0(9^^).  At  lower  levels  and  larger 
angles,  both  the  experimental  data  and  the  theoretical  results  from 
Eq.  (3.45)  indicate  a  broad  sidelobe;  however,  the  level  is  several 
decibels  higher  for  the  measured  data.  The  half -power  beamwidth  for 
this  broadband  parametric  array  is  about  5.3°  as  measured  from  the 
energy  beam  pattern  shown  in  Fig.  4.17. 

A  composite  of  the  frequency  and  level  at  the  spectral  peak 
for  various  directivity  ratios  or  source  piston  radii  is  given  in 
Fig.  4.18.  Data  are  included  for  each  directivity  ratio.  The  presenta¬ 
tion  of  the  data  as  a  function  of  angle  normalized  to  the  half-power 
angle  for  each  directivity  ratio  has  removed  the  effects  of  absolute 
angle  and  allows  the  characteristics  of  the  off-axis  transfer  function 
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FIGURE  4.18 

FREQUENCY  AND  LEVEL  AT  THE  SPECTRAL  PEAK  OFF-AXIS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 
AND  SEVERAL  SOURCE  PISTON  RADII 

f(t)  -  GAUSSIAN,  f,c  =  600  kHz,  r  =  23,6  m 
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to  be  studied  as  the  transition  from  the  absorption  limit  to  the 
diffraction  limit  progresses.  The  trend  of  decreasing  frequency  at  the 
spectral  peak  with  increasing  angle  off-axis  is  noted  in  each  case;  how¬ 
ever,  the  rate  of  frequency  shift  is  greater  for  smaller  directivity 
ratios  at  off-axis  normalized  angles  greater  than  about  2.5(02  ^g)-  At 
larger  normalized  angles,  the  trends  diverge  as  the  off-axis  transfer 
function  for  each  directivity  ratio  becomes  dominant.  As  the  transition 
from  the  absorption  limit  (small  directivity  ratios)  to  the  diffraction 
limit  (larger  directivity  ratios)  occurs,  the  trend  is  toward  frequency 
independence  as  a  function  of  off-axis  angle.  The  level  at  the  spectral 
peak  is  reasonably  well  behaved  at  all  angles  and  has  the  same  trend  for 
all  directivity  ratios  to  a  normalized  angle  of  about  1 . 6 ( 9  ) .  The 

relative  behavior  of  the  frequency  at  the  spectral  peak  and  the  level  at 
the  spectral  peak  indicates  the  higher  sensitivity  of  frequency  in  char¬ 
acterizing  the  secondary  sound  field. 

A  more  striking  comparison  of  data  for  various  directivity 
ratios  derived  in  part  from  different  source  piston  radii  is  given  in 
Fig.  4.19,  where  the  frequency  at  the  spectral  peak  is  given  as  a  func¬ 
tion  of  absolute  off-axis  angle,  and  in  Fig.  4.20,  where  the  level  at 
the  spectral  peak  is  given  as  a  function  of  absolute  off-axis  angle. 
Prominent  minor  peaks  in  frequency  at  the  spectral  peak  occur  at  off- 
axis  angles  of  about  4°,  6°,  and  11°  for  broadband  parametric  arrays 
formed  by  source  piston  radii  of  5.08  cm,  2.54  cm,  and  1.27  cm, 
respectively.  These  features  are  correlated  with  indications  of  low 
sidelobes  in  the  relative  level  at  the  spectral  peak  for  the  parametric 
arrays  formed  with  the  two  smallest  source  piston  radii  and  highest 
directivity  ratios. 

A  final  comparison  of  experimental  data  and  theoretical 
predictions  as  a  function  of  primary  directivity  is  shown  in  Fig.  4.21, 
where  the  secondary  half -power  beamwidth  is  given  as  a  function  of  pri¬ 
mary  half-power  beamwidth,  directivity  ratio,  and  source  piston  radius. 
The  Rutherford  beamwidth  is  indicated  from  1.4°  to  1.6°  due  to  the 
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FIGURE  4.21 

RELATIONSHIP  BETWEEN  PRIMARY  AND  SECONDARY  DIRECTIVITY  AND 
SOURCE  PISTON  RADIUS  FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

f(t)  -  GAUSSIAN,  f1c  =  600  kHz,  f2c  3s  66  kHz,  r  =  23.6  m 
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variations  in  absorption  imposed  by  temperature  changes  in  the  medium 
during  the  acquisition  of  experimental  data.  The  directivity  ratio  was 
computed  based  on  the  average  Rutherford  beamwidth  and  the  primary  pro¬ 
duct  beamwidth  was  computed  based  on  the  experimentally  determined  beam- 
widths.  The  agreement  between  measured  and  predicted  results  is  good. 
For  very  narrow  primary  beams  the  directivity  ratio  is  small  and  spatial 
performance  is  asymptotic  to  the  Rutherford  directivity,  which  is  pre¬ 
dicted  in  the  absorption  limit.  This  case  is  approximated  physically 
with  a  large  aperture  size  where  most  of  the  nonlinear  interaction 
occurs  in  the  Fresnel  diffraction  region  of  the  primary  source.  As  the 
primary  directivity  and  the  directivity  ratio  decrease,  the  nonlinear 
interaction  that  takes  place  in  the  Fraunhofer  diffraction  region 
increases  and  a  transition  occurs  in  the  spatial  characteristics  to  an 
asymptote  determined  by  the  primary  product  directivity,  which  is  pre¬ 
dicted  in  the  diffraction  limit.  This  behavior  is  the  same  as  noted 

37 

for  narrowband  parametric  arrays  with  intermediate  directivity. 

3.  Cross  Function  Analysis 

An  alternate  method  for  evaluating  the  off-axis  properties  of 
the  broadband  parametric  transmission  is  the  use  of  cross  functions. 

This  method  utilizes  coherent  signal  analysis  in  determining  relative 
signal  properties.  Experimental  results  were  obtained  by  evaluating 
cross  function  properties  of  the  on- .xis  signal  with  various  off-axis 
signals.  The  ARL  239-1  was  used  to  receive  the  on-axis  signals  and  the 
USRD-E27-130P  was  used  to  receive  the  off-axis  signals.  The  hydrophones 
were  oriented  at  various  angles  with  respect  to  the  source  and  separated 
in  range  from  about  0.3  m  to  0.6  m.  This  separation  in  range  allowed 
the  signals  to  be  transmitted  over  the  same  cable  with  a  150-350  nsec 
difference  in  arrival  time.  The  signals  were  separated  by  gates  and 
coupled  to  the  dual-channel  analyzer  where  the  time  delay  between  sig¬ 
nals  was  artificially  corrected.  Cross  functions  were  evaluated  by 
setting  the  time  delay  so  that  the  crosscorrelation  was  a  maximum  of 
time  shift,  t=0. 
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The  selectable  radius  circular  piston  projector  was  driven 
over  the  full  aperture  with  a  radius  of  8.89  cm,  and  at  a  primary  peak 
sound  pressure  level  of  195.2  dB  re  1  pPa  at  23.6  in.  The  crosscorrela¬ 
tion  of  the  on-axis  and  off-axis  signals  at  1.8°  (36.  ,  )  and  3°  (50. 

J  dB  J  dB 

are  shown  in  Fig.  4.22.  These  results  may  be  compared  with  the  auto¬ 
correlation  on-axis  in  Fig.  4.5.  The  duration  of  the  crosscorrelation 
signals  is  about  36  psec.  As  the  angle  is  increased  the  crosscorrela¬ 
tion  signals  become  asymmetric,  and  the  first  positive  and  negative 
peaks  diminish  in  amplitude  relative  to  the  other  peaks  in  the  signal. 

The  changes  that  occur  in  the  crosscorrelation  signals  are 
more  readily  evaluated  in  the  frequency  domain  with  the  crosspower  den¬ 
sity  spectrum.  Normalized  crosspower  amplitude  spectra  for  the  on-axis 
and  off-axis  signals  are  shown  in  Fig.  4.23.  These  curves  result  from 
the  product  of  two  secondary  signal  spectra  that  have  different  peak 
frequencies.  Since  these  are  low  Q2  transmissions  and  the  peak  fre¬ 
quencies  are  relatively  close,  the  resultant  crosspower  amplitude  spec¬ 
trum  has  a  single  peak,  which  is  somewhat  broader  than  that  associated 
with  the  secondary  power  density  spectra  in  Figs.  4.7  and  4.8.  The 
crosspower  amplitude  spectrum  for  0°  and  1.8°  has  a  peak  at  about 
60.9  kHz,  which  is  about  8%  lower  in  frequency  than  the  on-axis  value. 
The  crosspower  amplitude  spectrum  for  0°  and  3°  has  a  peak  at  about 
54.6  kHz,  which  is  lower  in  frequency  than  the  on-axis  value  by  about 
17%. 


The  crosspower  phase  spectrum  is  the  coherent  part  of  the 
crosspower  spectrum  and  provides  an  extremely  sensitive  comparison  of 
signals.  The  crosspower  phase  spectra  for  the  on-axis  and  off-axis 
signals  at  1.8°  and  3.0°  are  shown  in  Fig.  4.24.  The  phase  spectra  are 
very  similar,  both  relatively  constant,  and  nearly  zero  over  the  band. 
This  is  not  surprising  since  the  phase  spectra  were  derived  from  the 
crosscorrelation  function,  which  was  a  maximum  at  t=0.  This  condition 
implies  phase  coherence.  The  rapid  variations  in  phase  at  both  high 
and  low  frequencies  are  due  to  poor  S/N.  It  is  difficult  to  draw  any 
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FIGURE  4.22 

CROSSCORRELATION  OF  ON-AXIS  AND  OFF-AXIS  SIGNALS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

f(t)  -  GAUSSIAN,  f,c  =  600  kH*.  SPL,  =  195.2  dB  re  1  /uP a 
a  =  8.89  cm,  RQ  =  10.07  m,  r  =  23  6  m 
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(a)  ON-AXIS  CROSS  1.8°  OFF-AXIS 


FIGURE  4.24 

CROSS  POWER  PHASE  SPECTRA  OF  ON-AXIS  AND  OFF-AXIS  SIGNALS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

f(tl  -  GAUSSIAN.  f,c  =  600  kHz,  SPL,  =  195.2  dB  re  1  fiPa 
a  =  8.89  cm,  RQ  =  10.07  m,  r  -  23.5  m 
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further  conclusions  from  these  data  due  to  limitations  in  the  resolution 
of  the  analysis  equipment.  The  maximum  resolution  in  time  delay  of  the 
dual  channel  processor  is  3.91  psec,  which  was  extended  by  external  ana¬ 
log  delay  lines  to  0.5  psec.  This  time  delay  resolution  indicates  a 
maximum  phase  determination  accuracy  of  18.0°  at  100  kHz. 

C.  The  Broadband  Parametric  Array  with  High  Amplitude  Primary 
Waves 


As  the  amplitude  of  the  primary  waves  is  increased,  conditions 
are  reached  where  the  behavior  of  the  primary  signal  is  dominantly  non¬ 
linear.  Under  these  conditions,  the  primary  radiation  exhibits  finite 
amplitude  effects  and  the  broadband  parametric  transmission  is  subject 
to  finite  amplitude  attenuation.  The  primary  waveform  distorts  pro¬ 
gressively  with  propagation  range,  with  the  higher  amplitude  portions 
distorting  more  rapidly  than  lower  amplitude  portions,  resulting  in  wave 
steepening.  As  the  waveform  distorts,  energy  at  the  fundamental  is  con¬ 
verted  to  energy  at  harmonic  frequencies  which  suffer  higher  attenuation 
with  range  as  the  waves  propagate.  The  interaction  of  the  harmonic  com¬ 
ponents  generates  energy  in  the  secondary  band;  however,  the  conversion 
efficiency  is  poor  so  that  the  secondary  energy  that  results  from  inter¬ 
action  of  harmonic  frequencies  is  negligible. 

The  experimental  data  presented  here  are  high  amplitude 
extensions  of  conditions  that  were  evaluated  to  determine  low  amplitude 
performance  of  the  broadband  parametric  array. 

1.  The  Primary  Field 

The  Gaussian  pulse  primary  transmission  was  evaluated  at 
23.6  m  in  the  time  domain  as  a  function  of  amplitude.  The  on-axis 
acoustic  waveforms  are  shown  in  Fig.  4.25.  The  narrowband  presentation 
was  achieved  with  a  bandpass  filter  set  at  600  kHz  with  a  4-pole 
Butterworth  response.  The  broadband  presentation  was  achieved  with  a 
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FIGURE  4.25 

ON-AXIS  PRIMARY  WAVEFORMS  AS  A  FUNCTION  OF  SOURCE  LEVEL 

f(t)  -  GAUSSIAN.  flc  =  600  kHz,  a  =  8.89  cm.  RQ  =  10.07  m,  r  =  23.6  m 
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high  pass  filter  set  at  1  kHz  with  an  8-pole  Butterworth  response.  At 
low  level,  both  the  narrowband  and  the  broadband  waveforms  are  very 
similar.  As  the  level  is  increased,  waveform  distortion  is  observed. 

In  the  broadband  presentation,  the  secondary  signal  that  is  generated, 
the  harmonic  signals  that  are  converted  from  the  fundamental,  and  the 
remaining  primary  signal  are  all  contained  in  a  complex  waveform.  It 
should  be  noted  that  the  various  frequency  components  are  weighted 
according  to  the  response  curves  for  the  ARL  239-1  high  frequency  and 
low  frequency  receive  systems  (Appendix  B) .  In  the  narrowband  presenta¬ 
tion,  energy  is  converted  from  the  fundamental  to  harmonic  frequencies 
more  rapidly  at  the  higher  amplitude  peaks,  causing  a  rounding  off  of 
the  pulse  envelope.  This  rounding  off  or  change  in  pulse  envelope  shape 
effectively  increases  the  apparent  duration  if  the  pulse  transmission  is 
defined  by  the  0.17  amplitude  level  of  the  envelope. 

'*  The  on-axis  sound  field  was  evaluated  in  a  narrowband 

frequency  as  a  function  of  input  electrical  power  to  determine  the 
various  operating  regions  of  the  primary  transmission.  This  evaluation 
was  performed  as  a  function  of  primary  directivity  by  using  the  four 
piston  radii  that  were  originally  identified.  The  on-axis  peak  sound 
pressure  level  at  23.6  m  as  a  function  of  level  for  these  configurations 
is  shown  in  Fig.  4.26.  Linear  performance  lines,  which  are  drawn 
through  each  data  set,  indicate  relative  differences  in  level  cor¬ 
responding  to  differences  in  directivity  for  each  piston  configuration. 
As  the  input  electrical  power  is  increased,  the  amount  of  energy  that  is 
converted  from  the  fundamental  to  harmonic  frequencies  over  the  measure¬ 
ment  range  increases  due  to  progressive  finite  amplitude  effects.  This 
conversion  of  energy  at  the  expense  of  the  fundamental  component  indi¬ 
cates  the  onset  of  waveform  saturation  and  is  noted  in  Fig.  4.26  by  the 
departure  of  the  data  from  linear  behavior.  These  same  trends  were 
studied  by  Shooter,  Muir,  and  Blackstock^  for  the  radiation  from  a 
single  circular  piston  source. 
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PRIMARY  PEAK  SOUND  PRESSURE  LEVEL  -  dB  re  1  //Pa 


The 'departure  from  linear  behavior  of  the  fundamental 
component  becomes  progressively  greater  as  the  amplitude  is  increased; 
however,  a  rapid  departure  is  noted  in  the  data  at  input  electrical 
power  levels  greater  than  about  23  dB  re  1  W.  The  rapid  departure 
from  linearity  that  is  noted  in  the  data  was  identified  arbitrarily  as 
an  indicator  of  those  operating  regions  for  the  primary  source  which  are 
dominantly  linear  and  dominantly  nonlinear.  The.  low  amplitude  analysis 
discussed  in  previous  sections  involved  primary  source  amplitudes, 
which  were  dominantly  linear  according  to  this  criterion.  The  exact 
status  of  the  primary  source  for  any  data  presented  may  be  determined 
from  Fig.  4.26. 

The  data  indicate  that,  for  higher  directivity  source 
configurations,  the  departure  from  linear  behavior  occurs  at  higher 
peak  sound  pressure  levels,  at  slightly  lower  input  power  levels,  and 
increases  more  rapidly  with  increasing  input  power  levels.  These  data 
were  taken  in  the  farfield  of  each  source  configuration.  If  the  satura¬ 
tion  takes  place  predominantly  in  the  farfield  or  spherical  spreading 
region  of  the  sources,  then  the  peak  sound  pressure  level  at  which 
departure  from  linear  behavior  begins  would  be  higher  for  the  smaller 
pistons  or  lower  directivity  sources.  This  trend  is  based  on  the 

assumption  that  plane  waves  saturate  at  lower  amplitude  levels  than 
78 

spherical  waves,  and  is  not  noted  in  the  data.  The  data  support  the 

conclusion  that  finite  amplitude  effects  are  occurring  predominantly  in 

the  nearfield  region  of  the  sources,  i.e. ,  over  distances  ''ess  than  Rq. 

This  conclusion  may  be  evaluated  by  presenting  the  peak  sound  pressure 

level  measurements  as  peak  source  levels  at  Rq  foi  each  piston,  and  is 

given  in  Fig.  4.27.  If  the  simple  model  is  assumed  where  the  radiation 

in  the  nearfield  is  plane  wave  and  spherical  wave  in  the  farfield,  then 

the  data  in  Fig.  4.27  represent  the  level  at  the  end  of  the  plane  wave 

region,  ignoring  finite  amplitude  effects  that  may  have  occurred  in  the 

farfield.  With  this  presentation,  the  departure  from  linear  behavior 

occurs  at  higher  peak  source  levels  for  lower  directivity  sources. 

Since  the  lower  directivity  sources  have  shorter  values  for  R  and 

o 
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PRIMARY  PEAK  SOUND  PRESSURE  LEVEL  -  dB  re  1  tt Pa  at  R 


FIGURE  4.27 

PRIMARY  PEAK  SOURCE  LEVEL  AT  R0 

f(t>  -  GAUSSIAN,  f1c  =  600  kHz,  r  =  23.6  m 
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finite  amplitude  distortion  is  progressive,  this  trend  supports  the 
conclusion  that  the  finite  amplitude  effects  occur  predominantly  in 
the  nearfield  or  plane  wave  region  of  the  source.  The  effects  in 
the  farfield  or  spherical  region  are  noted  in  the  rate  at  which  the 
departure  from  linear  behavior  occurs  as  a  function  of  input  electri¬ 
cal  power.  In  the  farfield  region  for  the  lower  directivity  sources, 
spherical  spreading  slows  the  progressive  finite  amplitude  effects 
and  decreases  the  rate  of  departure  from  linear  behavior  relative  to 
the  higher  directivity  sources. 

Another  indication  of  finite  amplitude  effects  is  noted  in  the 
narrowband  energy  beam  pattern  at  the  fundamental  frequency.  On  the 
maximum  response  axis,  the  highest  amplitudes  are  reached  and  the  finite 
amplitude  effects  are  greatest.  The  conversion  of  energy  from  the 
fundamental  to  harmonic  components  causes  a  blunting  of  the  main  lobe 
and  apparent  increases  in  half-power  beamwidth  and  in  sidelobe  levels. 
These  phenomena  are  observed  in  the  primary  beam  patterns  as  a  function 
of  increasing  peak  sound  pressure  level.  The  source  directivity  of  the 
primary  transmission  as  a  function  of  source  level  from  the  selectable 
radius  piston  projector  with  radius  a  =  8.89  cm  is  shown  in  Fig.  4.28. 
The  energy  beam  pattern  in  Fig.  4.28(a)  was  taken  in  a  linear  operating 
region  according  to  Fig.  4.26  and  has  a  measured  half-power  beamwidth  of 
about  0.85°  with  maximum  first  sidelobe  levels  of  -13  dB.  The  energy 
beam  patterns  in  Figs.  4.28(b)  and  4.28(c)  were  taken  in  a  dominantly 
nonlinear  operating  region  and  indicate  a  broadening  of  the  half -power 
beamwidth  by  about  30%  and  increase  in  sidelobe  levels  of  about  2.5  dB. 
Source  directivity  data  for  the  other  configurations  of  the  selectable 
radius  circular  piston  projector  are  shown  in  Figs.  4.29,  4.30,  and 
4.31.  These  data  indicate  performance  similar  to  that  presented  in 
Fig.  4.28.  For  a  source  piston  radius  of  5.08  cm,  the  half-power  beam- 
width  broadened  by  about  20%  and  the  sidelobe  levels  increased  by  about 
3  dB.  An  approximate  11%  broadening  of  the  half -power  beamwidth  was 
noted  for  a  source  piston  radius  of  2.54  cm  and  the  sidelobe  levels 
increased  by  about  1  dB.  For  the  lowest  directivity  data  with  a  piston 
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FIGURE  4.28 

PRIMARY  SOURCE  DIRECTIVITY  AS 
A  FUNCTION  OF  SOURCE  LEVEL 
f(t)  -  GAUSSIAN,  f1c  -  600  kHz,  a  -  8.89  cm 
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(b)  SPL,  -  194.8  dB  re  1  fiPt 


FIGURE  4.29 

PRIMARY  SOURCE  DIRECTIVITY  AS 
A  FUNCTION  OF  SOURCE  LEVEL 
f(t)  -  GAUSSIAN,  f,c- 600  kHz,  a -5.08  cm 
R0“3.29  .1\  r  -  23.6  m 
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FIGURE  4.30 

PRIMARY  SOURCE  DIRECTIVITY  AS 
A  FUNCTION  OF  SOURCE  LEVEL 

f(t)  -  GAUSS'  f,c  -  600  kHz,  •  -  2.64  cm 

Rq  =  J.U2  m,  r  «  23.6  m 
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radius  of  1.27  cm,  the  broadening  of  the  half-power  beamwidth  was  about 
10%  and  the  sidelobe  levels  were  increased  by  about  0.5  dB.  A  summary 
of  primary  source  data  is  presented  in  Table  IV. 3. 

2 .  The  Secondary  Field 

The  secondary  field  characteristics  with  high  amplitude 
primary  waves  were  evaluated  for  several  primary  source  configurations 
and  directivity  ratios.  For  finite  amplitude  limited  parametric  trans¬ 
missions,  the  concept  of  directivity  ratio  is  somewhat  vague  since  the 
Rutherford  beamwidth  is  derived  from  low  amplitude  considerations.  The 
classification  of  the  radiation  from  a  parametric  source  as  plane  wave 
or  nonplane  wave  based  on  the  directivity  ratio  does  not  depend  on  the 
influence  of  finite  amplitude  effects.  The  use  of  the  directivity  ratio 
derived  from  low  amplitude  source  characteristics  may  be  used  to 
describe  finite  amplitude  limited  parametric  transmissions. 

The  broadband  parametric  array  with  high  amplitude  primary 
waves  and  high  primary  directivity  was  established  with  a  piston  radius 
of  8.89  cm  and  a  directivity  ratio  of  0.54.  The  secondary  waveforms 
both  on-axis  and  off-axis  are  shown  in  Fig.  4.32  for  primary  radiations 
with  a  peak  sound  pressure  level,  SPL^,  of  200.4  dB  re  1  uPa  at  23.6  m. 
This  value  of  primary  sound  pressure  level  is  clearly  dominated  by 
finite  amplitude  effects  as  about  5  dB  of  energy  at  the  fundamental  has 
been  converted  to  harmonic  frequencies  according  to  the  data  in 
Fig.  4.26.  A  comparison  of  the  on-axis  secondary  waveform  for  the  high 
and  low  (Fig.  4.4)  primary  amplitude  transmissions  indicates  several 
differences.  A  slight  increase  in  signal  duration  is  noted;  however, 
the  greatest  difference  is  in  the  reduction  of  amplitude  of  the  second 
positive  peak  in  the  waveform.  Both  of  these  features  result  from  the 
distortion  of  the  primary  envelope  due  to  finite  amplitude  effects.  As 
a  function  of  off-axis  angle,  the  relative  amplitude  of  the  second  posi¬ 
tive  peak  changes  little  over  angles  from  0.0°  to  3.25°.  This  behavior 
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FIGURE  4.32 

OFF-AXIS  ACOUSTIC  WAVEFORMS  FOR  THE  BROADBAND  PARAMETRIC 
ARRAY  WITH  HIGH  AMPLITUDE  PRIMARY  LEVEL 

GATED,  f(t)  -  GAUSSIAN,  f,c  =  600  kHz,  SPL,  =  200.4  dB  re  1  uPa 
a  =  8.89  cm,  RQ  =  10.07  m,  r  =  23.6  m 
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contrasts  with  the  decrease  and  final  disappearance  of  the  second 
positive  peak  in  the  data  for  low  amplitude  primary  levels. 

The  on-axis  and  off-axis  performance  noted  in  Fig.  4.32  was 
evaluated  in  the  frequency  domain  and  given  in  terms  of  the  normalized 
power  density  spectrum  at  each  angle  as  shown  in  Figs.  4.33  and  4.34. 
These  spectra  have  the  same  characteristic  Gaussian  shapes  as  noted  in 
the  low  amplitude  analyses.  The  frequency  and  level  at  the  spectral 
peak  were  extracted  from  these  spectra  and  are  shown  in  Fig.  4.35  with 
similar  data  for  other  primary  peak  sound  pressure  levels,  including  one 
that  is  considered  low  amplitude  as  well  as  the  theoretical  results  from 
Eq.  (3.45).  A  comparison  of  these  data  indicates  that  finite  amplitude 
effects  cause  the  on-axis  frequency  at  the  spectral  peak  to  decrease 
from  the  low  level  value.  For  the  data  shown,  the  decrease  is  from 
about  66  kHz  to  about  62  kHz,  or  about  6%.  As  the  off-axis  angle  is 
increased,  the  frequency  at  the  spectral  peak  under  finite  amplitude 
limitations  first  increases,  and  then  decreases  and  joins  the  low  level 
trend.  The  relative  level  at  the  spectral  peak  or  the  power  beam  pat¬ 
tern  indicates  a  decrease  in  directivity  as  the  primary  amplitude  is 
increased.  Measurements  of  the  energy  beam  pattern  as  a  function  of 
primary  peak  sound  pressure  level  are  shown  in  Fig.  4.36.  An  increase 
in  secondary  half -power  beamwidth  of  about  80%  is  noted  with  increasing 
primary  level.  The  decrease  in  frequency  at  the  spectral  peak  on-axis 
results  from  envelope  distortion  due  to  finite  amplitude  effects  in  the 
higher  amplitude  portions  of  the  Gaussian  pulse  transmission.  The 
decrease  in  secondary  directivity  is  due  to  the  finite  amplitude 
effects,  which  are  greatest  on  the  acoustic  axis.  Both  of  these  trends 
are  rooted  in  the  behavior  of  the  primary  field  and  will  continue  as  the 
amplitude  of  the  primary  waves  are  increased  toward  saturation. 

The  trends  established  in  the  analysis  of  high  amplitude 
primary  waves  with  high  primary  directivity  are  also  noted  in  the  data 
as  a  function  of  primary  directivity.  The  frequency  and  level  at  the 
spectral  peak  are  shown  in  Fig.  4.37  for  the  broadband  parametric  array 
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RELATIVE  LEVEL  AT  THE  FREQUENCY  AT  THE 
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NORMALIZED  OFF-AXIS  ANGLE,  0/03  dB  -  units 


FIGURE  4.37 

FREQUENCY  AND  LEVEL  AT  THE  SPECTRAL  PEAK  OFF-AXIS  FOR  THE 
BROADBAND  PARAMETRIC  ARRAY  WITH  HIGH  AMPLITUDE  PRIMARY  LEVEL 

f<t)  -  GAUSSIAN,  f,c  =  600  kHz,  a  -  5.08  cm,  R„  -  3.29  m,  r  =  23.6  m 
63dB  =0'80°  (MEASURED  LOW  AMPLITUDE) 


with  a  piston  radius  of  5.08  cm  and  a  directivity  ratio  of  0.89.  The 
frequency  at  the  spectral  peak  on-axis  is  again  observed  to  decrease 
from  66  kHz  for  low  amplitude  primary  waves  to  64  kHz  for  high  amplitude 
primary  waves.  The  relative  level  at  the  spectral  peak  indicates  a 
decrease  in  directivity.  This  trend  is  supported  in  the  measurements  of 
energy  beam  patterns  as  a  function  of  primary  peak  sound  pressure  level, 
which  are  shown  in  Fig.  4.38.  A  broadening  of  the  secondary  half-power 
beamwidth  of  about  66%  is  noted.  The  broadband  parametric  array  with  a 
directivity  ratio  of  1.74  was  formed  with  a  piston  radius  of  2.54  cm. 

The  frequency  and  level  at  the  spectral  peak  for  both  low  and  high 
amplitude  primary  transmissions  are  shown  in  Fig.  4.39.  These  data 
indicate  a  decrease  in  the  frequency  at  the  spectral  peak  on-axis  from 
about  68  kHz  to  66  kHz  or  about  3%,  and  a  decrease  in  directivity  at  the 
secondary  or  about  38%,  as  noted  in  the  energy  beam  patterns  shown  as  a 
function  of  primary  peak  sound  pressure  level  in  Fig.  4.40.  The  broad¬ 
band  parametric  array  with  the  lowest  directivity  and  a  directivity 
ratio  of  3.69  was  formed  with  a  piston  radius  of  1.27  cm.  The  frequency 
and  level  at  the  spectral  peak  are  presented  in  Fig.  4.41  where  the  fre¬ 
quency  at  the  spectral  peak  on-axis  does  not  decrease  under  the  influ¬ 
ence  of  the  high  amplitude  level  that  was  reached.  Only  a  slight 
broadening  of  6%  in  half -power  beamwidth  is  noted  in  the  energy  beam 
patterns  shown  in  Fig.  4.42;  however,  an  increase  in  the  level  far  off- 
axis  is  noted  in  Figs.  4.41  and  4.42. 

The  data  that  have  been  presented  cover  only  a  limited  range 
of  finite  amplitude  effects  and  directivity  ratios;  however,  certain 
trends  are  established  with  this  limited  set  of  data.  As  the  amplitude 
of  the  primary  waves  were  increased,  the  frequency  at  the  spectral  peak 
on-axis  generally  decreased  and  at  off-axis  angles,  the  downward  shift 
in  frequency  at  the  spectral  peak  was  reversed  or  inhibited.  These 
trends  result  from  waveform  induced  envelope  distortion  of  the  primary 
transmission.  The  secondary  directivity  generally  decreased  with  the 
onset  of  finite  amplitude  behavior  as  determined  by  an  increase  in  the 
half -power  beamwidth  and  level  far  off-axis.  The  magnitude  of  these 
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FIGURE  4.38 

SECONDARY  SOURCE  DIRECTIVITY  AS  A  FUNCTION 
OF  PRIMARY  PEAK  SOUND  PRESSURE  LEVEL 
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FIGURE  4.39 

FREQUENCY  AND  LEVEL  AT  THE  SPECTRAL  PEAK  OFF-AXIS  FOR  THE 
BROADBAND  PARAMETRIC  ARRAY  WITH  HIGH  AMPLITUDE  PRIMARY  LEVEL 

f(t)  -  GAUSSIAN,  f,  =  600  kHz,  a  =  2.54  cm,  RQ  =  0.82  m,  r  =  23.6  m 
63  dB  =  1  47°  (MEASURED  LOW  AMPLITUDE! 
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FIGURE  4.42 

SECONDARY  SOURCE  DIRECTIVITY  AS  A  FUNCTION 
OF  PRIMARY  PEAK  SOUND  PRESSURE  LEVEL 

fill  -  GAUSSIAN,  f,c  =  600  kHz.  a=  1.27  cm 
RQ  =  0,2 1  m.  r  =  23.6  m 
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trends  tends  to  diminish  as  the  primary  directivity  decreases  and  the 
directivity  ratio  increases.  This  is  supported  by  the  slower  rate  of 
departure  from  linear  behavior  for  the  primaries  with  decreasing  primary 
directivity. 

The  on-axis  performance  of  the  secondary  sound  field  with  high 
amplitude  primary  waves  is  considered  in  Fig.  4.43  where  the  secondary 
peak  sound  pressure  level  is  plotted  as  a  function  of  input  electrical 
power.  The  trends  in  this  figure  are  similar  to  those  in  Fig.  4.26, 
which  describes  the  on-axis  performance  of  the  primary  field.  The 
secondary  signal  level  is  proportional  to  the  square  of  the  low  ampli¬ 
tude  primary  signal  level,  which  accounts  for  the  steeper  slopes  in 
Fig.  4.43.  Departure  from  linear  behavior  is  a  maximum  for  the  broad¬ 
band  parametric  array  with  the  highest  directivity  and  varies  from  about 
10  dB  to  about  3  dB  for  the  data  shown. 

The  off-axis  performance  of  the  secondary  sound  field  with 
high  amplitude  primary  waves  is  considered  in  Fig.  4.44.  This  figure 
indicates  the  broadening  of  the  half-power  beamwidths  of  both  the  pri¬ 
mary  and  secondary  radiation  patterns  for  the  four  broadband  parametric 
arrays  that  were  formed.  For  each  array  formed,  the  secondary  beam 
experiences  greater  broadening  than  the  primary  beam.  This  is  another 
result  of  the  proportionality  of  secondary  level  to  the  square  of  the 
primary  level.  The  observed  half -power  beamwidths  at  low  level  are 

about  9%  broader  than  their  predicted  farfield  values  according  to  the 

25 

theory  of  Berktay  and  Leahy. 

A  cross  function  analysis  was  performed  with  the  secondary 
signals  resulting  from  the  transmission  of  high  amplitude  primary  waves. 
This  analysis  has  the  same  format  as  that  conducted  for  low  amplitude 
primary  wave  transmissions  in  Section  IV. B. 3.  This  analysis  allows  the 
determination  of  relative  signal  properties  using  coherent  signal 
analysis. 
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FIGURE  4.44 

PRIMARY  AND  SECONDARY  HALF  POWER  BEAMWIDtH 

f(t)  -  GAUSSIAN,  f,e  =  600  kHz,  f2e  3s 66  kHz,  r  =  23.6  m 
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The  selectable  radius  circular  piston  projector  was  driven 
using  the  full  aperture  with  a  radius  of  8.89  cm.  The  primary  peak 
sound  pressure  level,  SPL^,  was  200.4  dB  re  1  pPa  at  23.6  m.  The  cross¬ 
correlation  of  the  on-axis  and  off-axis  signals  at  1.8°(38»  )  and 

j  dB 

3°(59^(jg)  are  shown  in  Fig.  4.45.  As  the  angle  off-axis  is  increased, 
the  crosscorrelation  transitions  from  a  symmetric  signal  on-axis  to  a 
nonsymmetric  signal  off-axis.  The  first  positive  and  negative  peaks  are 
diminished  in  amplitude  relative  to  the  other  peaks  while  the  overall 
waveform  retains  the  duration  of  about  36  ysec.  This  behavior  is  simi¬ 
lar  to  that  observed  in  the  data  with  low  amplitude  primary  waves. 

The  signal  properties  were  evaluated  in  the  frequency  domain 
with  the  crosspower  density  spectrum.  The  crosspower  amplitude  spectra 
are  shown  in  Fig.  4.46  for  0.0°  and  1.8°,  and  for  0.0°  and  3.0°.  These 
amplitude  spectra  have  broad  peaks  that  result  from  the  product  of  two 
secondary  spectra.  The  spectral  peak  f rt  0.0°  and  1.8°  is  at  about 
60.4  kHz,  which  is  about  3%  lower  than  the  on-axis  value;  and  the  spec¬ 
tral  peak  for  0.0°  and  3.0°  is  at  about  57.1  kHz,  which  is  about  8% 
lower  than  the  on-axis  value.  The  phase  spectra,  which  indicate  the 
relative  phase  between  two  signals  as  a  function  of  frequency,  are  shown 
in  Fig.  4.47.  These  spectra  are  reasonably  constant  with  a  value  of 
phase  difference  near  zero  across  the  band.  This  performance  is  not 
surprising  since  the  data  were  derived  from  crosscorrelation  functions 
that  were  adjusted  to  be  maximum  at  t=0  thereby  implying  phase  coher¬ 
ence.  It  was  noted  in  the  analyses  of  low  level  performance  that  these 
observations  are  limited  in  their  contribution  to  conclusions  which  may 
be  reached  concerning  the  data.  The  limitation  results  from  equipment 
accuracy,  which  restricted  the  time  delay  resolution  between  signals 
and  established  the  accuracy  of  relative  phase  measurements  at  about 
18°  at  100  kHz. 
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6  jUsec/div 


(a)  ON-AXIS  CROSS  1.8°  OFF-AXIS 


6  pisec/div 


(b)  ON-AXIS  CROSS  3.0°  OFF-AXIS 


FIGURE  4.45 

CROSSCORRELATION  OF  ON-AXIS  AND  OFF-AXIS  SIGNALS 
FOR  THE  BROADBAND  PARAMETRIC  ARRAY 
WITH  HIGH  AMPLITUDE  PRIMARY  LEVEL 

f(t)  -  GAUSSIAN.  f,c  =  600  kHz,  SPL,  »  200  4  08  re  1  jzPa 
a  =  8.89  cm.  RQ  =  10  07  m,  r  =  23  6  m 
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CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 


The  parametric  acoustic  array  is  the  sound  field  of  the 
interaction  products  that  result  from  the  mixing  of  two  or  more  fre¬ 
quency  components  in  a  common  volume  of  an  acoustic  medium.  In  the 
context  of  this  investigation,  the  sound  field  of  the  baseband  inter¬ 
action  products  is  characterized.  An  inhomogeneous  wave  equation, 
which  is  derived  by  perturbation  of  the  exact  equations  of  acoustics, 
leads  to  analytical  descriptions  of  several  types  of  parametric  arrays. 
The  analytical  description  of  parametric  arrays  may  be  grouped  into 
four  basic  classifications  which  depend  on  both  spatial  and  temporal 
boundary  conditions.  These  groups  include  both  the  narrowband  and 
broadband  parametric  arrays  for  both  plane  wave  and  nonplane  wave 
fields.  Numerous  studies  have  been  conducted  on  the  narrowband  para¬ 
metric  array  for  both  plane  wave  and  nonplane  wave  fields;  however  only 
a  few  studies,  which  were  mostly  theoretical,  have  considered  the 
broadband  parametric  array.  The  characteristics  of  the  broadband 
parametric  array  for  nonplane  wave  fields  have  not  generally  been 
determined  and  are  the  subject  of  the  present  theoretical  and  experi¬ 
mental  investigation. 

The  broadband  parametric  array  for  plane  wave  fields  only 
approximates  the  physically  realizable  primary  source  configurations. 
The  analytical  treatment  for  the  broadband  parametric  array  for  plane 
wave  fields  was  extended  to  allow  detailed  analyses  of  the  sound  field 
in  both  the  time  and  frequency  domains.  The  analytical  descriptions 
were  cast  in  terms  of  a  linear  system  where  the  general  solution  off- 
axis  in  the  frequency  domain  is  given  by  the  product  of  the  on-axis 
solution  with  an  off-axis  transfer  function. 
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The  assumption  of  a  Gaussian  pulse  primary  transmission 
allows  detailed  evaluation  of  the  on-axis  and  off-axis  signal  character¬ 
istics  of  a  broadband  parametric  transmission.  A  closed  form  solution 
for  the  frequency  at  the  spectral  peak  on-axis  was  derived.  This  peak 
frequency  of  the  secondary  transmission  was  shown  to  be  directly 
proportional  to  the  primary  center  frequency  and  inversely  proportional 
to  Q  .  The  secondary  transmission  was  shown  to  have  a  constant  of 
1.2  regardless  of  the  primary  pulse  length  or  the  primary  or  secondary 
center  frequency.  It  was  determined  that  the  absolute  bandwidth  of  the 
secondary  transmission  is  approximately  the  same  as  that  of  the  primary 
transmission.  General  expressions  for  the  spectral  properties  off-axis 
were  derived  that  indicate  a  shift  toward  lower  frequencies  with 
increasing  off-axis  angle. 

The  broadband  parametric  array  for  nonplane  wave  fields 
describes  the  physically  realizable  parametric  array  that  results  from 
a  primary  transmission  from  a  real  source.  The  source  strength  density 
was  developed  for  a  pulsed  primary  transmission  that  is  spherically 
spreading  from  the  origin  and  which  has  an  off-axis  transfer  function 
that  is  valid  throughout  the  sound  field.  These  restrictions  are 
generally  met  in  the  farfield  of  a  planar  source.  A  general  solution 
for  the  secondary  pressure  was  given  in  a  form  that  involves  the  super¬ 
position  of  primary  frequency  field  components.  This  solution  involves 
a  four  dimensional  integration  over  space  and  frequency.  Using 
spherical  geometry,  the  general  solution  was  reduced  from  a  four¬ 
dimensional  integration  to  a  two-dimensional  integration  by  assuming 
symmetry  about  the  acoustic  axis  and  restricting  the  region  of  validity 
to  the  farfield  of  the  parametric  array. 

The  farfield  solution  was  expressed  for  a  Gaussian  pulse 
primary  transmission  from  a  circular  piston  source  and  limiting  forms 
of  this  solution  were  evaluated  as  a  function  of  directivity.  When  the 
primary  beams  are  extremely  narrow  and  act  as  spatial  impulse  functions, 
the  parametric  array  is  absorption  limited  and  the  broadband  parametric 

162 


array  for  plane  wave  fields  is  recovered.  When  the  Rutherford  beamwidth 
is  extremely  narrow  compared  to  the  primary  beamwidth  and  acts  as  a 
spatial  impulse  function,  the  parametric  array  is  diffraction  limited 
and  the  directivity  is  determined  by  primary  directivity  products  or 
off-axis  transfer  functions.  The  off-axis  spectral  properties  for  the 
diffraction  limited  case  were  evaluated  by  numerical  integration  of  the 
theoretical  off-axis  transfer  function  for  a  primary  signal  transmitted 
from  a  circular  piston  in  an  infinite  rigid  baffle  and  in  closed  form 
by  approximating  the  main  lobe  of  the  off-axis  transfer  function  by  a 
Gaussian  function.  As  the  characteristics  of  the  parametric  trans¬ 
mission  shift  from  the  absorption  limit  to  the  diffraction  limit,  the 
dependence  of  secondary  frequency  changes  from  quadratic  to  linear. 

With  this  change  in  secondary  frequency  dependence,  a  corresponding 
change  in  the  shift  toward  lower  frequencies  in  the  power  spectra 
occurs.  This  shift  is  characterized  by  a  reduction  in  the  on-axis 
frequency  at  the  spectral  peak  by  a  factor  of  /2  for  a  Gaussian  pulse 
primary  transmission.  In  the  absorption  limit,  the  characteristics  of 
the  secondary  off-axis  transfer  function  are  determined  by  the  absorp¬ 
tion  at  the  primary  frequencies  and  are  monotonically  decreasing  with 
angle  off-axis  and  secondary  frequency.  In  the  diffraction  limit,  the 
specific  primary  off-axis  transfer  function  is  dominant  in  determining 
the  characteristics  of  the  secondary  off-axis  transfer  function.  The 
characteristics  of  the  secondary  off-axis  transfer  function  generally 
follow  that  of  the  primary  off-axis  transfer  function  with  off-axis 
angle  but  are  extremely  insensitive  to  secondary  frequency. 

A  detailed  experimental  study  of  the  broadband  parametric 
array  provided  measured  data  with  variations  in  experimental  parameters 
and  emphasized  the  behavior  as  a  function  of  primary  directivity.  A 
selectable  radius  circular  piston  projector  was  used  as  the  primary 
source  at  600  kHz  to  transmit  Gaussian  pulses  from  14  to  17  usee  in 
duration.  Secondary  transmissions  with  a  near  1  and  an  on-axis 
peak  frequency  of  about  66  to  68  kHz  were  measured.  Data  were  taken 
for  directivity  ratios  (ratio  of  low  level  primary  half-power  beamwidth 
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to  the  Rutherford  beamwidth)  of  0.54,  0.89,  1.74,  and  3.69.  The 
transition  from  plane  wave  field  (absorption  limit)  to  nonplane  wave 
field  (diffraction  limit)  performance  mostly  occurs  for  directivity 
ratios  from  about  0.5  to  1.5. 

The  broadband  parametric  array  formed  with  a  directivity  ratio 
of  0.54  is  a  physical  approximation  to  the  performance  for  plane  wave 
fields.  This  parametric  array  was  evaluated  with  low  amplitude  primary 
waves  in  both  the  time  and  frequency  domains.  The  secondary  waveforms 
showed  the  characteristic  pulse  shape  of  the  second  time  derivative  of 
the  square  of  the  Gaussian  primary  envelope  on-axis  which  shifts  to 
the  first  time  derivative  off-axis.  This  is  the  trend  that  was  pre- 

9 

dieted  by  Westervelt.  In  the  frequency  domain,  power  density  spectra 
of  measured  signals  were  evaluated.  As  the  angle  off-axis  was 
increased  both  the  frequency  and  level  at  the  spectral  peak  decreased. 
The  frequency  at  the  spectral  peak  was  relatively  constant  at  small 
angles  with  a  steady  downward  shift  in  the  peak  frequency  at  angles 
greater  than  the  half-power  angle.  The  relative  level  at  the  spectral 
peak  was  smooth  and  monotonically  decreasing.  Measured  data  were 
compared  with  theoretical  results  for  the  farfield  absorption  and 
diffraction  limits  and  for  the  farfield  solution  from  Eq .  (3.45),  which 
was  evaluated  by  double  numerical  integration.  In  the  absorption  limit, 
the  predicted  peak  frequency  on-axis  is  higher  than  observed.  The  trend 
toward  a  downward  shift  in  frequency  at  the  spectral  peak  is  slow  for 
small  angles,  increases  to  a  rapid  shift,  and  becomes  relatively  con¬ 
stant  at  large  angles.  The  predicted  level  at  the  spectral  peak  in  the 
absorption  limit  is  slightly  lower  than  the  experimental  data.  In  the 
diffraction  limit,  the  predicted  frequency  at  the  spectral  peak  is  lower 
than  that  observed  on-axis  and  relatively  constant  out  to  the  half-power 
angle.  At  larger  angles,  secondary  maxima  in  the  frequency  at  the 
spectral  peak  are  noted  which  result  from  spatial  lobes  in  the  off-axis 
transfer  function.  Comparable  behavior  at  small  angles  is  noted  with  a 
Gaussian  approximation  of  the  main  lobe  of  the  off-axis  transfer 
function.  The  level  at  the  spectral  peak  decreases  very  rapidly 
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off-axis  for  both  diffraction  limit  representations.  The  theoretical 
results  obtained  with  the  farfield  solution  from  Eq.  (3.45)  indicate 
generally  good  agreement  with  the  experimental  data.  Excellent  agree¬ 
ment  was  noted  in  the  frequency  at  the  spectral  peak  at  small  angles. 

At  larger  angles  the  theoretical  results  are  in  general  agreement  but 
predict  slightly  lower  frequencies  at  the  spectral  peak.  This  differ¬ 
ence  is  attributed  to  the  higher  measured  sidelobe  levels  in  the  primary 
off-axis  transfer  function  than  those  assumed  in  the  theoretical  repre¬ 
sentation  of  a  piston  in  an  infinite  rigid  baffle.  The  higher  sidelobe 
levels  tend  to  retard  the  downward  shift  in  frequency  at  the  spectral 
peak  in  the  measured  data.  The  relative  level  at  the  spectral  peak  is 
in  excellent  agreement  with  the  theoretical  results  predicted  by 
Eq.  (3.45).  At  the  half-power  angle,  the  agreement  with  the  absorption 
limit  is  also  good.  This  agreement  was  noted  in  earlier  studies  as 
verification  of  the  Westervelt  model  for  the  parametric  array  for  plane 
wave  fields.  The  agreement  with  the  absorption  limit  degrades  at  large 
angles  and  supports  the  conclusion  that  this  parametric  array  is  only  a 
physical  approximation  for  plane  wave  fields. 

The  general  characteristics  that  have  been  described  in  the 

comparison  of  experimental  and  theoretical  results  are  also  noted  in  the 

data  for  the  other  source  configurations  which  transition  toward  the 

diffraction  limit.  These  other  source  configurations  are  characterized 

1  68 

by  a  shift  in  overall  weighting  of  the  secondary  frequency  from  u  '  to 

1.3  ^ 

0)2"  .  Recall  that  the  dependence  is  quadratic  in  the  absorption  limit 

and  linear  in  the  diffraction  limit.  In  each  data  set,  the  influence 
of  the  primary  off-axis  transfer  function  increased  as  the  diffraction 
limit  was  approached  and  a  trend  toward  frequency  independence  as  a 
function  of  angle  off-axis  was  observed.  The  experimental  data  for  a 
directivity  ratio  of  3.69  represent  the  most  diffraction  dominated  case 
that  was  evaluated.  These  experimental  data  indicated  a  secondary  maxi¬ 
mum  in  the  frequency  at  the  spectral  peak,  which  is  in  good  agreement 
with  the  predictions  of  Eq.  (3.45)  and  also  agrees  in  angular  placement 
with  the  diffraction  limit  of  Eq.  (3.45).  No  perceptible  sidelobes  were 
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predicted  from  Eq.  (3.45)  in  the  level  at  the  spectral  peak;  however, 
there  was  such  a  sidelobe  in  the  measured  data.  This  again  results 
from  the  higher  sidelobe  levels  in  the  measured  primary  off-axis  trans¬ 
fer  or  directivity  functions.  The  comparisons  between  the  experimental 
data  and  the  theoretical  results  indicate  that  the  frequency  at  the 
spectral  peak  is  a  more  sensitive  indicator  of  the  characteristics  of 
the  broadband  parametric  transmission  than  the  level  of  the  spectral 
peak  which  as  a  function  of  off-axis  angle  indicates  the  power  beam 
pattern. 

The  broadband  parametric  array  was  evaluated  with  high 
primary  levels  as  a  function  of  primary  directivity.  An  evaluation  of 
the  primary  field  indicated  the  influence  of  finite  amplitude  effects 
as  a  function  of  increased  levels  with  a  departure  from  linear  per¬ 
formance  of  the  narrowband  primary  transmission.  An  analysis  of  the 
on-axis  field  indicated  that  most  of  the  finite  amplitude  distortion 
occurred  in  the  nearfield  of  the  primary  source.  Finite  amplitude 
effects  compete  with  spherical  spreading  and  small  signal  absorption; 
therefore,  as  the  primary  directivity  decreased,  so  did  the  influence 
of  finite  amplitude  effects  on  the  primary  waves.  In  the  primary  field, 
finite  amplitude  effects  cause  envelope  distortion  of  a  narrowband 
primary  signal,  a  broadening  of  the  main  primary  beam,  and  an  apparent 
increase  in  sidelobe  levels. 

The  secondary  field  was  evaluated  in  the  time  and  frequency 
domains,  where  many  of  the  same  trends  were  indicated  as  observed  with 
low  amplitude  primary  waves.  The  frequency  at  the  spectral  peak,  which 
was  extracted  from  measured  power  spectra,  was  lower  on-axis  than  the 
value  with  low  amplitude  primary  waves  d.'e  to  the  envelope  distortion 
on-axis.  As  the  off-axis  angle  increased,  the  frequency  at  the 
spectral  peak  first  increased,  then  decreased  and  followed  the  low 
amplitude  trend.  A  broadening  of  the  secondary  beams  was  also  noted. 


A  comparison  of  the  experimental  data  for  each  source 
configuration  with  the  theoretical  results  of  Eq.  (3.45)  indicates 
several  deficiencies  in  the  analytical  modeling.  For  low  amplitude 
primary  waves,  the  on-axis  experimental  and  theoretical  results  agree 
very  well;  however,  the  off-axis  theoretical  results  agree  very  well 
for  the  level  at  the  spectral  peak  but  generally  predict  lower  values 
of  frequency  at  the  spectral  peak  than  observed.  This  difference 
results  from  the  assumption  of  theoretical  off-axis  transfer  functions 
for  a  piston  in  an  infinite  rigid  baffle  rather  than  an  exact  represen¬ 
tation  of  the  measured  off-axis  transfer  function  at  the  primary  fre¬ 
quencies.  There  were  also  slight  differences  in  the  measured  and 
predicted  values  of  SPL^  that  could  not  be  attributed  to  experimental 
error.  These  level  differences  result  primarily  from  the  incompleteness 
of  the  parametric  array  at  the  measurement  range.  To  account  for  these 
differences  the  farfleld  approximation  must  not  be  assumed  and  the 
range  integration  performed. 

The  correct  modeling  of  the  broadband  parametric  array  with 
high  amplitude  primary  waves  is  a  more  complex  problem  than  the  low 
amplitude  primary  wave  theoretical  treatment.  The  characteristics  of 
the  broadband  parametric  array  under  these  conditions  are  determined  by 
the  characteristics  of  the  primary  transmission  in  both  the  nearfield 
and  the  farfield  of  the  primary  source.  Finite  amplitude  effects  are 
cumulative  with  waveform  propagation;  therefore  the  characteristics  of 
the  primary  field  must  be  described  as  a  function  of  range  throughout 
the  nearfield  and  the  farfield  of  the  source.  This  characterization 
must  be  included  in  the  theoretical  modeling  and  the  range  integration 
performed . 
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APPENDIX  A 

RADIATION  OF  A  PULSED  CARRIER  IN  THE 
FARFIELD  OF  A  CIRCULAR  PISTON 


I 


The  radiation  in  the  farfield  of  a  circular  piston  set  in  an 
infinite  rigid  baffle  has  been  the  subject  of  many  theoretical  and 
experimental  studies.  The  symmetry  which  is  exhibited  in  the  radiation 
field  allows  straightforward  analyses  of  off-axis  characteristics.  In 
this  appendix,  the  radiation  in  the  farfield  of  a  circular  piston  will 
be  evaluated  to  determine  the  spectral  distortion  of  a  broadband  signal 
as  a  function  of  off-axis  angle.  This  study  will  allow  a  comparison  of 
the  off-axis  spectral  distortion  for  transmissions  from  a  conventional 
acoustic  source  and  a  parametric  array. 

A.  General  Solution 


The  time  independent  solution  for  the  pressure  field  produced 

by  radiation  in  the  farfield  of  a  circular  piston  set  in  an  infinite 

19 

rigid  baffle  is  given  by  Morse  and  Ingard.  Ignoring  losses  due  to 
dissipation,  the  solution  may  be  written  in  the  frequency  domain  in  the 
form 


R  -  j  k.  r 

P^r.e.f)  =  -f-  pQ  e  D(0,f) 


(A.  1) 


In  this  equation,  Pq  is  the  peak  pressure  amplitude  measured  at  Rq,  and 
D(0, f )  is  the  directivity  function  given  by 


D(6,f) 


2J^(k^a  sin6) 
k^a  sin0 


(A. 2) 


where  J^O)  is  the  Bessel  function  of  the  first  kind  and  first  order. 

An  analysis  of  the  spectral  distortion  of  a  broadband  signal 
which  is  transmitted  from  a  circular  piston  as  a  function  of  off-axis 
angle  is  facilitated  by  treating  the  on-axis  and  off-axis  properties  of 
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Eq.  (A.l)  as  the  input  and  output,  respectively,  of  a  linear  system 
with  the  form 


PjOr.e.f)  =  P1(r,0,f)  H10(f)  ,  (A. 3) 


where  Hin(f)  is  the  off-axis  transfer  function  for  the  first  order 
lo 

field.  By  comparison  with  Eq.  (A.l),  it  is  clear  that  the  off-axis 
transfer  function  in  the  farfield  of  a  circular  piston  set  in  an 
infinite  rigid  baffle  is  simply  the  directivity  function  given  by 


H10(f)  =  D(9’f) 


2Ji(kia  sin0) 
k^a  sin6 


(A. 4) 


2 

The  off-axis  power  transfer  function,  |H10(f)|  ,  is  shown  in  Fig.  A.l. 
Continuing  this  development,  the  solution  for  the  pressure  field  in  the 
time  domain  is  given  by 


p^r.e.t)  =  p^r.O.t)  *  h10(t)  ,  (A. 5) 


where  h 
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(t)  is  the  off-axis  impulse  response  function  given  by 


hieM  - 


(A. 6) 


where 


circ  |  X  | 


1  |X|  <  1 
0  f  X I  >  1 


(A. 7) 


Analysis  of  Eqs.  (A. 6)  and  (A. 7)  indicates  that  the  off-axis  impulse 
response  function  is  a  pulse  with  amplitude  and  duration  which  are 


functions  of  0.  On  axis,  the  impulse  response  function  becomes 


hi0(t)  =  <S(t) 


and  at  right  angles  to  the  acoustic  axis, 


hl*(t) 


(A. 8a) 


(A. 8b) 


The  duration  of  the  pulse  at  6-tt  is  t=2a/cQ  or  the  propagation  time  in 
the  medium  across  the  diameter  of  the  piston. 

The  separation  of  the  on-axis  and  off-axis  radiation 
characteristics  allows  the  analysis  of  signal  distortion  as  a  function 
of  off-axis  angle  in  both  the  time  and  frequency  domains.  This  analysis 
is  based  on  specified  or  measured  signal  characteristics  on  the  acoustic 
axis. 

B.  On-Axis  Radiation  of  a  Pulsed  Carrier  in  the  Farfield  of  a  Circular 
Piston 


The  on-axis  pressure  for  a  pulsed  carrier  signal  emitted  from 
a  circular  piston  set  in  an  infinite  baffle  may  be  written 


p^r.O.t) 


(r-R,) 


f(t  -  ■;-)  sin(„,lct-klcr) 

'  o' 


(A-9) 


where  is  the  absorption  coefficient  at  the  carrier  frequency  and 

a),  is  the  angular  carrier  frequency.  It  is  assumed  that  the  pulse 
lc 

envelope  function  f(t)  is  slowly  varying  with  respect  to  sin(w^ct) . 
Signal  level  dissipation  as  a  function  of  range  is  considered  with  the 
inclusion  of  absorption  at  the  carrier  frequency.  It  is  assumed  that 
this  signal  is  bandlimited  such  that  the  absorption  coefficient  at  each 


frequency  throughout  the  band  may  be  approximated  by  the  absorption 
coefficient  at  the  carrier  frequency.  In  the  frequency  domain,  the 
on-axis  pressure  signal  for  a  pulsed  carrier  is  given  by  the  Fourier 
transform  of  Eq.  (A. 9),  which  is  written 


Pj/r.O.f) 


(r-R,) 


(F 

-  t) 

*  JF 

sin(o>lct-klcr) 

\  O' 

.  (A. 10) 


This  equation  reduces  to 


R  -a.  (r-R  )  -jk1r 

Pl<r>°.t)  "if?  Po  e  1C  °  • 


r 

X 

SF 

f(o 

f=f1+f1  ■  ^ 

1  lc 

jf(t) 

f=frfic. 

(A. 11) 


Equation  (A. 11)  indicates  that  the  on-axis  frequency  spectrum  is  deter¬ 
mined  by  the  two  sided  Fourier  transform  of  the  pulse  envelope  function 
centered  at  the  carrier  frequency.  Further  study  of  the  spectral  dis¬ 
tortion  as  a  function  of  off-axis  angle  of  a  pulsed  carrier  signal 
emitted  from  a  circular  piston  set  in  an  infinite  rigid  baffle  requires 
the  specification  of  a  pulse  envelope  function.  The  assumption  of  a 
Gaussian  pulse  envelope  is  mathematically  tractable  and  a  reasonable 
approximation  to  pulse  envelopes  which  may  be  attained  in  practice. 


A  Gaussian  envelope  function  is  assumed  to  have  the  form 


f(t) 


2  2 
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(A. 12) 


which  has  a  Fourier  transform  given  by 


2 

=  F(f)  =  ^  exp^J)  f2j  ,  (A. 13) 

where  n  is  a  constant  that  determines  the  effective  duration  of  the 

pulse  and  therefore  the  effective  width  of  the  frequency  spectrum.  The 

spectral  characteristics  of  the  Gaussian  pulse  transmission  may  be 

specified  by  defining  the  constant  n  in  the  frequency  domain.  A  repre- 
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sentation  of  the  power  spectral  density  function  for  finite  duration 
signals  may  be  derived  from  the  energy  spectral  density  function  by 
weighting  the  energy  spectral  density  function  with  respect  to  the  dura¬ 
tion  of  the  signal,  or 
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(A. 14) 


where  P^^Xf)  is  the  power  spectral  density  function,  Exx(f)  is  the 
energy  spectral  density  function,  and  T  is  the  duration  of  the  signal. 
It  is  required  that  the  normalized  power  spectral  density  function  of 
the  on-axis  pulsed  carrier  transmission  be  equal  to  0.5  at  the  half¬ 
power  frequencies,  or 


PlxxN(f) 
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(A. 15) 


f=Fn  or  F12 


where  is  the  normalized  power  spectral  density  function  and 

and  are  the  lower  and  upper  half-power  frequencies,  respec¬ 
tively.  Substituting  Eqs.  (A. 11)  and  (A. 13)  into  Eq.  (A. 15)  and  solving 
for  n  gives 


n  =  2.668  |FU-F12| 


(A. 16) 
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where  (F^-F^l  t'^e  ^a^^"Power  bandwidth  of  the  pulsed  carrier 
transmission.  If  the  quality  factor  is  defined  by 


n  =  1c 

^1  Tf  -f 
1  11  12 


(A. 17) 


then  n  may  also  be  expressed  by 


n  =  2.668 


lc 


(A. 18) 


If  the  pulse  duration  is  defined  by  the  inverse  of  the 
half-power  bandwith,  then 


T 


(A. 19) 


and  the  pulse  duration  is  determined  between  the  0.17  amplitude  levels 
of  the  pulse  envelope  as  computed  from  Eq.  (A. 12)  for  t=T/2. 


Assuming  the  Gaussian  pulse  envelope  function  described  by 
Eq.  (A. 12),  the  on-axis  pressure  for  a  pulsed  carrier  signal  emitted 
from  a  circular  piston  set  in  an  infinite  rigid  baffle  is  found  from 
Eqs.  (A. 9)  and  (A. 12)  to  be  given  by 
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(A. 20) 


in  the  time  domain.  Tn  the  frequency  domain,  the  on-axis  solution  is 
found  from  Eqs.  (A. 11)  and  (A. 13)  and  is  given  by 
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C.  Off-Axis  Radiation  of  a  Pulsed  Carrier  in  the  Farfield  of  a 
Circular  Piston 


The  solution  for  the  off-axis  pressure  in  the  farfield  of  a 
circular  piston  set  in  an  infinite  rigid  baffle  is  given  in  the  time 
domain  by  the  convolution  of  the  on-axis  solution  with  the  off-axis 
impulse  response  function  and  in  the  frequency  domain  by  the  product  of 
the  on-axis  solution  with  the  off-axis  transfer  function. 


For  the  Gaussian  pulse  transmission,  the  off-axis  solution  in 
the  time  domain  is  found  by  substituting  Eqs.  (A. 6)  and  (A. 20)  into 
Eq.  (A. 5),  and  can  be  written  as 
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(A. 22) 


In  the  frequency  domain,  the  off-axis  solution  is  found  by  substituting 
Eqs.  (A. 4)  and  (A. 21)  into  Eq.  (A. 3)  and  can  be  written  as 
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Further  analysis  of  the  off-axis  spectral  distortion  is 
considered  for  the  set  of  parameter  values  given  in  Table  A.l.  These 
parameter  values  are  comparable  with  those  which  were  attained  in  the 
experimental  study  of  the  broadband  parameteric  array  discussed  in 
Section  IV.  The  results  that  are  presented  in  this  appendix  may  be 
compared  with  both  experimental  and  theoretical  results  given  in  Sec¬ 
tion  IV  to  determine  the  differences  between  linear  and  parametric 
implementations  of  broadband  sources  with  similar  characteristics. 

Figure  A. 2  shows  the  normalized  power  spectral  density  function  plotted 
as  a  function  of  frequency  for  various  off-axis  angles  and  a  Gaussian 
pulse  with  a  of  10.91.  This  value  of  is  often  attained  for  trans¬ 
ducers  used  in  underwater  applications.  On  axis,  the  spectrum  is  sym¬ 
metric  about  the  carrier  or  center  frequency  with  a  half-power  bandwidth 
of  6.05  kHz.  As  the  off-axis  angle  is  increased,  the  relative  spectral 
power  level  decreases  and  the  spectral  peak  shifts  downward  in  frequency 
for  off-axis  angles  up  to  1.50°.  It  should  be  noted  that  data  at  1.50° 
are  not  shown  in  the  figure.  This  results  because  the  level  at  the 
spectral  peak  at  this  angle  is  less  than  -20  dB.  Between  1.50°  and 
2.00°  the  relative  spectral  power  level  increases  and  the  spectral  peak 
shifts  upward  in  frequency  to  a  value  higher  than  the  carrier  frequency. 
As  the  off-axis  angle  increases  from  2.00°  to  2.75°,  the  spectral  peak 
again  shifts  downward  in  frequency  with  a  decrease  in  the  relative 
power  level.  These  characteristics  result  from  the  interaction  of  the 
relatively  narrowband  signal  with  the  lobe  structure  of  the  off-axis 
power  transfer  function  shown  in  Fig.  A.l. 

The  frequency  at  the  spectral  peak  is  determined  from  the 
normalized  power  spectral  density  function,  Eq.  (A. 25),  and  shown  as  a 
function  of  normalized  off-axis  angle  in  Fig.  A. 3.  The  angular  axis  is 
in  units  of  the  half-power  angle.  The  data  in  Fig.  A. 3  correspond  with 
that  in  Fig.  A. 2  and  initially  indicate  a  downward  shift  in  the  fre¬ 
quency  at  the  spectral  peak  for  increasing  off-axis  angle.  Included  in 
this  figure  is  the  relative  power  level  at  each  spectral  peak,  which  is 
also  computed  from  Eq.  (A. 25)  and  indicates  the  energy  or  power  beam 
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TABLE  A. 1 

EXAMPLE  PARAMETER  VALUES  FOR  RADIATION  FROM  A 
CIRCULAR  PISTON 


Carrier  or  center  frequency, 

Speed  of  sound,  c 

o 

Piston  radius,  a 

Half-power  beamwidth  at  66  kHz,  0  _ 

HP 

Primary  transmission  quality  factor, 
Half-power  angle  at  66  kHz, 


66.0  kHz 
1480  m/sec 
53  cm 
1.24° 

10.91  or  1. 
0.62° 
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NORMALIZED  OFF-AXIS  ANGLE,  0/03dB  -  units 


FIGURE  A.3 

FREQUENCY  AND  LEVEL  AT  THE  SPECTRAL  PEAK 
PREDICTED  OFF-AXIS  FOR  A  CIRCULAR  PISTON 

fit)  -  GAUSSIAN,  f1e  =  66  kHz,  Q,  -  10.91 
a -53  cm,  FARFIELD,  03dB=O.62° 
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pattern  for  the  transmission.  When  the  relative  power  level  at  the 
spectral  peak  has  decreased  3  dB,  the  frequency  at  the  spectral  peak 
has  shifted  downward  about  2.0%  to  a  value  of  approximately  64.7  kHz. 
This  occurs  at  an  off-axis  angle  of  about  0.62°  (63 jg) »  which  is  the 
half-power  angle  at  the  carrier.  The  frequency  at  the  spectral  peak 
shifts  downward  to  about  62.0  kHz,  or  about  6.0%  at  an  angle  of  1.49° 

(2.4  ®3(jg)»  where  the  frequency  shifts  discontinuously  upward  to  about 
69.0  kHz,  or  about  6.0%  with  respect  to  the  carrier  frequency.  The  dis¬ 
continuous  shift  upward  in  frequency  at  the  spectral  peak  results  from 
passing  through  a  null  in  the  off-axis  power  transfer  function.  The 
Gaussian  pulse  transmission  with  a  of  10.91  is  sufficiently  narrow- 
band  so  that  nulls  in  the  power  transfer  function  dominate  the  behavior 
of  the  off-axis  power  density  spectra.  The  frequency  of  the  spectral 
peak  again  decreases  from  1.49°  and  attains  the  value  of  the  carrier  fre¬ 
quency  at  an  angle  corresponding  to  a  peak  of  a  lobe  in  the  off-axis 
transfer  function.  The  frequency  of  the  spectral  peak  continues  to 
decrease  until  the  next  null  in  the  off-axis  transfer  function  becomes 
dominant.  The  trend  that  has  been  described  then  repeats  for  additional 
nulls. 

The  normalized  power  spectral  density  function  as  a  functior 

of  frequency  and  off-axis  angle  for  a  Gaussian  pulse  transmission  with 

a  of  1.2  is  shown  in  Fig.  A. 4.  This  value  of  is  difficult  to 

attain  in  practice  but  serves  to  demonstrate  the  characteristics  of  a 

broadband  linear  transmission  that  are  comparable  with  those  of  the 

broadband  parametric  transmission  w:  th  ~  constant  Qg  of  1.2.  On  axis, 

the  power  density  spectrum  is  symme ^  bout  the  carrier  frequency  with 

a  half-power  bandwidth  of  55.0  kHz.  As  -a  the  example  for  a  of 

10.91,  the  relative  power  levels  decrease  and  the  spectral  peak  shifts 

downward  in  frequency  as  the  off-axis  angle  increases.  This  trend 

continues  and  does  not  reverse  as  noted  for  the  higher  case.  The 

influence  of  lobes  in  the  off-axis  transfer  function  is  clearly  evident 

at  off-axis  angles  of  1.50°  (2.4  0_  ,_)  and  greater;  however,  the  Q, 

J  do  i 

of  the  Gaussian  pulse  transmission  is  sufficiently  low  that  the  relative 
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spectral  power  is  greater  at  lower  frequencies.  By  comparison  with  the 
high  transmission,  the  spectral  peak  shifts  downward  more  rapidly 
for  the  low  transmission,  but  the  decrease  in  relative  power  levels 
is  less  rapid  as  a  function  of  off-axis  angle.  These  trends  are  more 
explicit  in  Fig.  A. 5  where  the  frequency  at  the  spectral  peak  from 
Eq.  (A. 27)  and  the  relative  power  levels  at  each  frequency  from 
Eq.  (A. 25)  are  plotted  as  a  function  of  off-axis  angle.  For  the  lower 
case,  the  frequency  at  the  spectral  peak  and  the  relative  spectral 
power  level  both  decrease  with  increasing  off-axis  angle  to  about  3° 

(4.8  S^g)  where  the  frequency  at  the  spectral  peak  shifts  discontin- 
uously.  This  is  the  same  trend  as  noted  for  the  higher  transmission 
in  Fig.  A. 3;  however,  the  frequency  shifts  are  greater  and  the  level  at 
the  spectral  peak  is  monotonically  decreasing  with  increasing  off-axis 
angle.  This  trend  indicates  the  dominance  of  the  properties  of  the 
relatively  low  transmission  over  those  of  the  off-axis  power  trans¬ 
fer  function.  When  the  relative  power  level  at  each  spectral  peak  off- 
axis  has  decreased  3  dB,  the  frequency  at  the  spectral  peak  has  shifted 
downward  about  18.2%  from  the  on-axis  value  (66  kHz)  to  a  value  of 
about  54.0  kHz.  This  occurs  at  an  off-axis  angle  of  about  0.68°,  which 
is  about  10%  greater  than  the  half-power  angle  of  the  carrier.  At  a 
relative  power  level  of  -10  dB,  the  frequency  at  the  spectral  peak  has 
shifted  downward  about  56.0%  to  a  value  of  29.0  kHz.  This  occurs  at  an 
angle  of  1.7°,  which  is  greater  than  the  half-power  angle  at  the  carrier 
by  a  factor  of  about  2.7. 
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APPENDIX  B 

EXPERIMENTAL  EQUIPMENT 
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Details  regarding  the  design,  construction,  and  performance  of 
special  purpose  experimental  equipment  are  discussed  in  this  appendix. 
Included  in  this  discussion  are  individual  electronic  devices,  electro¬ 
mechanical  transducers,  and  subsystem  performance. 

A.  Transmit  Section 


A  block  diagram  of  the  transmit  section  is  shown  in  Fig.  B.l. 
The  signal  source  at  600  kHz  is  an  HP3330B  frequency  synthesizer  which 
is  input  to  the  shaped  pulse  generator.  A  detailed  description  of  the 
basic  components  of  the  transmit  section  follows. 

1.  Shaped  Pulse  Generator.  The  shaped  pulse  generator  allows 
the  production  of  pulsed  carrier  transmissions  with  selectable  envelope 
shapes  and  durations. 

A  block  diagram  of  the  shaped  pulse  generator  is  shown  in 
Fig.  B.2.  The  central  element  is  a  programmable  read  only  memory  (PROM) 
which  is  used  to  store  symmetric  but  otherwise  arbitrary  bandlimited 
envelope  shapes.  One  half  of  the  envelope  is  stored  in  the  PROM.  The 
complete  envelope  is  formed  by  sequentially  addressing  all  locations 
in  the  PROM  first  in  the  forward  direction  and  then  in  reverse.  The 
locations  are  addressed  by  a  set  of  up/down  counters  that  are 
incremented  by  a  clock,  which  is  part  of  a  timing  and  synchronization 
circuit.  The  frequency  of  the  shape  clock  determines  the  duration  of 
the  envelope,  which  varies  from  5  psec  to  10  msec.  Another  portion  of 
the  timing  and  synchronization  circuit  controls  the  pulse  repetition 
rate,  which  may  be  adjusted  from  10  pulses/sec  to  5000  pulses/sec,  or 
may  be  externally  adjusted.  The  pulse  repetition  rate  trigger  may  be 
synchronized  with  either  the  initiation  of  an  envelope  sequence  or  with 
the  carrier  frequency. 
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The  output  of  the  PROM  is  input  to  a  D/A  converter,  which 
provides  a  baseband  analog  envelope.  This  signal  is  coupled  to  a 
1.0  MHz  low  pass  filter  to  reduce  the  switching  noise  that  is  common  in 
high  speed  digital  circuitry.  The  filtered  baseband  envelope  signal  is 
multiplied  by  a  carrier  wave  to  produce  a  shaped  carrier  pulse.  This 
pulse  is  input  to  a  1.2  MHz  low  pass  filter  for  further  noise  reduction. 

2.  Power  Amplifier.  The  output  of  the  shaped  pulse  generator  is 
passed  through  a  bandpass  filter  for  further  noise  reduction  and 
coupled  to  an  LPPV5-4X,  2000  V*A  power  amplifier.  This  power  amplifier 
is  a  solid  state  design  with  eight  250  V*A  channels.  A  variable  impe¬ 
dance  output  transformer  is  achieved  through  various  series  and  parallel 
combinations  of  individual  output  transformers  from  each  channel.  The 
resultant  impedance  range  is  from  1.5  to  192  Q  in  eight  steps.  The  out¬ 
put  of  the  power  amplifier  is  coupled  to  the  projector  through  an  addi¬ 
tional  impedance  matching  transformer  and  tuning  elements  to  provide 
maximum  power  transfer  for  each  configuration  of  the  projector  array. 

3.  Source  Transducer.  The  source  transducer  (ARL  384-1)  was 
designed  to  allow  the  formation  of  parametric  arrays  for  both  plane 
wave  (approximately)  and  nonplane  wave  fields. 

a.  Transducer  Design.  The  active  transducer  elements  are 
12  concentric  circular  rings  and  a  center  piston  which  are  driven  in  the 
thickness  mode  and  may  be  selectively  driven  electrically  to  form  pis¬ 
tons  of  different  radii.  The  ceramic  rings  have  a  thickness  of  0.32  cm 
and  are  equally  spaced  with  a  radial  width  of  0.56  cm.  The  rings  and 
the  center  piston  were  constructed  from  four  polarized  squares  of  Chan- 
nelite  5400  ceramic,  8.89  cm  on  a  side.  The  individual  squares  were 
separated  by  0.08  cm,  arranged  into  a  larger  composite  square,  and 
mounted  on  a  flat  surface.  Using  circular  cutting  tools  with  different 
radii,  the  concentric  rings  and  the  center  piston  were  cut  by  grinding 
away  a  0.08  cm  thickness  using  carborundum  grit  mixed  with  water.  This 
technique  allowed  the  ceramic  to  be  cut  without  depolarization.  The 


ceramic  pieces  were  placed  face  down  on  a  flat  surface,  separated  by 
0.08  cm  of  co rprene,  and  reconstructed  electrically  into  concentric 
rings  and  a  center  piston.  This  array  was  then  bonded  to  a  0.64  cm 
thickness  of  corprene.  The  face  of  the  array  was  electrically  tied  in 
common  while  each  ring  and  the  center  piston  were  wired  separately  on 
the  back  side  of  the  array.  The  array  was  mounted  in  an  aluminum 
housing  with  an  opening  of  21.59  cm.  Electrically  separate  portions  of 
the  array  were  transferred  through  feedthrus  to  an  air  space  in  the 
housing  behind  the  array  where  they  were  coupled  to  15.24  m  of 
Alpha  1120,  28  ga,  stranded,  twisted  pair,  shielded  cable.  The  13  sep¬ 
arate  cables  were  run  through  15  m  of  1.91  cm  i.d.  Nalgon  tubing.  The 
array  space  in  the  housing  was  covered  with  an  acoustically  transparent 
window  0.32  cm  thick,  made  from  Products  Research  and  Chemical 
Corporation  No.  1527  polyurethane  and  filled  with  castor  oil.  Several 
views  of  the  selectable  radius  circular  piston  projector  are  shown  in 
Fig.  4.2.  Each  of  the  cables  which  represent  the  sections  of  the  array 
are  terminated  on  barrier  strips.  A  selected  radius  piston  is  driven  by 
summing  appropriate  array  sections  in  parallel. 

b.  Transducer  Performance.  The  performance  of  the  selectable 
radius  circular  piston  projector  was  evaluated  at  Lake  Travis  Test 
Station  (LTTS)  using  calibration  equipment  common  to  acoustical  measure¬ 
ment  facilities.  The  accuracy  on  amplitude  measurements  is  assumed  to 
be  at  least  ±1  dB. 

The  transmit  voltage  response  for  the  transducer  was  measured 
at  low  level  using  a  NRL/USRD  E-8  standard  transducer  as  a  hydrophone  at 
a  distance  of  11.7  m.  The  transmit  voltage  response  is  shown  in 
Fig.  B.3  for  the  center  piston  along,  and  for  the  large  piston  formed 
by  summing  all  array  sections  in  parallel.  For  the  center  piston  alone, 
the  maximum  response  is  about  188  dB  re  1  yPa/Vrms  at  1  m,  which  occurs 
near  the  thickness  mode  resonance  frequency  at  600  kHz.  Far  below 
resonance  the  response  generally  increases  at  6  dB  per  octave  with 
increasing  frequency;  however,  the  variability  in  the  curve  indicates 
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coupling  of  other  modes  in  the  ceramic.  This  performance  is  common  in 
high  frequency  underwater  transducers.  The  response  for  the  large  pis¬ 
ton  has  similar  behavior  except  for  a  peak  at  about  280  kHz  and  a  mini¬ 
mum  at  about  380  kHz.  These  features  are  caused  by  a  radial  mode  that 
results  from  the  equal  spacing  of  the  concentric  rings  and  which  couples 
to  the  thickness  mode  of  the  ceramic. 

The  directional  properties  of  the  selectable  radius  circular 
piston  transducer  were  evaluated  with  a  0.5  msec  pulsed  carrier  at 
6 00  kHz.  Measurements  were  made  at  a  range  of  11.7  m,  which  is  1.16  Rq 
for  the  piston  with  the  largest  radius.  Beam  patterns  as  a  function  of 
radius  are  shown  in  Figs.  B.4-B.7.  The  sidelobe  structure  in  the  beam 
patterns  is  uniform  and  symmetric.  There  are  grating  lobes  in  each  of 
the  patterns  due  to  the  2.5X  center-to-center  spacing  between  the  con¬ 
centric  rings,  but  these  spacing  lobes  are  not  a  problem  in  terms  of 
parametric  operation.  The  beam  patterns  were  also  recorded  as  a  linear 
function  of  angle.  The  half-power  beamwidths  are  given  in  Table  IV. 2 
with  the  corresponding  theoretical  values  assuming  a  circular  piston 
set  in  an  infinite  rigid  baffle. 

The  impedance  of  the  transducer  and  cable  was  measured  at  low 
level  underwater  loaded  conditions.  These  measurements  are  provided 
in  Table  B.l. 

B.  Receive  Section 

A  block  diagram  of  the  signal  conditioning  portions  of  the 
receive  section  is  shown  in  Fig.  B.8.  The  ARL  239-1  hydrophone  is  used 
for  all  of  the  single  channel  measurements.  Its  output  is  coupled  to  a 
preamplifier  directly  when  measuring  primary  signal  characteristics  and 
indirectly  through  a  low  pass  filter  to  remove  remaining  signals  at  the 
primary  frequency  when  measuring  secondary  signal  characteristics.  The 
USRL-E-27-130P  is  used  only  for  cross  signal  analysis  in  the  secondary 
band  and  is  coupled  to  a  preamplifier  through  a  voltage  follower  and  a 
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FIGURE  B.4 

600  kHz  PRIMARY  FREQUENCY  BEAM  PATTERNS  MEASURED  AT  11.7  m  FOR 
THE  SELECTABLE  RADIUS  CIRCULAR  PISTON  PROJECTOR  (ARL  384-1) 
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FIGURE  B.5 

600  kH*  PRIMARY  FREQUENCY  BEAM  PATTERNS  MEASURED  AT  11.7  m  FOR 
THE  SELECTABLE  RADIUS  CIRCULAR  PISTON  PROJECTOR  (ARL  384-1) 
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FIGURE  B.7 

)  kHz  PRIMARY  FREQUENCY  BEAM  PATTERNS  MEASURED  AT  11.7  m  FOR 
THE  SELECTABLE  RADIUS  CIRCULAR  PISTON  PROJECTOR  (ARL  384-1) 
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TABLE  B.l 


ELECTRICAL  PROPERTIES  AT  LOW  LEVEL  OF  THE  SELECTABLE 
RADIUS  CIRCULAR  PISTON  PROJECTOR  (ARL  384-1)  AT  600  kHz 


Radius 

(cm) 

Magnitude  of 
Impedance 

03) 

Phase  Angle 
(deg) 

1.27 

35.8 

-22.5 

1.91 

21.0 

-18.0 

2.54 

16.7 

-  2.0 

3.18 

13.8 

9.0 

3.81 

12.9 

22.5 

4.45 

12.5 

35.6 

5.08 

11.7 

48.0 

5.72 

11.5 

53.4 

6.35 

11.4 

56.4 

6.99 

11.4 

59.7 

7.62 

11.4 

62.4 

8.26 

11.4 

64.5 

8.89 

11.4 

68.1 
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FIGURE  B.8 
RECEIVE  SYSTEMS 


low  pass  filter.  The  outputs  of  the  two  hydrophone  channels,  which  are 
separated  in  time  by  the  physical  orientation  of  the  hydrophones,  are 
summed  and  driven  through  a  140  ■  length  of  RG58  coaxial  cable.  The 
hydrophones  are  calibrated  as  receive  systems  at  the  cable  output.  The 
calibrated  output  is  coupled  to  a  final  set  of  preamplifiers  and  filters 
that  are  detailed  in  Fig.  B.8(b)  and  complete  the  signal  conditioning. 

The  linearity  of  the  high  frequency  receive  system  was 
evaluated  by  introducing  a  calibrated  tone  at  600  kHz  at  the  input  to 
the  preamplifier  in  Fig.  B.8(a).  The  results  were  recorded  at  the  sys¬ 
tem  output  in  Fig.  B.8(h'.  These  data  are  shown  in  Fig.  B.9,  which 

indicates  good  linrari*  ;  jw  a  system  output  of  about  +5  dB  re  1  V  . 

P 

1.  Hydrophone;; .  probe  hydrophones  were  used  to  measure  the 

sound  field  produces  the  broadband  parametric  array.  The  frequency 
range  of  interest  is  from  about  10  to  150  kHz  in  the  secondary  band  and 
around  600  kHz  in  r*e  primary  band.  Hydrophones  which  were  selected 
are  small  in  aperture  so  that  the  angle  subtended  by  the  hydrophone 
connection  is  small  and  the  ability  to  determine  subtle  changes  in  sig¬ 
nal  characteristics  as  a  function  of  angle  is  retained.  These  hydro¬ 
phones  are  also  operated  well  below  resonance  in  the  secondary  band  to 
provide  a  frequency  response  that  is  reasonably  uniform. 

a.  ARL  239-1  Hydrophone.  The  ARL  239-1  hydrophone  is  the 
principal  receiving  array  used  to  obtain  experimental  data;  it  was 
used  for  all  single  channel  measurements  in  both  the  primary  and  secon¬ 
dary  frequency  bands. 

The  ARL  239-1  hydrophone  consists  of  five  radially  polarized 
Channelite  5500  ceramic  cylinders.  Each  ceramic  cylinder  is  1.14  cm 
o.d.,  0.09  cm  in  height,  with  a  wall  thickness  of  0.13  cm.  The  ceramic 
cylinders  are  stacked  with  0.01  cm  spacers  of  onion  skin  paper  and  are 
electrically  coupled  in  parallel.  The  ceramic  stack  is  acoustically 
loaded  on  each  end  with  corprene  and  mounted  on  a  0.44  cm  diam  stainless 
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steel  shaft.  The  assembly  is  covered  with  Emerson  and  Cummings  CPC-16 
polyurethane.  The  electrical  connections  are  run  through  the  stainless 
steel  shaft  and  coupled  to  12.19  m  of  Belden  4548  twisted  pair,  shielded 
cable. 

The  freefield  voltage  sensitivity  at  the  hydrophone  output 
for  both  the  low  and  the  high  frequency  bands  is  shown  in  Figs.  B.10  and 
B.ll,  respectively.  The  resonance  frequency  for  the  radial  mode  of  the 
ceramic  is  about  1  MHz.  The  sensitivity  in  the  low  frequency  band 
varies  about  5  dB  from  10  to  70  kHz  and  has  a  broad  peak  centered  at 
about  100  kHz.  This  peak  results  from  a  coupled  secondary  mode  in  the 
ceramic.  The  sensitivity  in  the  high  frequency  band  varies  about  10  dB 
from  350  to  600  kHz  but  is  reasonably  flat  from  600  to  u00  kHz.  Low 
frequency  and  high  frequency  calibrated  receive  systems  were  constructed 
according  to  Fig.  B.8(a).  The  freefield  voltage  sensitivity  for  the  low 
frequency  and  high  frequency  receive  systems  is  given  in  Figs.  B.12  and 
B.13,  respectively.  Subtle  differences  in  frequency  response  are  noted 
between  the  data  for  the  low  frequency  hydrophone  and  receive  system. 
These  differences  result  from  the  interaction  of  the  low  pass  filter 
and  the  hydrophone.  The  hydrophone  output  is  buffered  by  an  amplifier 
with  a  high  input  impedance  in  the  high  frequency  receiver;  therefore 
there  is  little  change  in  the  frequency  response  of  the  high  frequency 
hydrophone  and  the  receiver  system. 

A  farfield  directivity  pattern  at  150  kHz  is  shown  in 
Fig.  B.14.  This  pattern  was  measured  at  a  range  of  11.7  m  in  the  plane 
normal  to  the  axis  of  the  cylinder.  This  pattern  shows  only  a  1  dB 
variation  over  360°,  which  is  excellent  for  a  probe  hydrophone. 

b.  USRD-E27-13QP.  The  USRD-E27-130P  was  used  only  for  cross 
signal  experimental  analysis.  This  hydrophone  is  a  U.  S.  Navy  standard 
and  is  composed  of  seven  PZT  discs,  3.17  mm  in  diameter  and  1.52  ram  in 
thickness,  which  are  spaced  inside  an  11.1  mm  diam  circle.  The  discs 
are  cemented  directly  to  a  butyl  acoustic  window  and  air  backed. 
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FIGURE  B.11 

FREEFIELD  VOLTAGE  SENSITIVITY  FOR  THE  ARL  239-1  HYDROPHONE 

(HIGH  FREQUENCY) 
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FIGURE  B.14 

FARFIELD  BEAM  PATTERN  AT  150  kHz  FOR  THE 
ARL  239-1  HYDROPHONE  MEASURED  AT  11.7  m 
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The  freefield  voltage  sensitivity  for  the  USRD-E-27-13QP  is 
shown  in  Fig.  B.15.  The  frequency  response  varies  a  maximum  of  3.5  dB 
from  70  to  400  kHz.  Before  use,  the  USRD-E-27-130P  low  frequency 
receive  system  was  always  calibrated  with  respect  to  the  output  of  the 
ARL  239-1  receive  system. 
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FIGURE  B.15 

FREEFIELD  VOLTAGE  SENSITIVITY,  E27-130P 


APPENDIX  C 

ANALYSIS  OF  UNGATED  SECONDARY  OUTPUTS 
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The  analysis  presented  in  this  appendix  is  directed  toward 
the  determination  of  secondary  signal  characteristics  for  the  broadband 
parametric  array  based  on  a  measurement  at  a  hydrophone  output.  In  par¬ 
ticular,  an  evaluation  of  the  artificial  removal  of  energy  which  is 
characteristic  of  the  impulse  response  of  the  hydrophone  is  addressed. 
This  removal  of  energy  is  accomplished  by  the  use  of  a  receive  signal 
gate  and  is  termed  artificial  since  it  is  simply  ignored  or  replaced  by 
zero  for  purposes  of  secondary  signal  analysis. 

Given  the  experimental  parameters  of  Table  IV. 1  and  a  circular 
piston  of  radius  8.89  cm,  the  secondary  outputs  as  a  function  of  angle 
off-axis  measured  at  23.6  m  with  a  primary  peak  sound  pressure  level  of 
195.2  dB  re  1  pPa  are  shown  in  Fig.  C.l.  Although  the  duration  of  the 
outputs  is  over  about  60  psec,  only  the  first  20  psec  is  representative 
of  the  secondary  signal  produced  by  the  broadband  parametric  array  with 
Gaussian  pulse  carrier  transmission.  The  remaining  40  psec  of  the  out¬ 
put  result  from  the  impulse  response  of  the  receiving  hydrophone.  It 
should  be  noted  that  a  number  of  hydrophones  were  implemented  and  in 
each  case  the  signal  portions  of  the  outputs  were  very  similar  and  the 
remainder  of  the  outputs  (many  of  which  were  an  additional  40-100  psec 
in  duration)  were  very  different.  In  going  from  on-axis  to  off-axis 
angles,  the  signal  portion  of  the  output  changes  from  a  dependence  on 
the  second  derivative  with  respect  to  time  of  the  square  of  the  Gaussian 
envelope  function  to  the  first  derivative.  During  this  transition,  the 
portion  of  the  output  that  results  from  the  impulse  response  of  the 
hydrophone  retains  its  longer  duration  features.  Since  the  gain  at 
each  angle  was  adjusted  to  normalize  the  waveform  presentations,  subtle 
features  in  the  output  are  obscured  by  noise  as  the  angle  off-axis  is 
increased. 


A  more  detailed  analysis  of  the  transient  hydrophone  outputs 
may  be  considered  in  the  frequency  domain.  In  Figs.  C.2  and  C.3  the 
normalized  power  density  spectra  for  the  measured  hydrophone  outputs 
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FIGURE  C.1 

OFF-AXIS  ACOUSTIC  WAVEFORMS  FOR  THE  BROADBAND  PARAMETRIC  ARRAY 

UNGATED,  f(t)  -  GAUSSIAN,  f,c  =  600  kHz,  SPL,  =  195.2  dB  re  1  /uPa 
a  =  8.89  cm,  RQ  =  10.07  m,  r  =  23.6  m 
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are  indicated.  Each  spectrum  has  250  Hz  resolution  and  was  computed 
from  the  incoherent  average  of  100  spectra.  The  energy  at  60  Hz  and 
the  low  level  oscillations  on  the  spectr. ,  especially  at  low  frequen¬ 
cies,  are  artifacts  of  the  analysis  and  not  a  part  of  the  data. 

On-axis,  the  power  spectrum  of  the  ungated  hydrophone  output 
roughly  extends  from  about  20  kHz  to  above  100  kHz  with  prominent  peaks 
at  about  55  and  81  kHz.  As  the  angle  off-axis  is  increased,  the  overall 
level  decreases  and  the  decrease  in  level  of  the  higher  frequency  com¬ 
ponents  is  progressively  greater  than  for  the  lower  frequency  compo¬ 
nents.  This  behavior  is  the  effect  of  the  off-axis  transfer  function 
for  the  broadband  parametric  array. 

The  hydrophone  output  may  be  represented  by 

0(t)  =  i(t)  *  h(t)  ,  (C.l) 

where 

0(t)  is  the  hydrophone  output, 
i(t)  is  the  hydrophone  input,  and 

h(t)  is  the  hydrophone  impulse  response  (0  for  t  <  0) . 

When  a  receive  signal  gate  is  applied  to  the  hydrophone  output,  the 
relationship  becomes 

0  (t)  =  0(t)  U(t)  -  (U(t)*h(t))  (C.2) 

g 

where 

0  (t)  is  the  gated  time  function,  and 

g 

U(t)  is  a  gate  function 
0  0  <  t  <  T 

1  otherwise 

If  the  output  is  a  transient  that  is  shorter  in  duration  than  the 
impulse  response,  then  the  application  of  the  gate  function  may  be 
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interpreted  as  restrictions  on  the  duration  of  the  impulse  response  and 
Eq.  (C.2)  may  be  written  as 


0(t)  =  i(t)  *  h'(t)  ,  (C.2) 

where 

h(t)  0  <  t  <  T 

h'  (t)  = 

0  otherwise 

From  Eq.  (C.2),  the  gate  function,  is  shown  to  reject  energy  that 
results  from  the  signal  input,  but  is  part  of  a  bandlimited  or  signal 
dependent  impulse  response.  The  effect  of  the  gate  is  to  reduce  the 
energy  at  the  higher  peaks  in  the  bandlimited  frequency  response  of  the 
hydrophone. 

Transforming  to  the  frequency  domain,  the  relationships  are 

written 


SF 


=  ff|i(t)J  SFh’(t) 


(C.3) 

(C.4) 


From  Eqs.  (C.3)  and  (C.4)  it  is  shown  that  an  additional  effect  of  the 
gate  function  is  to  diminish  resolution  by  averaging  or  smoothing  in 
the  frequency  domain. 


The  effects  of  the  receive  signal  gate  are  shown  in  Fig.  (C.2) 
where  the  normalized  power  density  spectrum  for  the  on-axis  gated  signal 
is  indicated.  This  spectrum  may  be  compared  directly  with  the  on-axis 
ungated  spectrum.  It  is  observed  that  the  levels  at  the  frequency  peaks 
near  20  kHz,  53  kHz,  and  81  kHz  have  been  reduced  and  that  averaging  or 
smoothing  in  the  frequency  domain  has  produced  a  smooth  quasi-Gaussian 
shaped  curve. 


i 
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